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A  USERS'  MANUAL 


ABSTRACT 

Two  network  analyzer  programs,  MATHNET  and  RISCA,  which  allow 
the  analyst  to  simulate  a  aeneral  class  of  network  representat ions  are 
described  and  evaluated  for  the  potential  user. 

Network  concepts,  program  listings,  and  oroqram  flow  charts 
are  included  for  both  proqrams  in  addition  to  detailed  aescrip+ion  of 
input  preparation  and  output  i nterpretat ion  for  a  hypothetical  example. 
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I .  INTRODUCTION 


1 . 1  Background . 

Although  having  much  wider  application,  the  primary  use  of 
network  analysis  has  been  in  the  planning  and  control  of  R&D  projects. 
The  use  of  network  analysis  techniques  for  this  purpose  had  its  origin 
with  PERT  (Program  Evaluation  and  Review  Technique)  in  1958  in  the 
Polaris  Program.  Since  that  time,  PERT  and  many  other  network  analysis 
techniques  have  gained  wide  acceptance  both  in  the  Department  of 
Defense  and  in  private  industry. 

In  the  last  two  years,  two  network  analyzer  programs, 

MATHNET*  and  RISCA**,  have  gained  acceptance  within  AMC.  Under  PROMAP 
70  the  responsibility  for  instruction  in  risk  analysis  was  assigned 
to  the  Army  Logistics  Management  Center  (ALMC)  at  Fort  Lee,  Virginia. 

In  order  to  expedite  the  program,  a  contract  was  let  to  MATHEMATICA 
by  the  Army  Research  Office  to  develop  this  course  of  instruction. 
MATHNET,  developed  by  MATHEMATICA  as  a  teaching  aid  for  this  course 
was  modified  by  ALMC.  They  call  their  version  of  the  program  RISCA. 

The  version  of  MATHNET  currently  being  used  at  this  agency 
is  not  identical  to  that  developed  by  Mathematica.  The  original 
Mathematica  program  contained  several  logic  inconsistencies  with 
respect  to  the  time  and  cost  values  generated  in  various  types  of 
nodes.  These  logic  errors,  corrected  by  analysts  at  Picatinny  Arsenal, 
are  not  in  the  version  of  MATHNET  discussed  in  this  manual.  The  RISCA 
version  discussed  herein  is  also  free  of  these  logic  inconsistencies. 

Even  though  the  two  programs  are  accepted  and  the  results 
used,  there  does  not  exist  an  adequately  documented  users !  manual 
wjiich  compares  both  programs  at  this  time.  Recognizing  the  utility 
of  network  analysis  techniques  in  the  Materiel  Acquisition  Decision 
Making  Process  and  the  need  for  such  a  manual  both  within  AMSAA  and 
AMC,  a  project  was  initiated  to  prepare  this  users1  manual. 

1 .2  Organization  of  the  Report. 

The  remainder  of  the  first  section  is  devoted  to  defining 
the  various  character i st i cs  of  a  network  and  the  types  of  networks 
that  can  be  modeled  using  RISCA  and  MATHNET.  In  Section  2,  the 
various  capabilities  of  the  two  network  analyzer  programs  are  described. 
In  Section  3,  a  description  of  the  construction  of  a  network  and  the 
method  for  inputting  data  to  the  program  is  provided.  The  output  of 
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both  programs  is  described  in  Section  4.  In  the  last  section  of  the 
manual,  the  program  capabilities  are  compared  and  contrasted  and  a 
recommendation  is  made  concerning  which  program  to  use. 

I .3  Network  Concepts. 

Before  the  two  network  analyzer  programs  are  described  and 
compared,  the  more-  basic  concepts  of  a  graph,  a  node,  an  arc,  a  network 
and  a  path  must  be  defined. 

Figure  l.l  is  an  example  of  a  graph.  The  circles  represent 
nodes,  and  the  lines  joining  the  nodes  are  called  arcs.  Hence,  a 
graph  is  a  collection  of  two  or  more  nodes  joined  by  arcs.  Any  arc 
can  be  characterized  by  the  pair  of  nodes  that  it  connects.  For 
example  (1,2)  characterizes  the  arc  connecting  nodes  I  and  2  in 
Figure  l.l. 

The  only  difference  between  a  graph  and  a  network  is  that 
the  arcs  have  some  type  of  flow  in  them  (see  Figure  1.2).  One  example 
of  a  system  that  can  be  represented  by  a  network  is  a  development 
test-program..  The  nodes*  in  a  development  test-program  represent  the 
initiation  or  completion  of  various  tests,  the  arcs**  represent  the 
actual  tests  being  conducted  and  the  flow  in  the  arcs  is  time  and/or 
cost  involved  in  testing. 

Finally,  a  path  is  defined  as  a  sequence  of  arcs  connecting 
two  nodes.  For  example,  the  following  sequence  of  arcs  form  paths 
between  nodes  I  and  4  in  Figure  1.2: 

PATH  (M)  (1,2),  (2,4);  PATH  (0)  (1,4) 

PATH  (N)  (1,3),  (3,4) 

I  .4  Types  of  Network  Representations. 

Given  the  preceding  concepts,  it  is  now  possible  to  describe 
the  different  types  of  network  representations  and  network  analysis 
techniques  for  analyzing  them.  The  differences  in  the  networks  result 
from  assumptions  made  concerning  the  events  and  the  flows  in  the 
activities  being  modeled  in  the  project.  As  mentioned  previously,  the 
arcs  represent  activites,  the  nodes  represent  events,  and  the  flow  in 
the  arcs  usually  represents  time  and/or  cost. 


* 

Nodes  generally  refer 
#* 

Arcs  general ly  refer 


to  events. 

to  activities  or  jobs. 
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Figure  1.1  Example  of  a  Graph 


Figure  1.2  Example  of  a  Network. 
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Three  types  of  network  models  will  be  discussed.  They  are: 

a.  Deterministic  event  and  activity  time/cost  network 

b.  Deterministic  event  and  probabilistic  activity 
time/cost  network 

c.  Probabilistic  event  and  activity  time/cost  network 

It  should  be  pointed  out  that  this  discussion  will  be  very  general  and 
only  the  major  attributes  and  assumptions  about  the  networks  being 
modeled  will  be  discussed. 

The  differences  in  the  types  of  network  representat ions  are 
most  easily  described  by  analyzing  one  example  under  the  varying 
assumptions  made  concerning  the  events  and  flows  in  the  activities. 
Changing  the  oil  in  a  car  is  the  example  that  will  be  used  throughout 
this  d i scuss ion . 

The  first  type  of  possible  network  representat ion  for  describ¬ 
ing  the  events  and  activities  involved  in  changing  the  oil  in  a  car  Is 
illustrated  in  Figure  1.3.  Here  there  are  f i ve  m i I estone  events  in 
changing  the  oil  in  a  car.  For  this  network  representation,  as  well 
as  all  others,  it  is  assumed  that  the  events  must  be  completed  in  a 
particular  sequence  In  order  to  complete  the  project.  Further,  it  is 
assumed  that  all  events  must  be  completed  and  the  completion  times  are 
known  with  certainty  (i.e.,  the  events  and  completion  times  are  assumed 
to  be  deterministic).  For  most  programs,  these  last  two  assumptions 
are  not  thought  to  be  very  realistic.  For  rarely  are  the  events  and/or 
activity  times  known  with  certainty.  Even  though  the  first  type  of 
network  representat ion  is  not  realistic  for  R&D  projects,  it  is  real¬ 
istic  in  the  construction  industry  where  tasks  for  a  project  are  known 
with  certainty.  Further,  these  tasks  are  repetitive  so  the  assumption 
of  deterministic  activity  times  is  more  realistic.  The  Critical  Path 
Methodology  (CPM)  is  the  name  given  to  the  network  analysis  technique 
developed  by  DuPont  in  1958  to  handle  deterministic  event  and  activity 
time  networks  (Type  I).  This  technique  was  initially  used  to  find  an 
efficient  method  for  planning  the  construction  of  a  new  facility. 

In  the  second  network  type,  the  assumption  of  deterministic 
activity  times  is  replaced  by  the  assumption  of  probabilistic  activity 
times.  This  means  that  the  activity  times  are  not  known  with  cer¬ 
tainty  (i.e.,  there  exists  some  distribution  of  activity  times).  For 
instance,  the  activity  time  for  BC  in  the  changing  the  oil  example 
could  vary  due  to  random  i nterruptions  such  as  having  to  pump  gas  for 
cars  as  they  arrive  at  the  station.  For  instance,  the  most  likely 
time  for  completing  this  activity  might  be  ten  minutes  and  the  best 
and  worst  times  might  be  five  and  twenty  minutes  respectively. 

Therefore,  if  this  activity  time  is  assumed  to  be  distributed  as  a 
triangular  distribution,  with  a  minimum,  most  likely  and  maximum 
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Figure  1.3.  Changing  the  Oil  in  a  Car. 
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Figure  1.5.  Changing  the  Oil  in  a  Car. 
(Probabilistic  Network) 


ac+ivi+y  times  being  5,  10,  and  20  minutes  respectively,  then 
Figure  1.4  represents  the  distribution  of  this  activity  time.  The 
network  analysis  technique  developed  to  analyze  this  type  of  network 
representat ion  is  PERT  (Program  Evaluation  and  Review  Technique)  which 
was  developed  in  1958  for  planning  and  controlling  the  development  of 
the  Polaris  missile.  Even  though  PERT  is  more  realistic  than  the  CPM, 
it  still  is  deficient  for  modeling  many  R&D  projects  since  the  events 
are  assumed  to  be  deterministic. 

In  the  third  network  type,  all  events  and  activities  are 
modeled  probabilistically.  This  type  i s  the  most  realistic  for 
modeling  R&D  programs.  To  illustrate  this  network  type,  assume  that 
it  may  not  always  be  possible,  for  example,  to  replace  the  oil  filter 
because  the  type  of  filter  required  may  not  be  in  stock.  Assume  that 
there  is  a  .98  probability  that  it  is  in  stock  and  a  .02  probability 
that  the  filter  is  not  in  stock.  The  probabilistic  nature  of  this 
network  type  is  illustrated  in  Figure  1.5.  In  addition,  when  each 
of  the  activity  times  is  modeled  probabilistically,  this  network  type 
is  beyond  the  scope  of  PERT  and  CPM  and  due  to  the  additional  complex¬ 
ity  introduced  by  probabilistic  events,  the  development  of  RISCA  and 
MATHNET  was  prompted. 

It  should  be  noted  that  both  MATHNET  and  RISCA  also  allow 
one  to  analyze  the  first  two  network  types. 


2.  RISCA  AND  MATHNET  CHARACTERISTICS  AND  CAPABILITIES 

RISCA  and  MATHNET  are  computer  programs  that  al low  one  to 
analyze  systems  that  can  be  represented  by  a  general  class  of  networks. 
Since  the  events  and  activity  times  and/or  costs  can  be  modeled 
probabilistically,  a  simulation  process  is  utilized.  The  output 
consists  of  a  frequency  of  occurrence  distribution  for  each  of  all 
possible  terminal  events  and  correspond i ng  time  and/or  cost  distribution 
for  each  terminal  event.  In  addition,  the  distribution  of  time  and/or 
cost  weighted  over  all  possible  terminal  events  is  estimated. 

Many  of  the  benefits  derived  from  analyzing  a  network  result 
from  the  analysis  and  thinking  that  is  required  in  the  construction 
of  the  network.  Consider,  for  example,  the  development  of  a  tank 
where  there  are  several  alternative  designs.  Describing  the  sequence 
of  events  for  alternate  development  programs  for  each  design  provides 
insight  into  the  types  of  problems  that  one  is  likely  to  encounter  in 
each  program. 

The  oil  changing  example  (Figure  1.5)  from  the  introduction 
is  utilized  to  demonstrate  the  characteri sties  and  capabilities  of 
MATHNET  and  RISCA. 
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Since  in  this  example  the  event,  "removing  the  oil  filter," 
is  the  only  uncertain  event,  the  probabilistic  event  network  can  be 
described  in  terms  of  two  deterministic-event  sub-networks.  One 
sub-network  represents  the  events  and  activities  involved  in  changing 
the  oil  and  the  oil  filter  and  the  other  represents  the  events  and 
activities  involved  in  changing  only  the  oil.  In  this  simple  example 
there  is  only  one  path  in  each  sub-network;  however,  in  more  realistic 
problems  there  will  almost  certainly  be  several  possible  paths  in  a 
sub-network. 

Monte  Carlo  procedures  are  used  to  determine  which  determin¬ 
istic  event  sub-network  will  be  followed  in  the  probabilistic 
event  network.  Each  of  the  sub-networks  have  a  terminal 
event  whose  completion  time  is  determined  by  Monte  Carloing  all  the 
activity  completion  time  distributions  within  the  sub-network.  All 
of  the  potential  paths  in  the  sub-network  are  then  investigated  using 
these  sample  activity  times.  If  the  sub-network  chosen  in  the  oil 
changing  example  involves  removal  of  the  oil  filter,  then  the  activity 
time  distributions  in  this  sub-network  would  be  randomly  sampled. 

These  sample  values  would  then  be  summed  to  estimate  the  sub-network 
completion  time.  In  addition,  the  cost  of  all  activities  in  the 
sub-network  would  be  sampled  and  summed  to  estimate  completion  costs. 
However,  in  this  example  costs  were  not  considered. 

The  preceding  procedure  is  repeated  many  times,  and  the 
sampling  distributions  of  terminal  events  and  time  and/or  cost  are 
constructed.  It  should  be  pointed  out  that  in  reality  the  determinis¬ 
tic  event  sub-network  is  chosen  and  the  correspond i ng  time  and  cost 
estimates  for  each  activity  on  this  network  are  accumulated  as  the 
network  is  simulated. 

For  this  example,  assume  that  300  iterations  have  been  run. 
Both  RISCA  and  MATHNET  would  provide  the  frequency  histograms  of  the 
percentage  of  times  each  terminal  event  was  selected  as  shown  in 
Figure  2.1  and  the  completion  time  distributions  shown  in  Figures  2.2, 
2.3,  and  2.4.  In  addition,  RISCA  would  provide  a  cumulative  distribu¬ 
tion  of  time  for  each  completion  time  distribution  in  Figures  2.2,  2.3, 
and  2.4.  Further  discussion,  interpretation  and  comparison  of  the  two 
programs1  output  are  deferred  to  the  Output  Section. 

If  cost  is  considered,  there  are  two  options  available: 

(I)  the  cost  can  be  estimated  independent  of  time  by  running  a  separate 
simulation  or  (2)  the  cost  can  be  estimated  as  a  linear  function  of 
time  in  the  same  simulation,  i.e.,  cost  =  (fixed  cost)+( var i ab I e  cost)x 
(time) . 


In  addition  to  the  insight  derived  in  structuring  and  simulat¬ 
ing  the  network,  this  type  of  analysis  can  provide  a  framework  for 
evaluating  and  consolidating  relevant  information  for  decision  making 
purposes.  If  one  is  trying  to  select  an  alternative  system,  this  type 
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Figure  24  Frequency  Histogram  for  Project  Completion  Times  Weighted  Over 
Terminal  Events  G  and  G'. 


of  analysis  can  provide  a  framework  for  evaluating  and  consolidating 
relevant  information  for  decision  making  purposes.  If  one  is  trying 
to  select  an  alternative  system,  this  type  of  analysis  can  provide  a 
reasonable  framework  for  comparing  time  and  cost  trade-offs  in  light 
of  possible  terminal  events  for  each  possible  alternative.  Consider, 
for  example,  the  development  of  a  tank  where  there  are  alternate  designs. 
These  alternate  designs  could  be  compared  on  the  chances  of  program 
failure  (i.e.,  having  to  cancel  the  development)  and  the  chances  of 
developing  the  tank  for  a  fixed  number  of  dollars  in  a  fixed  period 
of  time.  In  addition,  if  the  chances  of  failure  are  great  with  a 
preferred  alternative,  then  this  analysis  might  serve  as  the  basis  for 
restructur i ng  the  development  program  to  allow  for  parallel  develop¬ 
ment  of  two  designs  thereby  increasing  the  chances  of  successful 
development. 

In  addition,  this  type  of  analysis  could  provide  a  framework 
for  examining  the  progress  in  a  program  and  provide  inputs  to  determine 
whether  or  not  to  continue  the  program  or  to  assess  the  impact  of 
program  changes.  Once  the  network  has  been  constructed,  it  is  a  rel~ 
atively  simple  matter  to  modify  the  basic  structure  or  update  the  data. 
Therefore,  this  technique  could  provide  an  up-to-date  statement  of  the 
chances  of  meeting  program  time  and  cost  objectives. 

Using  these  programs  is  no  substitute  for  analysis.  Major 
effort  still  must  be  applied  in  trying  to  rea I i st ica I  I y  mode  I  the 
system  as  a  network.  These  programs  are  only  tools.  If  a  conscien¬ 
tious  job  of  modeling  the  system  has  been  accomplished,  then  the  output 
of  these  programs  is  useful  for  both  planning  and  decision  making. 

The  uses  of  the  output  in  a  decision  making  context  will  be  discussed 
in  more  detail  in  the  Output  Section. 


3.  RISCA  AND  MATHNET  PROGRAM  INPUT 
3.1  I ntroduct i on . 


This  section  describes  how  the  elements  of  a  network  are 
coded  and  inputted  into  RISCA  and  MATHNET. 

Included  is  a  description  of  the  arc  and  node  notation,  the 
program  input  formats,  an  example  problem  and  a  comparison  of  the  two 
input  decks  used  in  the  example. 

It  is  important  to  note  that  all  of  the  arc  and  node  notation 
and  the  input  formats  described  herein  are  the  same  for  both  RISCA  and 
MATHNET.  MATHNET,  for  time  sharing  mode,  contains  several  additional 
control  cards.  The  nature  of  these  cards  and  their  use  will  be  fully 
explained  when  the  input  decks  for  the  example  are  compared. 
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3.2  Arc  and  Node  Notation. 


As  explained  in  the  introduction,  the  two  classes  of  symbols 
used  in  describing  a  network  are  arcs  and  nodes.  An  arc  is  used  to 
connect  two  nodes  and  represents  an  activity.  A  node  is  used  to 
represent  a  decision  or  the  initiation  or  termination  of  an  activity. 


Each  arc  inputted  to  the  RISCA*  network  analyzer  program  is 
characterized  by: 

a.  An  arc  name 

b.  The  name  of  the  initial  node 

c.  The  name  of  the  terminal  node 

d.  A  distribution  of  completion  times 

e.  A  distribution  of  activity  costs  (as  a  linear  function 

of  time) 


f . 

The  fol lowing 

a . 

b. 

c . 

d. 


The  probability  of  successful  completion 

time  distribution  types  are  available: 

Norma  I 

Tr iangu I ar 

Un i form 

Constant 


Arcs  may  exist  in  one  of  the  following  states: 

a.  Idle  -  the  activity  associated  with  the  arc  has  not 
been  initiated 


b.  Initiated  -  the  activity  associated  with  the  arc  has 
been  initiated,  but  its  outcome  has  not  been  determined 

c.  Completed  -  the  activity  associated  with  the  arc  has 
been  successfully  completed 

d.  Unsuccessful  -  the  activity  associated  with  the  arc  has 
been  initiated  but  not  successfully  completed 


This  also  applies  to  MATHNET 


2! 


e.  Utilized  -  the  node  which  terminates  the  arc  has 
been  satisfied. 

The  conditions  which  determine  in  which  of  the  above  states 
an  arc  will  be  are  discussed  following  a  description  of  the  node  types. 

Nodes  represent  major  events  or  decisions.  Each  node  inputted 
to  the  program  is  characterized  by: 

a.  a  node  name 

b .  an  i nput  ru I e 

c .  an  output  ru I e 

In  addition,  time  and/or  cost  values  are  also  assigned  to 
the  nodes.  However,  these  values  vary  (except  for  constant  distribu¬ 
tions)  with  each  iteration.  The  time  and/or  cost  values  assigned  to 
the  nodes  are  those  derived  from  the  arcs  entering  them.  In  the  cases 
where  more  than  one  arc  enters  a  node,  the  exact  values  assigned  to 
the  node  are  dependent  upon  that  node's  input  rule.  These  rules  are 
described  in  detail  in  this  section. 

In  order  for  an  activity  to  be  completed,  it  must  satisfy 
the  conditions  of  the  node  that  it  enters.  For  example,  assume  that 
an  activity  is  "developing  a  model"  and  the  terminating  event  is 
"completion  of  the  model".  If  the  model  is  completed,  the  arc  is 
said  to  be  successf u I  and  the  node  "satisfied". 

The  rules  which  determine  whether  or  not  the  nodes  of  the 
network  are  satisfied  are  the  input  rules.  They  are: 

a.  And  Input  -All  input  arcs  must  be  successfully  completed. 
The  time  value  assigned  to  the  node  is  the  maximum  of  all  input  arcs. 

The  cost  value  is  the  sum  of  the  costs  generated  by  each  input  arc. 

b.  Or  Input  -  At  least  one  input  arc  must  be  successfully 
completed.  The  time  value  assigned  is  the  minimum  of  all  input  arcs 
successfully  completed.  The  cost  value  is  the  sum  of  the  costs 
generated  by  each  input  arc  based  upon  the  smallest  time  value. 

When  an  activity  has  satisfied  the  input  rule  of  a  node,  the 
nature  of  its  exit  from  the  node  is  dependent  upon  the  node  output 
rules.  There  are  two  such  rules: 

a.  All  -  all  output  arcs  are  initiated  simultaneously. 

b.  Probabilistic  -  Associated  with  each  output  arc  is  a 
probability.  The  sum  of  the  probabilities  exiting  one  node  must  equal 
one. 
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One  important  point  should  be  stressed  at  this  time.  The 
probabilities  assigned  by  the  analyst  to  each  arc  emanating  from  a 
node  are  done  so  on  node  data  cards.  This  probability  that  an  arc 
will  be  initiated  is  not  to  be  confused  with  the  probability  that  an  arc 
will  be  successfully  completed.  This  latter  information  is  entered  on 
an  arc  data  card. 

There  are  five  special  nodes  which  do  not  conform  to  the 
input  and  output  rules  outlined  above.  They  are: 

a.  Initial  Node  -  used  to  initiate  the  simulation.  It 
possesses  no  input  arcs 

b.  Terminal  Node  -  used  to  terminate  simulation.  It 
possesses  no  output  arcs 


c.  I/I  (one-to-one  bar)  Node  -  This  node  possesses  N  input 
arcs  and  N+l  output  arcs.  Associated  with  each  input  arc  is  a  unique 
output  arc.  The  extra  output  arc  is  a  default  arc.  The  time  value 
assigned  to  this  node  is  the  maximum  of  the  successfully  completed 

input  arcs.  The  cost  is  the  sum  of  the  cost  values  for  each  of  the 

input  arcs  which  are  completed  successfully.  This  time  and  cost 
assignment  is  true  for  the  remaining  node  types  in  this  discussion. 
This  node  is  satisfied  in  one  of  two  ways: 

(1)  by  having  one  successful  input  arc  (all  input  arcs 

must  be  i n i t ia ted )  which  will  resu It  in  one  output  arc  being  i n i t i a ted 

(2)  by  initiation  of  the  default  arc.  This  is  done  when 
all  of  the  input  arcs  have  been  initiated  but  none  successfully 
comp  I eted . 


The  one-to-one  bar  node  is  used  to  simulate  the  time,  cost 
and  probability  of  completion  for  competing  activities.  For  example, 
if  two  contractors  submit  the  time,  costs  and  probabilities  associated 
with  producing  a  certain  prototype,  this  data  can  be  applied  to  two 
arcs  entering  a  I / I  Node.  By  examining  the  correspond i ng  output  arcs 
and  the  default  arc,  the  developer  can  gain  valuable  insight  from 
success  and  failure  probabilities,  as  well  as  the  time  and  cost 
involved  with  each  alternative. 

d.  I/I  (one-to-one)  Node  -  This  node  possesses  N  input  arcs. 
Associated  with  each  input  arc  is  a  unique  output  arc.  The  node  will 
be  satisfied  by  one  complete  input  arc.  One  output  arc  will  be 
initiated,  namely  the  one  associated  with  the  input  arc  satisfying 
the  node.  I /I  Nodes  are  used  to  preserve  the  time  and  cost  factors 
for  two  or  more  activities  experiencing  a  common  event  (  the  event 
being  represented  by  the  I/I  Node).  For  example,  suppose  two 
contractors  are  developing  blueprints  for  a  house  and  the  cost  and 
time  for  each  contractor  differs.  If  the  probability  that  both  will 
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be  successfully  completed  is  1.0,  each  of  these  activities  (arcs)  may 
enter  a  common  node  (blueprint  completion).  As  the  arcs  exit  the  node, 
the  time  and  cost  of  each  activity  is  preserved. 

Unlike  the  I / I  node,  there  is  no  default  arc  to  be  initiated 
when  activities  are  unsuccessful.  Therefore  all  arcs  entering  the 
I/I  node  must  have  a  probability  of  completion  1.0.  A  probability  of 
less  than  1.0  may  cause  a  premature  termination  of  the  simulation 
resulting  in  no  output. 

e.  Preferred  Node  -  This  node  possesses  N  input  arcs  and 
Ntl  output  arcs.  Associated  with  each  input  arc  is  a  unique  output 
arc.  The  extra  output  arc  is  a  default  arc.  The  preferred  node 
differs  from  a  I/I  node  in  that  the  input/output  arc  pairs  are  ordered 
by  the  users  preference.  This  node  is  satisfied  in  two  ways: 

(1)  If  all  the  input  arcs  are  initiated  and  at  least  one 
arc  has  been  successfully  completed  the  output  arc  associated  with 
the  most  preferred  of  the  complete  input  arcs  will  be  initiated. 

(2)  If  all  the  input  arcs  are  in  the  unsuccessful  state  the 
default  output  is  initiated. 

It  is  important  to  note  that  this  node  will  never  be  satisfied  unless 
all  input  arcs  have  been  initiated  and  are  in  either  the  complete  or 
unsuccessful  state. 

This  node  is  used  to  simulate  the  time,  cost  and  probability 
of  completion  for  competing  activities  when  the  activities  can  be 
listed  by  order  of  preference. 

A  standard  set  of  symbols  is  used  to  identify  the  various 
node  types  and  input/output  rules  in  a  schematic  network  representa¬ 
tion.  A  listing  of  these  node  types  and  the  symbols  used  to  represent 
them  is  found  in  Table  3.1. 

3.3  Input  Format. 

The  format  used  to  enter  the  arc  and  node  data  into  the 
RISCA  network  analyzer  program  will  be  discussed. 

The  Data  deck  consists  of  3  blocks: 

Block  I  -  Header  Card 
Block  I  I  -  ARC  Cards 
Block  III  -  Node  Cards 
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TABLE  3.1  SYMBOLIC  REPRESENTATION  OF  INPUT/OUTPUT  RULES 
AND  SPECIAL  NODES 


INPUT  RULES 

AND 

OR 

INITIAL 

OUTPUT  RULES 

ALL 

PROB 

TERMINAL 

SPECIAL  NODES 

ONE-TO-ONE 

ONE-TO-ONE  BAR 

PREFERRED 


SYMBOL 


A 

N 

D. 


O 

R 


I 

N 

i 

T 


A 

L 

L 


p 

R 

O 

B 

1 

& 

1/1 


1/1 


PREF 
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A  description  of  each  follows: 

Block  I:  The  header  card  is  one  card  used  to  describe  the 
network  being  analyzed.  All  80  columns  may  be  used. 

Block  II:  The  information  required  for  each  arc  appears  on 
one  card  which  is  formatted  as  follows: 


1  nfo 

Co  1 umns 

Fo rmat 

a . 

Arc  Name 

1-4 

A4 

b. 

Input  Node 

5-8 

A4 

c . 

Output  Node 

9-12 

A4 

d. 

Di stri bution 

13 

II 

1  n 

column  13,  only  integers  1 ,  2,  3,  and  4 

may  be  entered . 

These  integers  are  used  to  indicate 
types: 

one  of  four  possible  distribution 

1 nteger 

Distribution  Type 

1 

Norma  1 

2 

T  r iangu 1 ar 

3 

Un i form 

4 

Constant 

e . 

Fi rst  time 
d  i  stri bution 
argument 

14-23 

FIO.O 

f . 

Second  time 
d i stri bution 
argument 

24-33 

FIO.O 

g- 

Third  time 
d i stri bution 
argument 

34-43 

FIO.O 

argument 

standard 

When  the  distribution  is 
is  the  mean  and  the  second 
deviation.  The  third  time 

normal,  the  first  time  distribution 
time  distribution  argument  is  the 
distribution  argument  is  left 

blank. 
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When  the  distribution  is  triangular,  the  first  time  distribu¬ 
tion  is  the  most  optimistic  value,  the  second  time  distribution 
argument  is  the  most  likely  value  and  the  third  time  distribution 
argument  is  the  most  pessimistic  value. 

When  the  distribution  is  uniform,  the  first  time  distribution 
argument  is  the  optimistic  value,  the  second  is  the  pessimistic  value 
and  the  third  is  blank. 

In  a  constant  distribution,  the  first  time  distribution 
argument  is  the  constant  value  and  the  second  and  third  arguments  are 
blank. 


h. 

Fixed  cost 

44-53 

FIO.O 

i  . 

Variable  Cost 

54-63 

FIO.O 

Coef f i c i ent 

j  • 

Probabi 1 i ty  of 

64-73 

FIO.O 

successf u I 
comp  I et ion 

Both  the  fixed  cost  and  variable  cost  coefficients  are 
components  of  a  linear  equation  which  describes  the  cost  as  a  function 
of  time.  This  equation  is  written: 

Cost  =  fixed  cost  *l*  variable  cost  coefficient  x  t 
where  t  is  time  in  appropriate  units. 

The  probability  of  successful  completion  refers  to  the 
probability  that  an  activity  (arc)  will  reach  its  desired  end  (node). 

The  arc  cards  which  comprise  Block  II  may  be  inputed  in  any 
order.  However,  it  is  desirable  to  enter  the  arcs  in  the  order  in 
which  they  appear  in  the  network  to  avoid  confusion. 

The  end  of  Block  II  is  indicated  by  a  card  following  the 
last  arc  card  in  Block  II  with  "RETU"  in  card  columns  1-4. 

Block  III:  The  information  required  for  each  node  is 


formatted  as 

fo 1  lows : 

1  nfo 

Co  1  iimns 

Format 

a . 

Node  name 

1-4 

A4 

b. 

1 nput  ru 1 e 

5 

II 
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There  are  6  possible  integers  (input  rules)  which  appear 
in  column  5: 


I nteger  V a  I ue 


2 

4 

5 

6 


I nput  Ru I e 
And 
Or 

Initial 

I/I 

T7T 


7 


Pref erred 


c.  Output  rule  6  II 

There  are  6  possible  integer  values  which  can  appear  in 
column  6.  Listed  below  is  each  integer  and  the  indicated  output  rule 


I nteger  Va I ue 

1 

2 

4 

5 

6 
7 


Output  Rule 
A I  I 
Prob 
Term i na I 
I/I 

I/I  Bar 
Preferred 


When  the  node  card  has  an  input/output  rule  I/I,  i / I  or 
PREF,  a  second  card  with  the  following  information  is  required: 


1  nf  o 

a . 

No.  of 

Output  Arcs 

1-2 

12 

b. 

Input  Arc 

Names 

■  3-6,  11-14, 

19-22,  etc. 

A4 

c. 

Output  Arc 

Names 

7-10,  15-18, 

23-26,  etc. 

A4 
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Each  input  arc  name  and  its  corresponding  output  arc  name 
will  appear  as  pairs  on  this  card.  The  first  input  arc  name  will 
appear  first  in  columns  3-6  followed  by  the  output  arc  name  it  initiates 
(columns  7-10).  The  second  input  and  output  arc  pairs  will  appear  in 
columns  11-14  and  15-18  respectively,  and  soon  until  all  arc  pairs 
are  I i sted . 

Since  the  default  output  arc  has  no  corresponding  input  arc 
name  ( I / (  and  PREF),  "ZZZZ"  must  be  entered  as  the  input  arc  name. 

The  above  information  card  must  immediately  follow  the  node  card  it 
describes. 


When  the  node  card  has  an  output  rule  Prob  (2),  a  second 
card  with  the  following  information  is  required: 


1  nfo 

Co  1 umns 

Format 

No.  of  Output 

Arcs 

1-2 

12 

Output  Arc 

3-6,  13-16, 

A4 

Names 

23-26,  etc. 

Probabi 1 i ties 

7-12,  17-22, 

F6.3 

associated  with 

27-32,  etc. 

On  this  card,  the  name  of  each  output  arc  and  the  probability 
of  its  exiting  appear  in  pairs  similar  to  the  input  and  output  arc 
pairs  of  the  previous  card.  For  example,  the  name  of  the  first  output 
arc  will  appear  in  columns  3-6  followed  by  the  probability  that  it 
will  exit  in  col umns  7-12  and  so  on  until  all  arc  name  and  probabi I ity 
pairs  have  been  listed.  This  card  must  immediately  follow  the  node 
card  it  describes. 

The  node  cards  in  Block  III  may  be  i nputed  in  any  order. 
However,  be  certain  that  in  the  case  of  I/I,  T7l,  Preferred  and 
Probability  nodes  that  each  node  card  is  succeeded  by  its  appropriate 
information  card.  Following  the  last  node  card  in  Block  III  is  a 
card  with  "RETU"  in  card  column  1-4  marking  the  end  of  the  biock. 

In  the  next  section  a  hypothetical  problem  will  be  described, 
structured  and  inputs  prepared  as  a  network  for  RISCA. 

3.4  Example  Problem. 

The  use  of  the  program  will  be  illustrated  via  an  example 
problem  containing  a  variety  of  arc  and  node  types. 
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A  worker  is  confronted  with  a  new  tardiness  policy 
established  by  his  employer.  It  has  been  decided  that  an  employee 
who  is  late  for  work  more  than  10  percent  of  the  next  500  working  days 
will  be  docked  accordingly.  The  amount  that  the  worker  will  be  docked 
will  be  the  cumulative  time  lost  due  to  lateness. 

This  announcement  has  prompted  the  worker  to  closely  examine 
the  various  routes  and  hazards  facing  him  each  morning  in  order  to 
evaluate  his  chances  of  his  being  docked. 

His  analysis  of  the  paths*  for  driving  to  work  revealed  the 
following  critical  areas: 

a.  A  fork  in  the  road  about  12  minutes  from  home.  It  has 
been  his  experience  that  the  shortest  route  is  only  possible  90  percent 
of  the  time  due  to  hazardous  road  conditions.  10  percent  of  the  time 

he  must  travel  a  considerably  longer  route. 

b.  Fuel  Problems.  The  worker  is  assured  of  not  running  out 
of  gas  if  he  takes  the  shorter  route.  However,  if  he  takes  the  longer, 
alternate  route,  there  is  an  80  percent  chance  that  he  will  have  to 
make  a  5  minute  fuel  stop. 

c.  Rough  Road.  Because  of  the  extreme  punishment  to  his 
tires  along  a  stretch  of  road  beyond  the  service  station,  the  worker 
feels  that  he  has  a  5  percent  chance  of  getting  a  flat  tire  that  would 
take  15  minutes  to  repair. 

d.  Rider  Stop,  Beyond  the  stretch  of  bumpy  road  is  the  home 
of  the  worker’s  friend  who  rides  to  work  with  him  every  day.  Historical 
data  reveals  that: 

(1)  If  the  shortest  route  is  taken,  and  no  fuel  stops  or 
tire  troubles  are  experienced,  his  rider  will  be  waiting  outside  for 
a  ride  100  percent  of  the  time. 

(2)  If  the  alternate  route  is  taken,  and  no  fuel  stops  or 
tire  troubles  are  experienced,  his  rider  will  be  wai1Tny  outside  80 
percent  of  the  time. 

(3)  If  he  takes  the  alternate  route,  stops  for  gas  and  has 
no  tire  trouble,  his  rider  will  be  waiting  for  him  60  percent  of  the 
time. 


Recall  that  a  path  is  a  sequence  of  activities  conned ing  two  events. 

In  this  instance  the  events  are  departing  for  work  and  arriving  at  work. 
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(4)  If  he  has  a  flat  tire,  regardless  of  previous  routes 
or  stops,  his  rider  will  only  be  waiting  outside  40  percent  of  the 
time. 


In  each  case  when  the  rider  is  not  outside  waiting,  he  must 
stop  and  check  to  see  whether  he  is  inside  or  has  obtained  another 
rider.  This  is  a  5  minute  stop. 

Figure  3.1  illustrates  the  alternative  routes  and  hazards. 

In  order  to  determine  which  node  and  arc  type  should  be 
used  to  simulate  the  route  in  a  network,  each  critical  area  should 
be  exam i ned . 

a.  Fork  in  the  road.  This  is  clearly  a  probabilistic  node. 

The  probability  that  the  arc  representing  the  shorter  route  will  emanate 
is  .90  and  the  probability  that  the  arc  representing  the  longer  route 
will  emanate  is  .10. 

b.  Gas  Station.  A  I / I  node  will  be  used  to  simulate  this 
event  although  a  probabilistic  node  could  also  be  used.  The  arc  enter¬ 
ing  this  node  represents  the  longer  route  path  and  a  probability  of 
completion  of  .8  will  be  assigned  to  it.  In  the  event  that  the  arc  is 
not  completed,  the  default  arc  representing  a  gas  stop  will  be  initiated. 

c.  Rough  Road.  Three  I / I  nodes  will  be  used  to  represent 

this  hazard.  Three  nodes  are  used  instead  of  one  in  order  to  account 

for  this  event  in  each  of  the  three  deterministic  subnetworks  developed 
thus  far.  Which  I / I  node  is  used  in  each  iteration  is  dependent  upon 
which  of  the  following  alternatives  precedes  it: 

( 1 )  Short  route, 

(2)  Long  route,  no  gas  stop, 

(3)  Long  route,  gas  stop. 

An  arc,  representing  each  of  these  alternatives,  will  be 
assigned  a  probability  of  completion  of  .95.  In  each  case,  failure 
to  complete  the  node  (get  through  the  bumpy  road  without  a  flat)  will 
initiate  a  default  arc.  All  three  default  arcs  will  enter  the  fix 
flat  node  which  will  have  an  OR  input  rule. 

d.  Rider  Stop.  This  event  will  be  represented  by  three 
I/I  nodes  and  an  AND/ALL  node.  The  AND/ALL  node  represents  picking 

up  the  rider  when  the  short  route  is  taken  and  there  is  no  tire  trouble. 
Recall  that  in  this  case,  the  rider  is  waiting  outside  100  percent  of 
the  time. 
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Figure  31  Graphical  Representation  of  the  Driving  to  Work  Problem. 


In  the  other  three  instances,  the  probability  that  the  rider 
will  be  waiting  is  less  than  1.0.  These  probabilities  are 


Event 

Probab i 1 i ty 

(1) 

Long  route,  no  gas  stop, 
no  tire  troub 1 e 

.8 

(2) 

Long  route,  gas  stop, 
no  tire  trouble 

.6 

(3) 

Tire  trouble 

.4 

Separate  I / I  nodes  will  be  used  in  each  of  these  instances 
with  all  three  default  arcs  entering  an  OR/ALL  node  representing  a 
stop  to  check  on  the  rider. 

e.  Arrive  at  Work.  Five  AND/TERM  nodes  will  be  used  to 
represent  the  worker’s  arrival.  Their  time  and  probability  values 
will  represent  the  time  and  probability  associated  with  5  separate 
alternatives  which  are: 


(  1  ) 

Short 

route,  no  stops. 

(2) 

long 

route,  no  stops. 

(3) 

Long 

route,  gas  stop. 

(4) 

Stop 

to  fix  flat. 

(5) 

Stop 

for  rider. 

Figure  3.2  illustrates  the  network  representat ion  of  "Driving 
to  Work"  using  the  RISCA  and  MATHNET  arc  and  node  construction  symbols. 

Table  3.2  lists  each  of  the  time  distribution  arguments,  the 
events  to  which  they  correspond,  the  distribution  arguments,  the 
events  to  which  they  correspond,  the  distribution  type  involved  and 
the  arc  names  assigned  to  each  of  these  time  consuming  activities. 

Now  that  the  network  has  been  graphically  represented,  the 
arcs  and  nodes  have  been  labeled  and  the  time  distribution  argument 
has  been  determined,  the  input  deck  can  be  prepared. 

As  explained  earlier,  there  are  some  differences  in  MATHNET 
and  RISCA  with  respect  to  input  deck  content. 

The  contents  of  the  RISCA  deck  are  the  three  blocks  outlined 
previously.  That  is,  the  first  card  is  BLOCK  I  consisting  of  one  title 
card  describing  the  network  to  be  analyzed.  This  is  followed  by 


33 


ND10  ND16 


CM 

U 

oc 

< 


CO 

u 

QL 

< 


ClOCQcQ 


<zo 


a 


CM 

o 

o 

Z 


o 

< 

o 

CL 

Z 

CL 

o 


u 

CL 

< 


O 

O 

z 


<_l  -I 


-Z-K- 


UJ 

S 

o 

x 


34 


Figure  3.2  Network  Representation  of  the  Driving  to  Work  Problem 


TABLE  3.2  CHARACTERISTICS  OF  ARCS  USED  IN  EXAMPLE  PROBLEM 


Distribution 

Arc 

Time  Distribution 

Arguments 

Activity 

Type 

Name 

Optimistic 

Expected 

Pessimist 

From  Home  to 
Junction 

Triangular 

Arcl 

10.0 

12.0 

14.0 

Junction  to 

Gas  Station 
(short  route) 

Triangular 

Arc  2 

5.0 

6.0 

7.0 

Junction  to 

Gas  Station 
(long  route) 

Triangular 

Arc  3 

11.0 

13.0 

15.0 

Gas  Station 
Stop 

Constant 

Arc  8 

5.0 

Gas  Station 

Arc6 

3.0 

4.0 

5.0 

to  Bumpy 

Arc7 

3.0 

4.0 

5.0 

Road 

Triangular 

Arc9 

3.0 

4.0 

5.0 

Stop  to 

Change 

Flat  Tire 

Constant 

AR1 5 

15.0 

Bumpy  Road  To 

AR1 0 

4.0 

5.0 

6.0 

Rider 1 s 

AR1 2 

4.0 

5.0 

6.0 

House 

Triangular 

AR14 

4.0 

5.0 

6.0 

AR16 

4.0 

5.0 

6.0 

Stop  to 

Check  Rider 

Constant 

AR23 

5.0 

Rider's  House 

AR17 

8.0 

9.0 

10.0 

To  Work 

Triangular 

AR18 

8.0 

9.0 

10.0 

AR20 

8.0 

9.0 

10.0 

AR22 

8.0 

9.0 

10.0 

AR24 

8.0 

9.0 

10.0 
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BLOCK  II,  the  arc  cards,  the  last  card  being  a  "RETU"  and  indicating 
termination  of  the  block.  Then  comes  BLOCK  III,  the  node  cards,  the 
last  card  again  being  a  "RETU"  card  indicating  the  end  of  the  block. 

The  version  of  RISCA  described  in  this  manual  is  designed 
only  to  run  in  a  batch  mode  while  MATHNET  is  designed  to  run  in  either 
a  batch  or  time  shared  mode.  It  should  be  emphasized  that  adapting 
RISCA  for  a  time  sharing  environment  is  not  a  difficult  task.  However, 
it  would  probably  not  be  possible  to  run  these  versions  of  MATHNET  and 
RISCA  in  any  batch  processing  mode  without  first  adapting  the  program 
to  the  particular  system.  Similarly,  it  would  probably  not  be  possible 
to  run  this  version  of  MATHNET  in  any  time  sharing  system  without 
adapting  it  to  the  system.  Because  the  MATHNET  program  can  be  run 
in  either  mode,  the  input  deck  is  designed  to  give  operators  workinq 
in  time-share  the  capability  of  introducing  either  each  section  of 
the  input  deck  with  a  card  describing  the  section  or  a  card  describing 
an  operation  to  be  performed.  When  the  program  reads  the  introduction 
card,  it  is  directed  to  that  portion  of  the  program  which  performs 
the  indicated  operation. 

The  first  card  in  the  MATHNET  input  deck  is  a  card  with  either 
a  0  or  a  I  punched  in  column  I  (format  II).  The  0  card  indicates  that 
the  problem  will  be  run  in  a  time  share  environment,  and  a  I  indicates 
batch  mode. 


For  the  purposes  of  comparing  RISCA  and  MATHNET  inputs  in  this 
manual,  only  the  input  decks  for  a  batch  mode  environment  will  be 
utilized  (See  Figure  3.3). 


The  information  cards  mentioned  above  are  cards  with  an  integer 
punched  in  column  I  (format  II).  The  following  is  a  list  of  each  integer 
which  may  be  used,  and  the  information  card  it  indicates: 


I nteger 

1 

2 

3 

4 

5 

6 
9 


I nd i cates 
Node  cards  foil ow 
Arc  cards  fol low 
Net  iteration  number  is  set 
Net  is  scanned 
Net  is  run 

Identification  card  follows 

Run  is  to  be  completed  (Ends  Session) 
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Note  that  the  number  of  iterations  is  controlled  by  the  operator 
in  MATHNET.  In  RISCA,  the  number  of  iterations  is  fixed  at  500.  It  is 
conceivable  that  only  500  iterations  could  be  restrictive  in  a  network 
having  a  large  number  of  terminal  nodes.  However,  as  was  mentioned  above, 
the  program  can  be  easily  expanded  to  allow  for  more  iterations.  The 
card  following  the  information  card  with  a  3  in  column  I  is  a  card  with 
an  integer  value  in  columns  1-5  (format  15)  specifying  the  number  of 
iterations  (from  I  to  1000). 

The  information  cards  with  integers  4,  5  and  9  punched  in 
column  I  do  not  introduce  input  cards,  but  they  indicate  to  the  program 
that  one  of  three  operations  are  to  be  performed.  These  operations  are: 

a.  Scan  the  net  -  A  scan  of  the  net  means  that  al I  arc  and 
node  characteristics  are  printed  out  in  a  tabular  format  (see  Section  4). 

b.  Run  the  net  -  The  net  simulation  is  done. 

c.  Run  is  Completed  -  All  output  from  the  simulation  is 
printed  out.  (See  Section  4). 

The  arc,  node  and  identification  cards  are  formated  in  the 
same  manner  as  RISCA.  However,  the  order  in  which  they  are  inputted 
is  not  as  restrictive  as  RISCA.  The  data  cards  must  be  entered  before 
the  operation  cards,  but  the  order  in  which  the  data  cards  are  entered 
is  not  important.  The  operation  cards,  however,  must  be  entered  as 
fo I  I ows : 

a.  Scan  the  net 

b.  Run  the  net 

c.  Complete  the  run. 

Figure  3.3  provides  a  comparison  of  the  RISCA  and  MATHNET 
input  decks  used  in  the  example  problem.  This  figure  illustrates  the 
fact  that  with  the  exception  of  fhe  MATHNET  control  cards  discussed 
above,  the  two  decks  are  identical. 

Since  the  RISCA  input  deck  has  already  been  discussed  in 
detail,  and  the  fact  that  its  MATHNET  counterparts  are  identical  has 
been  established,  only  an  explanation  of  the  MATHNET  confrol  cards  in 
the  example  problem  input  deck  remains  to  be  covered. 

The  first  control  card  indicates  that  the  problem  will  be 
run  in  a  batch  mode  environment  (I  in  column  I).  The  second  card 
indicates  that  the  problem  identification  card  follows  (6  in  column  I). 

The  third  control  card  (3  in  column  I)  indicates  that  the  number  of 
iterations  follows.  The  next  card  sets  the  number  of  iterations  at 
500.  The  following  card  indicates  that  the  arc  cards  follow  (2  in 
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column  I).  The  next  control  card  appears  at  the  end  of  the  arc  card 
block  and  indicates  that  the  node  cards  follow  (I  in  column  I). 

The  last  3  cards  in  the  MATHNET  input  deck  control  the 
running  of  the  simulation.  The  first  card  initiates  the  scanning  of 
the  net  (4  in  column  I).  The  second  card  indicates  the  simulation 
is  to  be  run  (5  in  column  I).  The  last  card  indicates  that  the  run 
is  to  be  completed  and  the  results  printed  out  (9  in  column  I). 

It  is  important  to  note  that  the  blanks  which  appear  in  the 
listing  of  the  RISCA  input  deck  in  figure  3.3  were  included  for  ease 
of  comparison  only  and  are  not  required. 

It  was  pointed  out  above  that  there  are  certain  restrictions 
as  to  the  number  of  iterations  permissable  in  these  versions  of  MATHNET 
and  RISCA.  In  addition  to  these  iteration  restrictions,  there  are 
certain  other  inherent  restrictions  which  deserve  mention.  For  a 
complete  listing,  see  Table  3.3.  These  limits  may  be  increased  (through 
minor  program  modifications)  to  the  upper  bound  of  the  computer  memory 
core . 


4.  RISCA  AND  MATHNET  PROGRAM  OUTPUT 
4.1  I ntroduction. 


In  this  section,  the  output  from  MATHNET  and  RISCA  for  the 
goi ng-to-work  example  will  be  described  and  compared.  In  addition, 
the  potential  usesof  the  output  for  decision  making  purposes  will  be 
discussed  both  in  and  out  of  the  context  of  this  example. 

4.2  MATHNET  Output. 

The  output  of  any  MATHNET  run  consists  of  a  detailed  input 
listing,  frequency  histograms  for  completion  time  and  cost  for  each 
of  the  terminal  nodes,  frequency  histograms  of  completion  times  and 
cost  weighted  over  all  terminal  nodes,  and  frequency  histograms  of  the 
percent  of  time  each  terminal  node  was  selected  in  the  simulation. 

A  complete  output  listing  for  MATHNET  for  the  example  is  provided  in 
Appendix  IV. 

The  first  four  pages  of  output  consists  of  a  printout  of 
the  dialogue  between  the  user  and  the  program,  and  a  detailed  description 
of  the  network  representation  being  simulated.  (See  the  Input  Section). 
This  dialogue  is  really  only  applicable  to  the  terminal  user.  It 
provides  him  with  a  detailed  set  of  instructions  for  inputting  his  network 
representation.  This  example  was  run  in  a  batch  mode,  therefore  the 
first  data  card  inputted  into  the  program  had  to  have  a  I  in  the  first 
column  of  the  first  card.  On  the  two  pages  following  the  dialogue  there 
is  a  detailed  listing  of  the  arc  and  node  input  data.  This  is  particularly 
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TABLE  3.3  Restrictions  to  MATHNET  and  RISCA 


MATHNET 

RISCA 

Maximum 

number 

of 

nodes 

100 

100 

Maximum 

number 

of 

arcs 

500 

100 

Maximum  number  of 
into,  or  out  of  a 
node 

arcs 

single 

10 

10 

Maximum 

nodes 

number 

of 

initial 

10 

10 

Maximum 

nodes 

number 

of 

terminal 

30 

30 

Maximum 

number 

of 

iterations 

1000 

500 
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useful  information  for  checking  the  network.  The  format  of  the  arc  and 
node  data  will  not  be  discussed  since  this  is  thesame  format  that  was 
discussed  in  the  Input  section. 

Next,  for  each  of  the  terminal  nodes  (ND  16,  17,  ND  18,  ND  19, 

ND  20)  a  frequency  histogram*  of  the  completion  time  is  provided.  In 
this  example,  costs  were  not  considered,  but  if  they  were  then  there 
would  be  a  frequency  histogram  of  cost  for  each  terminal  node.  In 
addition,  each  terminal  node  represents  the  event  "arriving  at  work" 
however,  in  each  instance  the  path  (i.e.  sequence  of  events)  is 
different.  For  example,  the  first  histogram  is  for  the  completion 
times  for  ND  16.  ND  16  is  the  short  route  to  work  with  no  time  delay 
associated  with  hazards,  stopping  for  gas,  or  picking  up  a  rider.  The 
horizontal  axis  is  the  frequency  of  occurrence  and  the  vertical  axis 
is  time  in  minutes.**  Looking  at  this  graph,  the  probability  of  arriv¬ 
ing  at  work  within  35.6  and  35.9  minutes  given  this  route  is  taken  is 
.087.  The  interpretation  of  all  the  remaining  frequency  histograms 
is  the  same  with  two  exceptions.  One  is  that  if  the  terminal  cost 
frequency  histograms  were  given  then  vertical  axis  would  be  cost. 

The  other  is  that  the  final  pair  of  frequency  histograms  for  completion 
times  and  cost  (where  considered)  is  weighted  over  all  possible 
terminal  nodes  (routes). 

The  final  graph  provided  is  the  frequency  plot  of  the  percentage 
of  times  each  terminal  node  was  selected  in  the  simulation.  The 
horizontal  axis  again  represents  frequency  of  occurrence  and  the 
vertical  axis  lists  the  possible  terminal  nodes.  Looking  at  the  graph 
in  Appendix  IV,  the  probability  of  arriving  at  work  having  taken  route 
ND  17  is  .046,  or  ND  17  was  chosen  in  4.6  percent  of  the  simulations. 


4.3  RISCA  Output. 

The  output  of  any  RISCA  run  consists  of  a  detailed  input 
listing,  frequency  and  cumulative  frequency  histograms  for  completion 
times  and  cost  for  each  of  the  terminal  nodes,  frequency  and  cumulative 
frequency  histograms  of  completion  times  and  cost  weighted  over  all 
terminal  nodes,  and  a  frequency  histogram  of  the  percent  of  times  each 
terminal  node  was  selected  in  the  simulation.  A  complete  listing  of 
RISCA  for  this  example  is  provided  in  Appendix  V. 


# 

On  each  time  and  cost  frequency  histogram  the  following  statistics 
are  provided:  mean,  variance,  median,  and  mode. 

** 

The  time  and  cost  units  are  selected  by  the  analyst  and  depend 
upon  the  problem. 
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The  first  page  of  output  consists  of  a  detailed  listing  of 
the  arc  and  node  input  data.  This  is  particularly  useful  information 
for  checking  the  network.  Once  again,  the  format  of  the  arc  and  node 
data  will  not  be  discussed  since  this  is  the  same  format  that  was 
discussed  in  the  Input  Section. 

Next,  for  each  of  the  terminal  nodes  (ND  16,  ND  17,  ND  18, 

ND  19,  ND  20)  a  frequency  and  cumulative  frequency  histogram*  of  the 
completion  times  and  cost  is  provided.  It  is  emphasized  again  that 
costs  were  not  considered.  In  addition,  each  terminal  node  represents 
the  event  arriving  at  work,  but  in  each  instance  the  path  Ci.e.,  sequence 
of  events)  is  different.  Since  the  i nterpretation  of  the  frequency 
histograms  is  the  same  as  for  MATHNET,  only  the  cumulative  frequency 
histograms  will  be  discussed.  For  example,  the  first  histogram  is  for 
the  completion  times  for  ND  16.  It  is  followed  by  the  cumulative 
frequency  histograms  for  completion  times  for  ND  16.  The  vertical  axis 
is  time**  to  completion,  and  the  horizontal  axis  is  the  probability 
that  the  true  time  is  less  than  or  equal  to  the  time  on  the  vertical 
axis.  For  instance,  the  probability  of  arriving  at  work  in  less  than 
35.647  minutes  given  that  this  path  was  taken  is  .401.  The  interpreta¬ 
tion  of  all  of  the  remaining  frequency  and  cumulative  frequency  histo¬ 
grams  is  the  same  with  two  exceptions.  One  is  that  for  cost  frequency 
histograms,  the  vertical  axis  is  cost**.  The  other  is  that  the  final 
pair  of  frequency  and  cumulative  frequency  histograms  for  completion 
times  and  cost  is  weighted  over  all  possible  terminal  nodes  (routes). 

The  final  graph  provided  is  the  frequency  plot  of  the  percent 
of  time  each  terminal  node  was  selected  in  the  simulation.  The 
i nterpretat ion  of  this  histogram  is  identical  to  the  MATHNET  interpreta¬ 
tion  of  the  same  histogram. 


4 . 4  MATHNET  and  RISCA  Output  Comparison. 


There  are  only  two  differences  in  the  output  of  these  programs. 
One  is  that  RISCA  provides  cumulative  frequency  histograms  for  time  and 
cost,  for  terminal  nodes,  and  MATHNET  does  not.  Clearly,  the  cumulative 
frequency  histogram  is  requ i red  information  in  most  applications  as  for 
example,  in  estimating  the  probability  of  meeting  program  time  and  cost 
goals.  Consequently,  there  is  an  advantage  in  using  this  version  of 
RISCA  over  this  version  of  MATHNET  since  it  el im inales  the  need  for 
manually  generating  cumulative  frequency  histograms. 

* 

On  each  time  and  cost  histogram  and  following  stalistics  are  provided: 
mean,  variance,  and  standard  deviation. 

** 

The  time  and  cost  units  are  selected  by  the  analyst  and  depend  upon 
the  problem. 
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The  other  difference  is  in  the  output  statistics  computed 
for  each  time  and  cost  distribution.  In  MATHNET,  the  mean,  variance, 
median  and  mode  are  computed  while  in  RISCA  the  mean,  variance  and 
standard  deviation  are  computed.  On  first  glance  it  may  appear  that 
one  is  given  more  information  with  the  MATHNET  statistics,  but  this 
is  not  the  case.  Since,  if  one  is  given  the  frequency  and  cumulative 
frequency  histogram  for  any  terminal  node,  determining  the  mode  and 
median  is  a  minor  operation. 

4 . 5  Use  of  the  Output. 

While  there  are  more  ways  to  use  the  output  of  a  network 
analysis  than  will  be  discussed  in  this  section,  the  two  decision 
problems  that  are  discussed  probably  represent  the  more  popular 
applications.  Before  describing  these  problems  and  how  to  use  the 
output  tor  decision  making  purposes  several  general  comments  should 
be  reiterated. 

First,  the  greatest  benefit  to  be  derived  in  using  a 
network  analyzer  program  comes  from  the  effort  that  is  put  into 
model i ng  the  project  or  system  as  a  network.  Using  this  type  of 
tool  forces  one  to  examine  all  possible  events  and  the  interaction 
of  these  events  in  the  program.  Further,  representing  a  system  as 
a  network  allows  one  to  handle  all  of  the  relevant  decision  information 
in  a  systematic  and  composite  fashion  and  to  evaluate  the  impact  of 
interactions  that  would  be  otherwise  impossible. 

Next,  no  matter  what  the  decision  is,  there  will  always  be 
a  need  for  an  overall  framework  for  consolidating  all  the  information. 
Network  analyzer  programs  such  as  MATHNET  and  RISCA  provide  such  a 
f ramework. 


Finally,  once  the  network  has  been  structured,  future 
modifications  to  the  network  are  a  relatively  simple  matter.  In  many 
instances,  this  type  of  analysis  should  not  be  a  one  time  effort.  The 
network  representat ion  and  time  and  cost  estimates  should  be  modified 
on  a  periodic  basis  because  as  time  passes,  more  information  is  gained 
and  the  initial  network  may  no  longer  realistically  represent  the 
system. 


In  the  remainder  of  this  section,  the  use  of  network  analysis 
for  two  general  decision  problems  and  the  going  to  work  example  will 
be  discussed. 

The  first  decision  problem  either  is  one  where  the  program* 
status  must  be  evaluated  periodically  to  determine  whether  or  not  to 

* - 

Program  refers  to  a  development  program  for  a  system. 
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continue  or  one  where  the  imoact  of  a  program  change  must  be  evaluated. 
The  discussion  that  follows  will  only  be  in  terms  of  the  decision  to 
continue  for  brevity  purposes,  although  network  analysis  could  be  used 
to  help  in  a  program  change  decision.  A  good  example  of  a  program 
would  be  a  weapon  system  development.  Clearly,  the  decision  to 
continue  depends  on  the  chances  of  successfully  developing  the  system 
within  time  and  budget  constraints. 

A  network  representat ion  for  a  development  program  would 
probably  consist  of  two  sets  of  outcomes.  One  set  would  constitute 
successful  development,  and  the  other  set  would  constitute  failure. 

In  this  case,  the  chances  of  successful  development  equals  the  sum 
of  the  percent  of  times  each  successful  terminal  was  selected  in  the 
simulation.  For  any  successful  outcome  the  probability  of  developing 
the  system  within  time  and  cost  constraints  may  then  be  evaluated. 

This  information  is  taken  directly  from  the  RISCA  output  (i.e.  frequency 
histogram  of  the  percent  of  times  the  terminal  nodes  were  selected  and 
the  appropriate  terminal  node  cumulative  time  and  cost  histograms).  In 
addition,  without  a  great  deal  of  effort,  the  frequency  and  cumulative 
frequency  histograms  for  time  and  cost  weighted  over  all  successful 
outcomes  are  estimated.  This  is  done  by  either  modifying  the  network 
representat ion  so  that  all  successful  terminal  events  feed  into  one 
terminal  event  (successful  development),  or  it  may  be  constructed  f  rom 
the  detailed  frequency  histograms  of  the  successful  outcomes. 

This  type  of  information  should  give  the  decision  maker  a 
foundation  upon  which  to  decide  whether  to  continue  with  the  program 
in  light  of  the  risks.  Of  course  this  doesnft  relieve  the  burden  of 
decision  making,  but  it  should  enable  the  decision  maker  to  make  a 
more  informed  decision.  Further,  as  the  program  continues,  the  network 
could  be  updated  periodically  to  evaluate  the  program  status  for  any 
future  program  decisions. 

The  other  decision  problem  is  one  in  which  a  choice  between 
alternative  systems  for  meeting  a  particular  set  of  requ i rements*  must 
be  made.  For  this  type  of  problem,  a  network  could  be  structured  for 
each  alternative.**  Each  of  these  networks  could  then  be  simulated 
and  the  chances  of  successfully  developing  each  alternative  system  and 
the  corresponding  development  costs  and  completion  times  could  be 
estimated.  Given  this  information,  the  analyst  has  several  relative 
measures  for  comparing  the  alternatives,  and  the  decision  maker  should 

This  assumes  that  "meeting  a  particular  set  of  reau i rements . . . " 
generates  equal  and  required  effectiveness. 

** 

It  is  assumed  that  a  system  i  s  to  be  developed. 
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be  able  to  systematically  consider  time  and  cost  trade-offs  in  light  of 
the  risks. 


Even  though  these  two  hypothetical  decision  problems  are 
over-simplified,  they  are  representat i ve  of  many  of  the  decisions 
that  must  be  made  within  AMC.  In  both  of  these  decision  problems 
there  is  a  great  need  for  an  overall  framework  for  the  evaluation.  All 
too  frequently  this  framework  seems  to  be  lacking.  Network  analysis 
provides  an  overall  framework  for  synthesizing  a  large  portion  of  the 
decision  information  in  a  systematic  fashion  that  takes  into  account 
the  known  uncertainties.  It  is  not  meant  to  imply  that  the  network 
analysis  is  the  answer  to  all  of  the  decision  makers  problems,  but 
reasonable  applications  are  a  step  in  the  right  direction. 

Returning  now  to  the  going  to  work  example:  How  can  this 
analysis  help  the  worker  in  achieving  his  objective?  First,  it  is 
assumed  that  the  objective  is  not  to  be  docked  for  tardiness.  Since 
in  all  instances  the  worker  arrives  at  work,  the  distribution  of 
interest  is  the  arrival  time  distribution  weighted  over  all  possible 
outcomes.  Assuming  the  workman  continues  to  leave  his  house  at  0705 
every  morning  and  work  starts  at  0745,  the  probability  of  his  being 
late  is  .157  or  1.0  minus  the  probability  of  arriving  at  work  in  less 
than  or  equal  to  40  minutes  (.843)  (See  the  cumulative  frequency 
histogram  for  arrival  times  weighted  over  all  outcomes  in  Appendix  V.) 


If  the  probability  of  being  late  on  any  given  day  that  the 
worker  departs  at  0705  is  .157,  the  risk  of  being  late  over  10  percent 
of  the  next  500  working  days  is  the  probability  of  being  late  51  or 
more  times.  Let  t  be  the  random  variable  representing  the  number  of 
late  arrivals  and  p  the  probability  of  being  late  on  a  given  day,  then 


P[5I  <  l  <  500]  = 


500 

l 

1=  51 


5°° 

'  l  } 


l 

P 


(  I  -p) 


500-1 


is  the  risk  of  being  late  in  excess  of  10  percent  of  the  time. 

If  the  worker  leaves  at  0705,  the  above  calculation  results 
in  a  risk  of  . 99(p= . I  57) . 

Using  this  method  it  is  possible  to  derive  a  r i sk-of-bei ng- 
docked  profile  for  the  worker  as  a  function  of  departure  time  (See 
F i gure  4.1). 

What  should  the  worker  do?  The  answer  is  not  clear  cut. 

It  depends  upon  the  worker’s  assessment  of  the  value  of  additional 
rest  versus  the  potential  financial  impact  of  being  docked.  However, 
given  the  risk  profile,  the  worker  is  in  a  position  to  consider  the 
trade-offs.  For  example,  if  the  worker  cannot  possibly  afford  to  be 
docked,  he  can  select  a  departure  time  from  the  risk  profile  where  the 
risk  is  approximately  zero  (i.e.,  leave  at  0654  or  earlier). 
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Thus,  it  is  seen  in  this  simpie  example  that,  through  the 
modeling  of  a  network  using  RISCA  or  MATHNET,  the  analyst  gains 
significant  insight  into  the  risks  inherent  in  a  system  composed  of 
various  alfernatives.  This  insight  into  the  risks  provides  the 
information  for  making  decisions  concerning  which  alternative  to 
choose,  or  decisions  involving  the  various  cost  and/or  time  parameters 
(in  this  example  departure  time  is  the  deciding  factor). 

Of  course,  much  more  complex  systems  containing  a  myriad  of 
activities,  events  and  alternatives  may  be  modeled  using  these  modeling 
tools.  The  scope  and  complexity  of  each  network  is  only  constrained 
by  the  number  of  arcs  and  nodes  that  can  be  inputted. 


5.  SUMMARY 

MATHNET  and  RISCA  are  programs  for  simulating  a  class  of 
networks  where  both  the  events  and  activities  can  be  modeled  probabil¬ 
istically.  These  programs  should  provide  the  analyst  with  a  tool  for 
evaluating  the  status  of  existing  programs,  for  evaluating  the  impact 
of  proposed  program  changes  on  the  total  program,  for  deciding  whether 
to  continue  the  program,  and  for  deciding  between  alternative  systems. 
These  decisions  may  not  be  mutually  exclusive. 

Since  RISCA  is  a  modification  of  MATHNET,  there  are  very  few 
differences  in  the  input  and  output  of  the  programs,  and  there  are  no 
differences  in  the  method  of  simulating  the  networks.  As  seen  in  the 
Output  Section,  the  additional  features  of  this  version  of  RISCA's 
output  to  this  version  of  MATHNET's  output  are  not  really  significant. 

Throughout  this  report  R&D  applications  have  been  emphasized, 
however,  these  techniques  obviously  have  broader  applications. 


47 


Next  page  is  blank. 


BIBLIOGRAPHY 


A  Course  of  Instruction  in  Risk  Analysis,  ALM-33I5-H,  United  States 
Army  Logistics  Management  Center,  Fort  Lee,  Virginia,  UNCLASSIFIED. 

MATHNET  Mathemati ca 1 s  Network  Analyzer  Program,  Preliminary  Edition, 
August  1970,  Mathematica,  Princeton,  New  Jersey,  UNCLASSIFIED. 


49 


Next  page  is  blank. 


APPENDIX  I 


GENERAL  FLOW  CHART  AND  LISTING  OF  MATHNET 
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This  appendix  includes  a  general  processing  flow  chart  for 
MATHNET,  but  it  does  not  include  a  description  or  flow  chart  of  the 
individual  subroutines.  A  detailed  description  and  flow  chart  of  the 
individual  subroutines  are  provided  in  Appendix  iii.  in  addition  to 
the  flow  chart,  a  complete  program  listing  is  provided.  This  version 
of  MATHNET  is  designed  to  be  run  in  either  a  batch  or  time  sharing 
mode. 


MATHNET  is  written  in  FORTRAN  IV.  Even  though  FORTRAN  IV 
is  considered  to  be  standard  language,  adapting  the  program  for  a 
particular  computer  will  probably  require  minor  program  modifications. 
These  modifications  generally  result  from  pecularities  of  the  given 
system. 


There  is  one  definite  modification  that  the  user  must  make 
before  MATHNET  can  be  run  on  his  computer.  To  run  the  simulation  it 
is  necessary  to  generate  uniform  random  numbers.  Within  the  sub¬ 
routine  RANDU  these  random  numbers  are  generated  by  calling  the 
library  subroutine  peculiar  to  the  particular  machine.  in  this 
listing  the  library  subroutine  RANSET  is  called  to  generate  the 
uniform  random  numbers.  Therefore,  the  user  must  either  call  the 
appropriate  uniform  random  number  generator  for  his  machine  or  the 
analyst  must  write  his  own  uniform  random  number  generating  routine 
within  the  RANDU  subroutine.  The  latter  was  done  by  ALMC  in  their 
modification  of  MATHNET.  The  interested  reader  is  referred  to 
Appendix  II  where  the  ALMC  Program  is  listed. 
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U  O  n  l)  U  O  O  U  L)  U  O  I  3  u  u  u  u  u  u  u  U  U  O  l)  O  U  U  U  (J  u  o  oo  .  oo  u  u 


PROGRAM  MATHW  (  I NPUT , TAPF 5= I NPUT , OUTPUT , TAPE6* OUTPUT , 


♦  PUNCH, TAPE7=PUNCH,TAPEl,TAPE2,TAPE3,TAPE4  ) 

** *** ************ * *** *  *  *** ** **********************************  ****t*hj 000 20 

********************************************************* *********M AT00010 

MAT000  30 

MATHNET  WAS  DEVELOPED  BY  MATHEMATICA  INC.  FOR  THE  ARMY 
RESEARCH  OFFICE  UNDER  CONTRACT  NUMBER  DAHC  04  70  COC25. 

ALTHOUGH  THIS  PROGRAM  HAS  BEEN  EXTENSIVELY  TESTED  WITH  A  LARGE 
NUMBER  OF  N tTW ORKS , EXHAUST  I VE  TESTING  OF  ALL  POSSIBLE  CLASSES 
OF  NETWORKS  IS  COMPUTATIONALLY  I MP RACT I C AL . M ATHEMA T I CA  THEREFORE 
MAKES  NO  GUARANTEES  CONCERNING  THE  ACCURACY  OF  THE  OUTPUT 
BY  A  RUN  OF  THE  MATHNET  PROGRAM 

IT  WOULD  BE  APPRECIATED  IF  ANY  NETWORKS  WHICH  PRODUCE  ERRONEOUS 
OUTPUT , uR  ANY  OTHER  SYSTEM  ERRORS,  WOULD  BE  REPORTED  TOO 
MR.  STEPHEN  ROBINSON 
MATHEMATICA 
ONE  PALMER  SQUARE- 
PRINCET  N ,  NEW  JFRSEY  08540 


SUGGESTIONS  CONCERNING  ADDITIONAL  OUTPUT 
CHARACTERISTICS, OR  NODE  CAPAB I L I T I ES , ARE 


FEATURES, ARC 
ALSO  WELCOMED. 


MATOOO  40 
MAT0005C 
MAT00060 
MATOOO  70 
MAT00080 
GENERATED MATO 00 90 
MATOOIOO 
MATOOI 10 
MAT001 20 
MATOOI 30 
MAT00I40 
MAT00150 
MA  TOO l 60 
MATOOI 70 
MAT00180 
MATOOI 90 
MAT0020C 

******* ********************************** **************** *********MAT002 10 
**************************************************** *******«******MAT002 20 

MAT00230 

MATHNET  -  VERSION  1  MODIFICATION  LEVEL  0  -  AUGUST  1,  1970  MAT00240 

MAT0C250 

**************************** ******************************** ******MAT002 60 
****************************************************************** MAT0027C 

MAT002  80 
MAT00290 
MAT00300 
MAT003 10 
MAT00320 
MAT00330 
MAT00340 
MAT00350 
MAT00360 
MAT00370 
MAT00380 
MAT00390 
MAT  00400 
MAT00410 
MAT00420 
MAT  004  30 
MAT00440 
MAT00450 
MAT  00460 
MAT  00470 
MAT00480 
MAT00490 
MAT00500 
MATOOslO 
MAT00520 
MAT00530 
MAT00540 
LISMAT00550 
MAT00551 
MAT  00560 
MAT00570 
MAT00580 


THIS  IS  THE  MAIN  ROUTINE 
COMMON/RANC/I XX 
COMMON/ INDEV/  INDV 

INITI AL IZE  RANDU 
IXXX=65 j39 

CALL  RANDU! IXXX, IXYY.RVAL) 
I XX= I XYY 
WR IT  E ( 6 , 20 ) 


FIRST 
IF  WE 


I  WILL 
ARE  IN 


READ  IN  A  RECORD  FROM 
BATCH  MODE  OR  NOT 


DATA  SET  5  TO  DETERMINE 


YOU  ARE  RUNNING  THIS  FROM  A  TERMINAL) 

IF  RUNNING  BATCH  YOU  SHOULD) 


ENTER  A  1, 


20  FORMAT! IX,  39HIF 
WR ITE ( 6  ,21 ) 

21  FORMAT!  IX,  45HPLEASE 
WR I TE ! 6 , 22 ) 

22  FORMAT!  IX,  28HHAVE  ENTERED  A  CARD  WITH  A  0) 

WR I T  E ! 6 , 23 ) 

23  FORMAT!  IX,  1 2HF0RMAT  IS 
READ! 5, 24)  INDV 

24  FORMAT!  ID 

:  NOW  TEST  MODE,  IF  BATCH 

IF! INDV.EQ.O)  GO  TO  7 
WR ITE ! 6 , 30 ) 

30  FORMAT!  IX,  53HY0U  ARE  NOW 
CT  ) 

WR ITE ! 6 ,31 ) 

31  FORMAT!  IX,  35HSELECT  THE  MODE  YOU  WISH  TO  GO  INTO) 
WR I TE ! 6 , 32 ) 


II  ) 


SKIP  EXPLANATIONS  FROM  NOW  ON 


IN  MONITOR  MODE  ,  FROM  THE  FOLLOWING 
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o  o  o 


32 

FORMAT * IX,  16H1 

ENTER  NODES) 

MAT00593 

WRITE*6,33) 

MAT00600 

33 

FORMAT * IX,  15H2 

ENTER  ARCS) 

MAT00610 

WR I T E  *  6 , 34 ) 

MAT00O20 

34 

FORMAT* LX,  25H3 

SET  ITERATION  NUMBER) 

MAT  006  30 

WR  1 1 E  *  6 , 35 ) 

MAT0C64 

35 

F  ORM A T *  IX,  24H4 

SCAN  THE  NET  SO  FAR) 

MAT0C65  ) 

WRITE! 6,36) 

MAT00660 

36 

FORMAT * IX,  L2H5 

RUN  NET) 

MAT00670 

WR I T E  *  6 , 37 ) 

MAT00fc>83 

37 

TCRMAT *  IX,  25H6 

ENTER  RUN  IDENTIFIER) 

MATO 06  90 

WRITE*6,38) 

MAT00700 

38 

FORMAT*  IX,  L6H9 

END  SESSION) 

MAT0G7  10 

GO  TO  7 

MAT00720 

77 

WRITE*  6,50) 

MAT0073 J 

WRITE*6,78) 

MAT00740 

78 

F  ORMA  T (  1 X  ,  33HYOU  HAVE  RETURNED  TO  MONITOR  MODE) 

MA  T 007  50 

WRITE(6»79) 

MAT00760 

79 

F0RMAT(1X,  65HSFLECT  THE  MODE  YOU  WISH  TO  GO  INTO  AS 

I NO I CA  T  FD  BEFMAT0077C 

CORE  ) 

MAT00771 

7 

WRITE  *  6,50) 

MAT00790 

50 

FORMAT*  2H  * ) 

MATOOBOO 

R  E AD  *  5 ,  1  )  MOOT YP 

MATOOR  10 

1 

FORMAT*  11  ) 

MAT00820 

IF ( MOOT YP«  EO* 1 ) 

CALLNODIN 

MAT00830 

I F  *  MOOT YP  *  E  0#  2 ) 

CALL  ARCIN 

MAT00840 

I F  *  MOOT  YP • ECU  3 ) 

CALL  REPSET 

MAT00850 

I F  *  MOOT  YP  *  E  0#  4 ) 

CALL  SCAN 

MAT00860 

IF*M0DTYP.E0#5) 

CALL  RUNSYS 

MAT008  70 

I F  *  MOOT YP#  EQ#  6 ) 

CALL  IDIN 

MAT008  80 

I F  *  MODT  YP #  E  C#  9  ) 

CALL  EXIT 

MAT00890 

GO  TO  77 

MAT0090U 

END 

MAT  009 10 

subroutine  arcin 

THIS  ROUTINE  WILL  SERVE  TO  READ  IN  ARCS  AND  ASSOCIATED  DATA 

COMMON/ INDEV/INDV 
COMMON/ PAR A /NNODE. N ARC 

COMMON  /ARC2/  TI ME ( 500 ) , I S TA T ( 500 ) , PR  OB < 500 ) 

COMMON/ ARC!/  ARC  I  5  00) , INODE!  500 ) , ONODE I  5  00 ) , I T I  MET <  500 ) , 
CTARGK  5  00),TARG2(500)  ,  TARG3I  500)  ,COSTC(500  )  .COST  VI  500 
COMMON/ NODI/  NODE (100)  , 1  ARC (  100 . 1 0 )  .  OARC ( 100 . 10 ) » PPJ A( 100, I  0 ) , 
Cl  ARC  I ( ICO), OARC I ( 100) , I  TYPE ( 100) ,OTYPE ( 100 ) , MNI ND< 100 ) , 

CTIMEN( ICO) 

REAL  NODE 
INTEGER  ONODE 
INTEGER  OARC I , OARC 
DATA  RETC/4HRE  TU/ 

C 

C 

I F  (  INDV.EQ.'i)  GO  TO  50 
WR  ITE  (  6, 60) 

60  FORMAT (  IX ,  29HY0U  ARE  NOW  IN  ENTER  ARC  MODE) 

150  WR ITE ( 6 ,6 1 ) 

61  FORMAT ( 2H  * ) 


MAT00920 
MAT  00930 
MAT00940 
MAT00950 
MAT00960 
MAT00990 
MAT01000 
MAT00970 
MAT00971 
MATO  10  10 
MATO  1J11 
MATO  1 0  1 2 
MAT01040 
MATO  1080 
MAT  0 10  50 
MATO  1060 
MATO  1070 
MAT  01090 
MATO  1 1 00 
MAT011 10 
MATO  1 1 20 
MATO  1 1 30 
MAT01140 


56 


WR ITE ( 6, 70)  PAT  0 1 1 50 

70  FORMAT! IX,  48HENTER  ARC  NAME, INPUT  NODE  NAME, OUTPUT  NODE  NAME,/  MAT01160 
C  IX,  57HTIME  DISTRIBUTION  TYPE, TIME  DISTRIBUTION  ARGUMENTS  1 , 2, 3 , MATO  1 1 6  1 
C  /IX,  71 MCONST  ANT  COST  CO EFF I C I  ENT ,COEFF I C I  ENT  OF  TIME  TERM  IN  MAT01162 
C  COST  TERM, //IX,  4QH PROBABILITY  OF  SUCCESSFUL  ARC  COMPLET I ONMATO 1 1 63 

C.)  MATO  1 1 64 


WR ITE ( 6,71 ) 

MATO  1200 

7  1 

FuRMAT ( IX ,  23HF0RMAT  IS  3A4, I 1 , 6F 1 0, 0 ) 

MATO  1 2  10 

WR ITE 1 6, 6  1 ) 

MAT  012  20 

WRITE(6,72) 

MATO  1 2  30 

77 

WR I TE l 6 , 61  ) 

MATO  1 240 

. 

72 

FORMAT!  IX,  38HT0  RETURN  TO  MONITOR  MODE  ENTER  RETU  ) 

MATO  1 2  50 

so 

READ! 5,1)  ANAME , A  INODE , AON  ODE , I D I  ST , P 1 , P2 , P3  ,C1 ,C2  ,D 

MAT  012  60 

1  FORMAT! 3A4, 11, 5F10.0.F9.0) 

I F( ANAME.EO.RETC  )  RETURN 

MATO  1 2  80 

NARC=NARC+1 

MAT  0 1 2  90 

ARC( NARC) =ANAME 

MA  TO  1 3  00 

I T I  ME  T ( NARC )  =  IDI$T 

MATO  13  10 

fARGl(NARC)=Pl 

MAT  0 1 3  20 

T  ARG2 ( N  ARC ) =P2 

MA  TO  1 3  30 

T  ARG3 ( N  ARC ) =P 3 

MATO  1 3  40 

CO  STC ( N  ARC ) =C 1 

M  ATO 13  50 

CGS  T  V ( NARC ) =  C2 

MATO  1 360 
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O  o 


PROB I NARC )  =  D 


MATO  L  3 70 


9 

10 


1  L 


1  2 


20 

C 


30 


40 


IF(NNODE.EO.O)  GO  TO  10 
PO  9  I  -  1 , NN  DE 
I S  AVE  =  I 

ISAVE  SIMPLY  KEEPS  VALUE  OF  I  TO  USE  OUT  OF  00  LOOP 

IFI AIN0DE.EU*N0DEII) )  GO  TO  II 

CONTINUf- 

NNODE  =  NN“»DE  +  I 
NODEINN  JOE) =  A I  NODE 
INODE(NARC) =  NN  00  E 
OARC ( NNODE t 1)=NARC 
L ARC  I (NNODE )  =  L 


GO  TO  12 

INODE l N ARC ) - I  SAVE 

OARC I  I ISAVE ) - 9 ARC  I  (  ISAVE )  ♦  1 

L M  =  0 ARC  I (  ISAVE ) 

OARC (  ISAVEt  LM)*NARC 

CONTINUE 

DO  20  I  = 1 » NNODE 


I S  AVE  =  I 

IF(AONODE.EU#NODEI  I  )  )  GO  TO  3'* 
CONTINUE 

WE  HAVE  A  NEW  NODE 
NN0DE=NNCDE+1 
NODE(NNCDE)=AONODE 
ONODE l N ARC ) -NNODE 
IARC(NN°DE, 1)=NARC 
I  ARC  I  INNODE  )  =  1 
GO  TO  4 

ONODE ( N  ARC ) = I S AV  E 
I  ARC  I (  I  SAVE )  =  I  ARC  I  (  I  SAVE)* l 
L  M= I  ARC  I  (  ISAVE) 

I  ARC  I  I  SAVE  t  LM ) =N  ARC 
CONTINUE 
GO  TO  77 
END 


MATO  1  3  80 
M  A  T  0 1 3  90 
MAT0140 J 
M  A  T  0  1 4  1 0 
MAT  0 14  20 
MATO  1 4  30 
MAT0144G 
MATO  14  50 
MAT  014  60 
MATO  1470 
M ATO 1480 
MAT01490 
MATO  1  SCO 
MAT  C 15  1C 
MATO  1520 
MATO  1 5  30 
MATO  154C* 
MATO  1 5  50 
M A  TO  1 5  60 
MATO  15  70 
M  ATO 1 5  80 
MATO  1 5  90 
MAT0160C 
MATO  1 6  1C 
MAT01620 
MA TO  1 6  30 
MATC164u 
MAT01650 
MATO  1 6  60 
MAT  0 1670 
MA  TO  1 6  80 
MAT  016  90 
MATO  17  00 
MATO  1 7 10 
MATO  1720 


SUBROUTINE  NODIN 

THIS  ROUTINE  WILL  READ  IN  NODE  NAMES  AND  DATA 
COMMON/ INDEV/  INDV 

COMMON /NODI/  N OD E ( 1 00 ) , I AR C (  100 , 10 ) t OARC 1 100  , 10) ,PPOAI 100,10), 

Cl  ARC  I ( 190), OARC I ( 100), I  TYPE (  100) ,0 TYPE (100) ,MNIND( 100) , 

CT  I  MEN  I  100) 

REAL  NODE 

DIMENSION  DDUM ( 10) 

COMMON/PARA/  NNODE, NARC 

COMMON/ TERN  N /NQDN ( 30 ) ,NODl ,TIMEZ( 1000)  ,COSTZ (10^0) ,NODEZ I  1000) , 
CNCOUNTI 30 ) 

COMMON  /ARC2/  T I  ME ( 500 ) , I S TA T ( 500 ) , PROB ( 5 00 ) 

COMMON/ ARC  1 /ARC ( 500) ♦ I NODF (500) , ONODE ( 500 ) ♦ I TI MET  I  500 ) ♦ TARGl ( 500 ) 
CTARG2( SCO ), TARG3 ( 500 ) » COSTC ( 500) ,C0STV(500) 

D  I  MENS  ION  ZNAM  I  1 0 )  ,  ZPROB 1 1 0)  ,  IZNAM(IO)  ,YI  (10)  ,  Y  'J  ( 1 0  ) 

INTEGER  OARC , 0 ARC  I 
INTEGER  OTYPE 

COMMON/ TERN  1/ I  NOD  I ( 10  )  ,  IN9DT 
DATA  CZZ/4HZZZZ/ 

DATA  RETC/4HRETU/ 


MATO  1 7  3C 
MATO  1 7  40 
MATu  17  50 
MATO  1 7  60 
MAT  017  70 
MATO  1 77 1 
MA  TO  1 7  72 
MATO 1 8  00 
MATO  1 8  10 
MATO  1 8  20 
M  ATO 1 830 
MAT  018  31 
MATO  1 8  50 
r  MATO  1 8 60 
MATO  1 8  6 1 
MATO  1 8  80 
M A  TO  1 900 
MATO  1930 
MATO  1920 
MATO  18  90 
MAT  019  10 


58 


r* 

C 

r 

1 

50 

51 

52 

53 

540 

54 

55 

56 

57 

58 

59 

20 

30 

40 

100 

101 

C 

222 

C 

C 

60 

r 

300 

301 


MA  TO 1 940 


TEST  IF  BATCH,  IF  SO 

SKIP  GOBBELDYGOOK 

MATO  1950 
MAT01963 

I F ( INDV.EQ.O)  GO  TO 

20 

MATO  1970 

TERMINAL  MODE  PRINT 

INSTRUCTIONS 

HAT01980 

WR  ITE ( 6 , l ) 

MA  TO  1 990 

FORMAT!  IX,  26HY0U  ARE  IN  ENTER  NODE  MODE) 

MAT02000 

WR I T E  <  6  »  50 ) 

MAT02010 

FORMAT!  IX,  38HENTER 

NODE  NAME, INPUT  RULE, OUTPUT  RULfc) 

MAT02020 

WR I T  E  <  6 , 5  1 ) 

MAT02033 

FORMAT!  IX,  18HF0RMAT 

IS  A4,I 1, I I  ) 

MAT  02  040 

WRITE(6,52) 

MAT020  50 

FORMAT (  IX,  37HINPUT 

AND  OUTPUT  RULES  ARE  AS  FOLLOWS) 

MAT02060 

WR I T  E ( 6 , 53 ) 

M  A  T  0  2  j  7  0 

FORMAT!  10X,  11HRULE 

NUMBER, IOX,  LOHINPUT  RULE , 5 X ,  11H0UTPUT 

RULE)  MAT02080 

WR I T  E ( 6 , 540 ) 

MAT02090 

FORMAT! 9X , 1 3 ( 1H, ),8X 

,12(1H.)  ,13X,13(1H.) ) 

MATO  2 1 00 

WRITE<6,54) 

MAT021  10 

FORMAT!  1 5  X  t  1H1 , 1 8X , 

3HAND ,2  OX ,  SHALL  FIRE) 

MAT  0  2120 

WRITE(6,55) 

MAT02130 

FORMAT!  1 5  X , 1H2«18X« 

2H0R, 21 X ,  IOHPROB.  FIRE) 

MAT92140 

WR  I T  E ( 6 , 56 ) 

MATO  2 1 50 

FORMAT! 15X, 1H4, 18X, 

7HINITIAL, 16X,  8HTERMINAL) 

MAT02160 

WRITE(6,57) 

MAT  0  2170 

FORMAT!  15  X , 1H5 , 18X, 

3H 1 / 1 , 2  OX ,  3HI/1) 

MATO  2  1  80 

WRITE(6,58) 

MATO  2190 

FORMAT!  15X, 1H6, 18X, 

7H1/I  BAR.16X,  7H1/1  BAR) 

MAT02200 

WRITE(6,59) 

MAT  02210 

FORMAT!  15  X , 1H7, 18X, 

9HPREFERRED,14X,  9HPRE  FE  RREO ) 

MAT02220 

WRITE!6,30) 

MAT02230 

FORMAT! 2H  * ) 

MAT02240 

READ! 5,40)  ADUM, IDUM 

, I0DUM 

MAT02250 

FORMAT!  A4,  I  1,  ID 
I F ( ADUM.EQ.RETC)  RETURN 
DO  100  l -1, NN9DE 
NDUM  = I 

I F ( NODE ( I  ) • FQ. ADUM )  GO  TO  10  1 
CONT INUF 
NNODE*NNODE ♦ 1 
NDUM=NNODE 
ITYPE(NDUM)*IDUM 
0  TYP  E ( NDUM )  =  I  ODUM 
I F ( ITYPE!NDUM).E0.4)  GO  TO  90 
I T YP E  =  4  MEANS  INITIAL  NODE 
I F ( ITYPE(NDUM) .E0.6)  GO  TO  400 
I  TYPE-6  MEANS  WE  HAVE  A  1/1  NODE 
I F ( ITYPE(NDUM).EQ.5)  GO  TO  60 
I  TYPE  -  5  MEANS  A  1/1  NODE  WHICH  MEANS 
I  F  (  0TYPE(NDUM).EQ.2)  GO  TO  70 
OTYP E  =2  MEANS  PROBABILISTIC  FIRINGS. .WE 
IF(0TYPE(NDUM).E0.4)  GO  TO  80 
IF(0TYPE(NDUM).E0.7)  GO  TO  600 
GO  TO  2C 
CONTINUE 

INSERT  READS  HERE  FOR  l/l  NODES 
WR  ITE ( 6  »  300 ) 

FORMAT! IX,  29HY0U  HAVE  INDICATED  A  1/1 
WRITE(6?301 ) 

FORMAT!  IX,  52H  INPUT  NUMBER  OF  ARCS  INPUT 
C) 


MAT022  60 
MAT022  70 
MAT022  80 
MAT02290 
MAT02300 
MAT023 10 
MAT  0  2320 
MAT02330 
MAT02340 
MAT02350 
MAT02360 
MAT02370 
MAT02380 
MAT02390 
MAT  0  2400 
MAT02410 
MAT  02420 

MUST  READ  IN  PRO  BAB  I L I TYS  MAT02430 

MAT  02440 
M AT02450 
MAT02460 
MAT02470 
MAT02480 
MAT  02490 
MAT02500 
MAT02510 

ARC  NAMEtOUTPUT  ARC  NAMES MAT025 20 

MAT02521 


WITH  A  NEGATIVE  INPUT 


WE  MUST  READ  IN  ORDERINGS 


NODE) 


59 


302 


?  0  3 
203 


306 

305 

70 


20  3 
20  1 

202 

220 


WR l T  E ! 6  »  302  )  MAT025’ 

FORMAT! IX,  23HF0RMAT  IS  I2,10!A4,A4  )  MAT0254E 

WR ITE ! 6, 30)  MAT0255G 

R  E AD ( 5 ,  303 )  MM, (YI(I ) ,Y0< I)  t 1  =  1, MM)  M  A  T  C  2  6  fc 

FORMAT ( 12, 1  ! A4, A4) )  MAT02C70 

FORMAT!  12,1  ! A4,F 6.3) )  MATC20  8 J 

DO  305  L  =  1,MARC  MAT025Or 

DO  3 06  K=1,MM  MAT02‘.0( 

I F ( ARC ( L ) .EO. Y I ( K )  )  I  ARC ( NDU M , K ) =L  MAT026I0 

I  F(  ARC!  L  )  .ElJ.YO!  K  )  )  PARC  !  NDUM  ,K  )  =L  MAT0262C) 

CONTINUE  M  A  T  0  2  ^  3  0 

CONTINUE  MAT  0264  3 

GO  TO  2  MAT0265) 

CONTINUi  MAT026FO 

INSERT  READ  FOR  <’R  OBArt  I  L  I  T  I  F  S  MAT0267 

IF1INDV.E0.  )  GO  TC  220  MAT02680 

WR  ITE  <  6 , 200  )  MAT02r,90 

FORMAT!  IX,  49HY.U  HAVE  INDICATED  A  NODE  I*  I TH  STOCHASTIC  OUTPUTS)  MAT027G 
WRITE!6,291)  MAT0271C 

FORMAT!  IX,  83H  INPUT  NUMBER  OF  OUTPUT  ARCS  NA^C  OF  UllTPljT  ARC ,  PRC  H  .M  A  T02  72r 


206 

207 

205 

80 


90 

400 

401 

402 

403 

420 


C.  ) 

WR  ITl ! 6  »  292 ) 

FORMAT!  IX,  25HF0PMAT  IS  12 , L C  !A4,F6. 3 )  ) 

WR I T  E ! 6 , 30 ) 

READ! 5, 203)  NN , ! /NAM ! I ) , ZPRC B ( I ) , I =1 , NN ) 

I F ( NN« N E. OARC I ( NUUM ) )  CALL  T  ERM ! 99 ) 

DO  205  1=1, NN 
DO  206  J  =  I , NN 
JJ  =  J 

LM=0 ARC ( NDUM, I ) 

IF!ZNAMIJ).EQ.ARCILM)  )  GO  TO  207 

CONT INUF 

CALL  TERM! 100) 

PPOA ( NDUM, I  )  =  Z  PR j  B ( J  J ) 

CONTINUE 
GO  TO  2 
NODI =NOP I +i 
NOON!  N0I3I  )  =NDUM 
GO  TO  2 

I  NOD  T= I N  3DT  +  1 
INODI! I NOOT ) =NDUM 
GO  TO  222 

IF! INDV.EQ.u)  GO  TO  420 
WR I T  E ( 6 , 40  I ) 

FORMAT! IX,  33HY0U  HAVE  INDICATED  A  1/1  BAR  NODE) 
WR I T  E ! 6 , 402 ) 


MAT 0272 1 
MAT02730 
MATC2740 
MATC275C 
MAT027  60 
MAT 02  7  70 
MAT02  7R0 
MAT  02790 
MAT  C28CT 
MATO 2  8  1  r 
MAT0282C 
MATO2830 
MATO 2 8 40 
MATC2850 
MAT 020 60 
MAT02870 
MAT  02  8  80 
MAT  02890 
MAT02900 
MAT02910 
MAT  0292  0 
MAT02930 
MAT  02940 
MAT02950 
MA  T02960 
MAT02970 


406 

405 


FORMAT!  IX,  5 1 H INPUT  NUMBER  OF  ARCS  INPUT  ARC  NAME, OUTPUT  ARC  NAME ) MAT02980 
WR I T  fc ! 6, 403  )  MAT02990 

FORMAT!  IX,  47H INPUT  A  NAME  OF  ZZZZ  FOR  THE  NO  INPUT  CONDITION)  MAT03000 
WR  IT  E ! 6 , 302 )  MAT030I0 

WR  I T  E ! 6 , 30 )  MAT03020 

READ ! 5 , 303 )  MM , ! Y I! I ) , YO !I) , I =i , MM )  MAT03030 

DO  405  L=  1 , N ARC  MAT03040 

DO  406  K= 1 , MM  MATU3050 

IF(ARC!L).EO.YI!K)  )  I  ARC  I NDUM.K ) =L  MAT03060 

IF(ARC!L).EO.YO(K) )  OARC ! NDUM, K ) =L  MAT03070 

IF(YIlK).EQoCZZ)  IARC1NDUM.K )=500  MAT03080 

CONTINUE  MAT03090 

CONTINUE  MAT  03 100 

I  ARC  I ( NDUM) =1  ARC  I  I NDUMJ  +  1  MAT031 10 


60 


GO  TO  20 

MAT03 1 20 

600 

I  F ( INDV.EQ.O)  GO  TO  666 

WRITE(6,601  ) 

MAT  03 1 30 
MATO  3 1 40 

601 

FORMAT (  IX ,  35HY0U  HAVE  INDICATED  A 
WR IT  E ( 6  *  302 ) 

WR ITE ( 6  ,602  ) 

PREFERRED  NODE) 

MAT03150 
MATO  3 1 60 
MAT03170 

602 

FORMAT! IX,  43HI/0  ARC  PAIRS  SHOULD 
WR I TE ( 6 , 403 ) 

WR IT  E ( 6 , 30 ) 

BE  IN  PREFERENCE  ORDER) 

MAT031 80 
MAT03190 
MAT03200 

666 

READ! 5, 303)  MM , ! YI! I ) , YO ! I ) , I =1 , MM ) 

no  605  L=  1 , NARC 

DO  606  K= 1 , MM 

IFIARCIL).EQ.YI(K)  )  I  ARC ! NDUM, K ) =L 
IF1ARC1D.EQ.Y01K)  )  0  ARC  !  NDUM,  K  )  =L 
IF! Y  I ! K  ). EQ.CZ Z  )  I  ARC! NDUM, K ) =500 

MATO  32  10 
MATO  32  20 
MAT03230 
MAT0324C 
MAT  032  50 
MAT03260 

606 

CONTINUE 

MAT03270 

605 

CONT INUE 

I  ARC  I 1 NDUM) =1 ARC  I! NDUM) ♦! 

GO  TO  20 

END 

MAT032  80 
MAT03290 
MAT  033  00 
MAT03J  10 

SUBROUTINE  IDIN 

MAT03320 

COMMON/IDD/  RUN  I D ( 2 0 ) 

M AT03330 

WRI TE! 6, 1) 

MAT03340 

1 

FORMAT (IX,  47HENTER  A  RUN  IDENTIFIER  OF  8C  CHARACTERS  OR  LESS) 

MAT03350 

R  EA  D ( 5 , 2 )  RUNID 

MAT03360 

2 

FORMAT! 20A4 ) 

M AT03370 

RETURN 

MAT03380 

END 

MAT03390 

c 


1 

2 

3 

4 

5 


SUBROUTINE  REPSET 
COMMON/  ITER A/ I  TER 
COMMON/ IN DEV/ I NDV 

TEST  IF  IN  BATCH  MODE f I F  SO  SKIP  GOBBELDY  GOOK 
I F ( I NDV • EQ«  0 )  GO  TO  4 
WR ITE ( 6 , 1 ) 

FORMAT! IX,  40HY0U  CAN  NOW  SET  THE  NUMBER  OF  ITERATIONS) 

WR I T  E ( 6 , 2 ) 

FORMAT! IX,  42HENTER  A  5  POSITION  INTEGER,  RIGHT  ADJUSTED) 
WR I T  E ( 6 , 3 ) 

FORMAT! IX, 1H^) 

READ! 5, 5 )  ITER 
FORMAT!  15) 

RETURN 

END 


MATO  3400 
MAT034 10 
MAT03420 
MAT03430 
MAT03440 
MAT03450 
MAT03460 
MAT03470 
MAT03480 
MAT03490 
MAT03500 
MAT035 10 
MAT03520 
MAT03530 
MAT03540 


C 


SUBROUTINE  SCAN 

THIS  ROUTINE  WILL  PRINT  OUT  THE  NET  TO  DATE 
COMMON/RANC/  IXX 

COMMON  /ARC2/  T I ME ( 500 ) , I STAT I  500 ) , PROB ! 500 ) 


MAT03550 

MAT03560 

MAT03570 

MAT03580 


COMMON/ ARC1 /ARC! 500) , INODE (5  00) , ONODE ( 500 ) , I T I  MET ( 50C ) , T ARG 1 ( 500 ) ,MAT03  590 


61 


1)0  0 


2 


50 

51 


61 

60 


CTARG215C0),  TARG3( 50^) f COSTCl  500)  ,C0STV(50C  ) 

COMMON  / PARA/  Nf\PDE»NARC 

COMMON/ N001/  NODE( 1  DO ) t I  ARC ( 1 00 1 1 0 ) ,DARC( ICO , 10 ) , PP  )  A ( 100 , 10  )  , 
C I  ARC  I ( 100),  JARCI ( 100) , ITYPEC 100) ,0TYPE(10G) ,MNIND( luC) , 

CT IMEN( 100 ) 

REAL  NODE 


INTEGER  ONODE, 

OTYPF 

WR  ITE ( 6, 1 ) 

FORMAT (  IX ,47H 

ARC 

INPUT  NODE 

JUTPUT  NODE 

time  DIST  , 

1 56H  ARG1 

DC  2  I =  1 , NARC 
LM=INODF (  I  ) 

LN  =  QNODE(  I  ) 

ARG2 

ARG3 

COST 

PROD  ) 

WRITE (6 ,3)  ARC! I ) , NODE ( LM ) , NODE ( LN ) , I T I  MET  (  I  ) , 

1TARG1  (  I  ) , TARG2  ( I  )  ,TARG3(  I  )  tCOSTC  (  I  )  ,COSTV(  I  )  , PKJ (3 (  I  ) 
3  FORMAT( 2X, A4f  6X, A4, 11X,A4, 11 X,I 1 ,5X,F7.2 ,3 (3X,F7.2 ) , 
1  2  H  +,F7.2, JX,F6.3) 


M  A  T  0  3  5  9 1 
MATO  36 1 0 
MAT  0362C 
MAT03621 
MATC  3  622 
MAT 0 36  50 
MATO 36 60 
MAT0367C 
MATO 36 80 


MAT037  K 
MAT  03  7  2  J 
M A  TO  37  3^ 


CONTINUE 

NOW  WRITF  OUT  ARCS  ANO  ASSOCIATED  DATA 

WRITE(6,50) 

FORMAT ( 1H  ) 

WR  I TE ( 6 , 5 1 ) 

FORMAT!  IX,  53H  NODE 
CS,  29H  INPUT  TYPE 
DO  60  I=1,NN0DE 

WRITE (6,61)  NODE ( I ) , I  ARC  I ( I)  ,OARC 1(1), ITYPE ( I )  ,  0  TYPE  (  I  ) 
FORMAT ( 3X,A4f 15X, 12, 15X,I2,20X, I2,15X,I2) 

CONTINUE 

RETURN 

END 


NO.  OF 
OUTPUT  TYPE) 


INPUT  ARCS 


NO.  OF  OUTPUT 


MAT0378C 
MA  TO  3  7  90 
MA  TO  38  00 
MAT  0  33 10 
MATO  38  20 
MAT  038  30 
MAT  038  40 
ARCMAT03850 
MAT  03851 
MAT  0  3870 
MA  TO  38  80 
MAT  038  90 
MAT  0  3900 
MA  TO  39  10 
MAT03920 


C 

c 


c 

c 

c 


SUBROUTINE  RUNSYS 

THIS  ROUTINE  CONTROLS  THE  RUNNING  OF 


A  NET 


COMMON/IRRZ/IREPIT 

COMMON/KKIND/KIND 

COMMON/ ARC  1/ ARC ( 500) , INODE (500) »0N0DE( 500 ) , I T I  ME T ( 5U0 ) , T ARG1 ( 500) 
CTARG2( 500 ), TARG3(500) ,COSTC( 500 ) ,C0STV(500) 

COMMON/ ITERA/I TER 

COMMON  / ARC2/  T I  ME ( 500 ) , I S TA T ( 500 ) , PROB ( 500 ) 

COMMON/RUNER/  ITERMS 
COMMON/ M I  NT /SM TIM, INSM 

COMMON/ TERNN/NODN( 30) ,NODI ,T IMEZ(IOOO) ,COSTZ (1000) ,NGDEZ (  1000  )  , 
CNCOUNT ( 30) 

COMMON/ NODI/  NODE ( 100) , IARC(  100 , 1 0 ) , 0 ARC ( 1 00 , 1 0 ) , PP J A ( 1 DO , 1 0  )  , 

Cl  ARC  I ( 100), OARCI ( 100), I  TYPE (  1  CO ) , 0 T YPE ( 100 ) » MN I ND ( 1 OC ) , 

CT  IMEN ( 100  ) 

COMMON/ PARA/NNODE,N ARC 
CCMMON/TERNI/  INODI (10), INODT 
INTEGER  OARC, OARCI ,OTYPE 

SET  NUMBER  OF  ITERATIONS  SO  FAR  TO  0 
IREP IT*0 


MATO  39  3! 
MAT  0  3940 
MAT039  50 
MATO  3960 
MATO  3970 
, MAT  0  3980 
MATO  39  8 1 
MAT04000 
MAT  040  10 
MAT04020 
MAT  04030 
MAT  04040 
MAT  04 J41 
MAT04C60 
MAT  040  61 
MAT04362 
MAT04G90 
MAT04100 
MAT  0  41 10 
MAT04120 
MAT  04 l 30 
MAT04140 
MAT04150 


62 


oon  r>  n  o  oc>  o  r>  o  o  ci 


SMTIM=999999.0 

MAT 04  L  60 

ABOVE  INITIALIZES  SMALLEST  TERMINAL  FINISH  TIME  TO  A  LARGE 

NUMBER  MAT 041 70 

MAT  04 1 80 

THE  FOLLOWING  CODE  SETS  THE  INDICATER  OF  WHETHER  OR  NOT 

MAT0419C 

A  TERMINAL  NODE  HAS  BEEN  FILLED  TO  INDICATE  NO. 

MAT04200 

MAT042  10 

I T  ERM  S=  0 

MAT04220 

NOW  FIRE  INITIAL  NODES 

MAT04230 

MAT04240 

DO  10  1=1,1 NOOT 

MAT0425C 

LM= INOD I ( I ) 

MAT0426C 

T I  MEN ( L  M ) =0 . 

MAT04270 

IF(OTYPE(LM).EQ. 1)  CALL  ALLFIR(LM) 

MAT04280 

IF(OTYPh(LM).EO»2)  CALL  PROFIR(LM) 

MAT04290 

CONTINUE 

MAT  04300 

CALL  ARCCHK 

MAT  043  10 

CALL  NODCHK 

MAT04320 

I F( K IND.EO. J )  GO  TO  60 

MAT04330 

I  F ( I  TER  MS. E 0.0)  GO  TO  l 

MAT043  40 

CALL  END  I T ( KEY ) 

MAT04350 

I  F ( KEY. EO.O )  GO  TO  1 

MAT0436G 

CALL  PTERM 

MAT  04  3  70 

I REP  I T=  IREP IT  + 1 

MAT04380 

I  F (  IREPIT.EQ.ITER)  CALL  SGRAPH 

MAT  04390 

I  F ( IREP  IT.EU. ITER)  RETURN 

MAT  04400 

DO  40  I =1 ,NARC 

MAT04410 

TIME! I ) =0. 

MAT04420 

I  STATU  )  =  0 

MAT04430 

CONTINUE 

MAT  04440 

DO  50  1=1, NNODE 

MAT04450 

T IMENI I  )  =  0. 

MAT  04460 

CONTINUE 

MAT  04470 

GO  TO  21 

MAT04480 

COME  HERE  IF  NO  NODES  HAVE  FIRED 

MAT  044  90 

I F ( I  TERMS • EQ. 0 )  CALL  TERM!  3030 ) 

MAT  04500 

DO  ABOVE  IF  NO  TERMINAL  NODES  HAVE  BEEN  FILLED 

MAT04510 

GO  TO  61 

MAT  045  20 

DO  ABOVE  IF  THERE  IS  A  TERMINAL  NODE-END  ITERATION 

MAT045  30 

END 

MAT04540 

SUBROUTINE  NODCHK 

SEE  WHAT  NODES  ARE  READY  TO  FIRE, FIRE  THOSE  THAT 
ARE  READY 

COMMON/ PARA/NNODE,N ARC 

COMMON/NODI/  NODE! IOO) , I  ARC! 100, 10 ) ,OARC ( 100 , 1 0 ) , PPOAI 100, 10 ) , 
Cl  ARC  I ( 100), OARC  I ( 100) , I  TYPE! 100 ) ,OTYPE ( 100 ) , MN I ND( 100) , 

CT IMENI 100) 

COMMON  /ARC2/  T I  ME ( 500 ) , I S TA T ( 500 ) , PROB ( 500 ) 

INTEGER  OTYPE , OARC I »  OARC 
COMMON  /KKIND/KIND 
REAL  NODE 
K  IND=0 

DO  1  1=1, NNODE 

I  L=  I 

J=0 

I  F  (  ITYPEI  D.NE.l)  GO  TO  40 


MAT04550 
MAT04560 
MAT04570 
MAT04580 
MAT  04590 
MAT04600 
MAT04601 
MAT  04602 
MAT04630 
MAT  04640 
MAT04650 
MAT04660 
MAT04670 
MAT04680 
MAT04690 
MAT04700 
MAT047 10 
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CALL  AN  DTST (  I  L  ,  J  ) 

MAT 047  20 

GO  TO  100 

MAT  0  4 / 30 

40 

I  F  (  ITYP  (  I  KNE.2)  G  0  TO  41 

M A  TO  47  40 

CALL  GRTST(  IL  ,  J) 

MAT0475C 

GO  TO  1  *C 

MAT  047  60 

41 

CONT INUt 

MAT04  /  70 

IFdTYPH  I  )  •  EO  •  5  )  GO  TO  52 

MAT0473C 

IFdTYPi  (  I  )  «  EO  •  6  )  GO  TO  53 

MAT04790 

I  F  (  I  TYPE  ( I  UFQ.7)  GO  TP  54 

M  A  T  0  4  5  0  0 

100 

IF(J.EQ.O)  GO  TO  1 

MAT049  l  J 

K I ND  =  K I  NO  ♦  1 

MAT04520 
M  A  T  0  4  M  3  > » 

NOW  ZER'  OUT  INPUT  ARCS 

MA  T 04 8 40 
MAT04R50 

IRKM  ARCI  (  I  ) 

MAT  04  9  6 

DO  22  I J=lt IRK 

MAT04H70 

LM=I ARC d  ,  I  J  ) 

MAT04BB0 

I F (  ISTAT(LM).EQ.O)  GO  TO  22 

MAT 04 8  90 

ISTAT ( L  M ) -4 

MAT0490u 

22 

CONTINUE 

MA  T  04  0  1  ; 

I  F  (  OT  YP  F  (  IJ.NE.l)  GO  TO  50 

MAT  0492 

CALL  ALLFIRd  ) 

MAT04930 

GO  TO  1 

MAT  0494 j 

50 

I  F ( 0 T YP  E (  I  )  • NE  •  2  )  GO  TO  51 

MAT  049  5 

CALL  PR  F  IR ( I  ) 

M A  TO 49 60 

GO  TO  1 

MAT  04970 

51 

I F ( OT Y P n (  I ) • NE • 4  )  GO  TO  52 

MAT  049  80 

CALL  I T  ALL (  I ) 

MAT  04990 

GO  TO  1 

MA  T  05000 

52 

IF(OTYPc( I ) • NE • 5 )  CALL  TERM11210) 

MAT05010 

CALL  ONFONE ( IL, J) 

MAT  05020 

K  I  NO  =  K IND+J 

MATO  50  30 

GO  TO  1 

MAT  0  5  040 

53 

IF(OTYPc( I ) • N E • 6 )  CALL  TERM( 1211) 

MAT05050 

CALL  ON EB AR (  IL, J ) 

MAT05J60 

K  I  ND  =  K I ND  + J 

MAT  050  70 

GO  TO  1 

M  A  T  0  5  9  80 

54 

IF(OTYPE( I) .NF.7 )  CALL  TERM(1212) 

MATO  5  090 

CALL  PREFER ( IL,J ) 

MAT  0  5100 

K I ND=K I ND+J 

MAT051  10 

1 

CONTINUE 

MATO  5 1 20 

RETURN 

MAT  05 1 30 

END 

MATO  5 1 40 

SUBROUTINE  ANDTST ( I • J )  MAT05150 

C  THIS  ROUTINE  TESTS  AND  NODES  MAT05160 

COMMON/ NODI/  NODE! 1 00 ) . I  ARC (  100 . 1C ) . OARC ( 1 00 , 10 ) , PPO A(  100 , 10  ) .  MAT05I70 
C I  ARC  I ( 100), OARC I (100) , I  TYPE! 1 00 )  ,0T YPE ( 1 00 ) , MN I N0( 100 ) ,  MAT05171 


CT I  MEN ( l  JOl 
REAL  NODE 

COMMON  /ARC2/  T I  ME ( 500 ) , I STAT ( 500 ) 
I  I  =  I  ARC  I  (  I  ) 

T IM-0« 

DO  1  K«l»  II 
KK  =  I  ARC (  I  » K  ) 

I F ( ISTATIKK ) • NE • 2  )  GO  TO  2 
IF(TIMEIKK).GT.TIM)  T I M=  Tl ME ( KK ) 


MAT05 172 
MAT05200 

PR0BI500)  MAT05210 

MAT05220 
MAT05230 
MAT0524C 
MAT05250 
MAT  05260 
MAT05270 
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1  CONTINUE  MAT052BC 

T  I  MEN (I  )  =  TIM  MAT05290 

J=I  MAT0530C 

RETURN  MAT053 1C 

2  J=0  MATQ532G 

RETURN  MAT  05330 

END  MAT05  i40 


1 


SUBROUTINE  )RTST(I,J) 
THIS  ROUTINE  TEST  OR 
COMMON/NODI /  NDDE(ICO) 
C I  ARC  I  ( 1 00 ) f  OARC  I  (  100)  » 
CT IMENI 100) 

REAL  NODE 

COMMON  / ARC2/  TIMEI500 
COMMON/ ARC  1 /ARC!  500) , I 
CTARG2I 500), TARG3I 500) , 
I  I = I  ARC  I  (  I ) 

T IM=1003000. 

J=0 


MAT053  50 

NODES  MAT05360 
, I  ARC! ICO, 10) , OARC ( 100,10) ,PPOA( 100,10) ,  MAT0537U 
I  TYPE!  100) ,OTYPE ( 100 ) , MNI ND( 100 ) ,  MATO  5 3 71 

MAT05372 
MAT  054CC 

) ,ISTAT(500) ,PR0B(500)  MAT05410 
NODE (5 00) , ONODE ( 500 ) , ITI ME T ( 500 ) , T ARG l ( 500 ) , M AT0542C 
C0STCI500) ,C0STV(500)  MAT05421 

MAT05440 
MAT054  50 
MAT05460 


DO  1  K*  1 » 1 1 

MAT05470 

KK  =  I  ARC ( I » K ) 

MAT  054  80 

IF(  ISTAT (KK ) • NE«  2 ) 

GO  TO  1 

MAT  054  90 

J=1 

MAT05500 

IF(TIME (KK) .LT.TIM] 

1  T  I  M=TI  ME  (  KK  ) 

MAT05510 

CONTINUE 

MATO  5520 

IF(J.EQ.O)  RETURN 

MAT05530 

TIMENt  I ) *  T I M 

MATO  55  40 

RETURN 

MATO  55  50 

END 

MAT  0  55  60 

SUBROUTINE  PROFIR(I) 

C  FIRE  NODE  I  USING  STOCHASTIC  CONSIDERATIONS 

COMMON  /ARC2/  TI ME ( 500 ) , I STAT ( 500 ) , PROB ( 500 ) 

COMMON/NODI/  NODE! 100) , I ARC! 100 , 10 ) , OARC ( 100 , 10 ) ,PPOA( 100,10 ) , 
Cl  ARC  I  ( ICO),  OARC I ( 100) , I  TYPE!  100) ,OTYPE( 100) ,MNIND( 100) , 

CTIMENI 100) 

REAL  NODE 
COMMON/RANC/  I XX 
INTEGER  OARC.OARCI »OTYPE 
C 

IXXX=IXX 

CALL  RANDUI IXXX, IXY.RVAL) 

I XX= I XY 
AHIGH=0. 

I  I  =OARC  I  (  I  ) 

DO  1  K* 1,11 
KK*0 ARC (  I  ,X  ) 

ALOH=AHIGH 
AHIGH=ALOW*PPOA( I,K) 

IF((RVAL.GE.AL0H).AND«(RVAL.LE«AHIGH) )  GO  TO  2 

1  CONTINUE 
CALL  TERM! 1  ) 

2  I STAT ( KK ) =1 


MAT0557C 
MAT05580 
MAT05590 
MAT05600 
MATO  560 1 
MAT05602 
MAT05630 
MAT05640 
MAT05650 
M AT05660 
MAT05670 
MAT05680 
MAT05690 
MAT05700 
MAT05710 
MAT05720 
MAT05730 
M AT05740 
MATO  5750 
MAT05760 
MAT05770 
MAT057  80 
MAT05790 
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RETURN 

END 


MAT  05800 
MATOS 8  10 


SUBRCUT  IN"  ALLF  I  R (  I  ) 

FIRE  ALL  OUTPUT  ARCS 

COMM  .  N  /ARC?/  TIMF(  500),  ISTATI500)  ,PR0B(50C») 

COMMON/  NODI  /  NCDCI  100)  ,  I  ARC  (  100, 10 ) . OARC ( ICO , 1 0 ) , PPOA ( K 0 »  10  )  , 
Cl  ARC  I ( 100 ) ,  ARC  I  ( 1 Oi  )  , I TYPFI  ICO)  ,OTYPE ( IOC ) , MN I  NO t IOC ) , 

CT I  MEN ( 100) 

REAL  NODE 

INTEGER  OARCinARC 
I  I  =  0  A  R  C  I  (  I  ) 

DO  1  K=  l ,  I  I 
K  K  =  0  A  R  C  (  I  ,  K  ) 

I STAT ( KK ) =1 
CONT INUL 
RETURN 
END 


MAT05?  20 

MAT0583C 

MAT0584C 

MAT05850 

MA  TO  5°  c  I 

MATO  58  52 

MAT058P0 

MATC589C 

MAT05900 

MAT05910 

MAT0592. 

M  A  T  0  5  9  3  0 

MAT05949 

MAT0595C 

MAT059n 


SUBROUTINE  ONEONE ( ILK, J )  MAT05970 

COMMON/ P AR A /NNODE.N ARC  MAT 05980 

COMMON/NODI/  NODE ( KO) , IARCI 100. 1 0) ,OARC( 100,10) ,PP)A( 100,10  )  ,  MATO 5990 
C I  ARC  I ( l 90), OARCI I  100) , I TYPE( 100) ,0  TYPE (100 ) ,MNIND( 1 JO) ,  MATO  5991 

CTIMEN(IOO)  MAT  05992 

COMMON/ ARC  1/ ARC! 500 ), INODE  1500) , ONODEI 500) ,  ITIMET ( 5  )i ) , TARGl < 500 ) .MAT 06020 
C  T  ARG2 I  5  00 ) , TARG3I500),C0STC( 500)  ,COSTV(500)  MAT06021 

COMMON  /ARC  2/  T I  ME ( 50  ') , I S TA T ( 500 )  , PROB ( 50C )  MAT06040 

INTEGER  OTYPE.nARC  I  ,OARC  MAT06j50 

REAL  NODE  MAT06060 

j  =0  MAT06070 

I  I  *  I  ARC  II  ILK)  MAT  06080 

T I M= 1 00  0000  *  MAT06090 

DO  1  K  =  1 ,  II  MAT0610C 

KK= I  ARC  I  I LK , K )  MAT06110 

IF! ISTAT(KK).NE.2)  GO  TO  1  MAT0612H 

IF(TIMEIKK).GE.TIM)  GO  TO  1  MAT06130 

T I M=  T I  ME  I KK )  MAT0614C 

j=K  MAT06150 

CONTINUE  MAT06160 

IF(J.EQ.O)  RETURN  MAT06170 

ZERO  OUT  ALL  INPUTS,  FIRE  J  TH  OUTPUT  ARC  MAT06180 

DO  2  K=  1, I  I  MAT 061 90 

L M= I  ARC ( I LK , K  )  MAT06200 

I F ( ISTAT(LM).EQ.O)  GO  TO  2  MAT0621C 

ISTAT(LM)  =  4  MAT  06220 

CONTINUE  MAT06230 

LM=OARC I  ILK, J )  MAT06240 

ISTAT(LM)=1  MAT  062  50 

TIMENI ILK)=TIM  MAT06260 

RETURN  MAT  062  70 

END  MAT06280 
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SUBROUTINE  ONEB AR <  I LK , J ) 

MAT06290 

COMMON/ PARA/NNODE,N ARC 

MAT06300 

COMMON/ NODI/  NODE ( 100) » I  ARC (  100, 10) , OARCI  100 « 10 )  t PPOA (100, 10 ) , 

MAT0631C 

C I  ARC  I ( 100),  ^ARCI (100) , I  TYPE ( 100)  ,OTYPE< 100 ) , MNI ND( 100) , 

MAT06311 

C  T I  MEN ( 100) 

MAT063 12 

COMMON/ ARC 1/ ARC ( 500 ) , I  NODE < 500 ) , ONODE ( 500 ) , I T I  ME T ( 500 ) , T ARG 1 (500  ) 

, MAT06340 

CTARG2( 500), TARG3( 500) ,COSTC( 500) ,C0STV(500) 

MAT06341 

COMMON  /ARC2/  T I  ME ( 500 ) , I S TA T ( 500 )  , PROB ( 500 ) 

MAT  06360 

INTEGER  OTYPE, OARCI ,OARC 

MAT063  70 

REAL  NODE 

MAT06380 

J  =0 

MAT06390 

I  1  =  1  ARC  I (  ILK) 

MAT0640C 

T I M= 1 000000* 

MAT0641C 

DO  1  K*l,  II 

MAT06420 

KK  =  I  ARC (  I LK  *  K  ) 

MAT06430 

I F(KK.EQ*500)  I  B AR  =  K 

MAT06440 

IF(KK.E0.500)  GO  TO  1 

MAT06450 

I F ( I S  TA  T ( KK ) • NE*  2 )  GO  TO  1 

MAT  064  60 

I F ( T I ME ( KK ) •  G  F  ♦  T  I  M  )  GO  TO  1 

MAT06470 

T I M=  T I M E ( KK  ) 

MATO  64  80 

J  =K 

MAT  064  90 

1 

CONT INUE 

MAT06500 

IF(J.EQ.O)  GO  TO  5 

MAT065  1C 

c 

ZERO  OUT  ALL  INPUTS,  FIRE  J  TH  OUTPUT  ARC 

MAT06520 

DO  2  K=l,  II 

MAT0653C 

LM= I  ARC (  I LK , K ) 

MAT06340 

I F ( I STAT(LM) .EG.O)  GO  TO  2 

MAT  065  50 

I  ST  AT ( LM )  =  4 

MAT  06560 

2 

CONTINUE 

MAT06570 

LM=OARC( ILK, J ) 

MAT065  80 

ISTAT(LM)=1 

MAT  06390 

T  IMEN (  I LK ) - T I M 

MAT  06600 

RETURN 

MAT06610 

5 

T I M  =0# 

M AT06620 

DO  6  K  =  l,  II 

MAT06630 

KKM  ARC  (  I  LK  ,  K  ) 

MAT  06640 

I F ( K K *  E  Q«  500 )  GO  TO  6 

MAT06650 

I  F ( I$TAT(KK).NE.3)  GO  TO  7 

MAT06660 

IF( TIME(KK)#LE«TIM)  GO  TO  6 

MAT06670 

T I M=T I  ME ( KK ) 

MAT06680 

6 

CONTINUE 

MAT06690 

LM=0 ARC ( ILK, IBAR ) 

MAT06700 

I STA  T ( L  M ) =1 

MAT06710 

J  =  1 

MAT06720 

DO  12  K*1 , I I 

MAT06730 

L  M  =  I ARC ( I LK , K ) 

MAT06740 

I  ST  AT ( LM ) =4 

MAT06750 

12 

CONTINUE 

MAT06760 

TIMEN( ILK)=TIM 

MAT06770 

7 

CONTINUE 

MAT067  80 

RETURN 

MAT06790 

END 

MAT06800 

SUBROUTINE  PREFER (  ILK » J  )  MAT06810 
COMMON/ PARA/N NO DEtNARC  MAT06820 
COMMON/ NODI/  NODE! IOC ) 1 1  ARC!  100» 10) , OARCI 100 , 10 ) , PPOA ( 100. 10 ) .  MAT0683C 
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C I  ARC  I ( K  0 ),  ARC  I ( 10(  ) , I TyPE( 100) ,OTYPt ( 100) »MNINO( IOC  >  »  MAT0603 1 

C  r  I  ME N  (  1  r0  )  MAT06R32 

COMMON/ ARC  1  /  ARC  (  500  )  ,  I  NODE  (500,  ONODE(  500  )  ,  ITIMrT  (  5  )0  )  ,  T  ARG1  (5O0),MATO6R6'' 
CTARG2(  500),  T  A0  G"1  (  50C  )  ,COSTC(  500)  ,COSTV  (500)  MATO  680 

COMMON  /ARC//  T  I  *E  (  50  )  ,  I  S  TA  T  (  500  )  ,  PROB  (  50C  )  MAT0688G 


INTEGER  GTYPE , 9ARC I ,OARC 
REAL  NODE 
J  =0 

I  1= I  ARC  I  ( ILK  ) 

T IM=0. 

DO  1  K - 1 ,  I  1 
KK  = I  ARC (  I LK, K  ) 

IF(KK. -0.5)0)  G  TO  1 
I F  (  I  ST  AT ( KK ) «  EC • 0 )  GO  TO  3J 
I F ( ISTAT(KK).EO.l)  GO  TO  30 
I F ( ISTAT(KK).F0.4)  GO  T0  3. 

IF( T IME(KK ) .GT.T IM)  TIM  =  TIME(KK  ) 
CONTINUE 
J  =  1 

IF  WE  GFT  HCRF  TlK  MODE  WILL  RE  FIRED 
FIRE  FIRST  ARC  RAIR  WITH  ?  STATUS, 

I  I  1= I  I-  l 
DO  2  K  =  1 ,  I  11 
KK= I  ARC (  I L K , K  ) 

KKK  =K 

I F (  I  ST  AT ( KK ).E 0.2)  GO  TO  4 
CONTINUE 

IF  WE  GFT  HERE  FIRE  RAR  ARC 
LM=OARC(  ILK,  I  I  ) 

GO  TO  5 

FIRE  THE  KKK  ARC 
LM=OARC (  ILK, KK*  ) 

LN= I  ARC ( ILK, KKK) 

T I M  =  T I  ME ( LN ) 

ISTAT ( LM  )=  1 
DO  40  K= 1 ,  II 
LM= l ARC (  ILK , K  ) 

I  ST  AT ( LM  )  =4 
CONTINUE 
CONTINUE 
T I  MEN ( ILK  )  =  TIM 
RETURN 
END 


MAT  06990 
MA  T 0690  5 
MAT  069  10 
MAT  06920 
MAT0693C 
MAT0694C 
MAT069*'0 
MAT 06960 
MAT0697C 
MAT  06980 
MAT06990 
MAT0700C 
MAT070  10 
MAT  07  0  20 
MAT07030 

THERE  IS  ONE  MAT07040 

M  A  T  07  J  5'j 
MAT07U60 
MAT  D707C 
MAT  07'..  80 
MAT  07090 
MAT07100 
MAT071 10 
MaT07I20 
MAT  07130 
MAT  07 140 
MAT07150 
MAT  07160 
MAT  07  1 70 
M«T07 1 80 
MAT07190 
MAT07200 
MAT072  10 
MAT07220 
MAT  07230 
MAT07240 
MAT07250 
MAT  07260 


SUBROUTINE  ARCCHK  MAT07270 

IF  INIT IATED8CHECK  PROBABILITY  OF  COMPLETION  MAT07280 

IF  C0MPLETE08CALCULATE  TIME  AND  COST  MAT07290 

M AT0730C 

COMMON/RANC/  I  XX  MAT07310 

COMMON  /ARC 2/  T I  ME ( 500 ) , I STA T ( 500 ) , PROS ( 500 )  MAT07320 

COMMON/ ARC  1/ ARC ( 500) , INODE (500) ,ONODE( 500) , I T I  ME  T ( 50C ) ,TARGl (500  )  ,MAT07330 
CTAR32( 500), TARG3 ( 500) , COSTC( 500),CCSTV(500)  MAT07331 

COMMON/PARA/  NNODE , NARC  MAT07350 

COMMON/ NODI/  NODE (100) , I  ARC (  100 , 10 ) , 9 ARC ( 1 00 , 10 ) , P PO A ( 100 , 10  )  ,  MAT07360 

Cl  ARC  I ( 100 ), EARCI ( 100) , I TYPE( 100) ,OTYPE (100) ,MNIND( 100) ,  MAT07361 

CTIMEN(IGO)  MAT  07362 

REAL  NODE  MAT07390 
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DO  1  I  =  1 »  NARC 
I F ( I  ST  AT (  I ) »NE# 1 )  GO  TO  l 
IXXX=IXX 
A  =  TARG1 ( I  ) 

B  =  TARG2 ( I  ) 

C  =  TARG3( I  ) 

I F ( I T I  MET ( I )  • EQ#  1 )  CALL  GAUS S ! I  XXX , B , A , TVAL ) 

I F ( I  TIMET ( I  ).EQ. 2>CALL  TRI AMG (I  XXX , A  ,  B ,C , TV AL ) 

I F ( I T I  MET ( I  )  .EQ.  3)  CALL  UN  I FRM ( I  XXX, A ,B ,C , TV AL  ) 

IXX=IXXX 

LM= I N0DE( I ) 

T I  ME (  I ) =TVA  L  +  T I MFN ( LM) 

I F ( PROB ( I ) . EQ. 1 )  GO  T  3  2 
IXXX=  IXX 

CALL  RAMDU!  I  XXX, I XYY.RVAL) 

I XX= I XYY 

IF(RVAL,LE.PROB! I  )  >  GO  TO  2 
ISTAT ( I  )  =3 
GO  TO  l 

2  ISTAT ( I  )  =  2 

1  CONTINUC 

RETURN 
END 


mat  07400 

MAT  07410 
MAT07420 
MAT074  30 
MAT07440 
MAT074  50 
MAT074  60 
MAT07470 

MATC748C 
MAT  07  490 
NAT07500 
MAT07510 
MAT 07b  20 
MAT  076  30 
MAT07  54C 
MAT07550 
MAT  07  560 
MAT07  5 70 
MAT07580 
MAT07190 
MAT0760C 
MAT076 10 


SUBROUTINE  I  TALL (II) 

THIS  ROUTINE  WILL  HANDLE  A  TERMINAL  NODE  BEING  FILLED 
IT  WILL  SEE  IF  THE  TIME  IS  SMALLER  THAN  ANY  OTHER  TERMINAL 
NODE  TIME  AND  IF  SO  SWAP  TIME  AND  COST  INDICATORS 

COMMON/MINT/SMTIM, INSM 
COMMON/RUNER/  ITERMS 

COMMON/ NODI/  NODE! 100) , I  ARC!  100, 10) ,OARC( LOO , 10) ,PPOA< 100,10 ) , 
C  I  ARC  I! 1 30), 3 ARC  I ! 100) , I  TYPE! 100 ) ,OTYPE ( LOO ) , MNINDI LOO  )  , 

CT IMENI 100) 

I TERMS=  1 

ABOVE  INDICATES  A  TERMINAL  NODE  HAS  BEEN  FILLED 
I F ( T IMENI I  I ).GE.SMTIM)  RETURN 
SMTIM  =  T  IMENI I  I  ) 

I NSM=  1 1 

RETURN 

END 


MAT07620 
MAT  07630 
MAT07640 
MAT07650 
MATO  7660 
MAT  07670 
MAT07680 
MAT07690 
MAT  07700 
MAT07701 
MAT  077C2 
MAT07730 
MAT  07740 
MAT07750 
MAT  07760 
MAT07770 
MAT077  80 
MAT07790 


SUBROUTINE  ENDIT(KEY) 

THIS  ROUTINE  CHECKS  TO  SEE  IF  THERE  ARE  COMPLETED  ARCS  WITH  TIMES 
SMALLER  THAN  THE  SMALLEST  TERMINAL  NODE 

IF  SO  S FT  KEY=0,  IF  NOT  SET  KEY=l  AND  TERMINATE  THIS  RUN 
COMMON  /ARC2/  T I  ME ( 500  )  , I  ST AT ( 500 ) , PROB ( 500 ) 

common/mint/  smtim.insm 

COMMON/ PARA/NNODE, NARC 

DO  l  1=1, NARC 

I F  (  ISTAT!  D.EQ.O)  GO  TO  l 

IF( ISTAT! II.E0.3)  GO  TO  l  69 


MAT07800 
MAT078L0 
MAT07820 
MAT07830 
MAT07840 
MAT  07850 
MAT07860 
MAT07870 
MAT07880 
MAT  07890 
MAT07900 
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I  F (  I  ST  AT (  I )  • EQ* 4 )  GO  TO  i  MAT07910 

I F  C  T I  ME ( I  I.LT.SMTIM)  GO  TO  2  MAT0792G 

L  CONTINUE  MAT07930 

K  EY=  L  HAT  07940 

RETUP  N  M AT0795C 

2  KEY  =  0  M  A  T  0  7  9  6  r 

RETURN  MAT  07  97  0 

END  MAT07980 


SUBROUTINE  PTCRM 
COMMON/PARA/  NN  *  DE , NARC 
COMMON / I RR Z /  I  REPIT 
COMMON/MINT/  S  MT I ^ »  I  NS  M 

COMMON/MODI/  NODE < l C 0 ) , I  ARC <  100 1 10 > , OARC < 1 00 , 10 ) , PP J A < 1  JO » 1 J ) , 

Cl  ARC  I ( 100),  ''ARC I ( 10. ) , I  TYPE!  100)  ,OTYPE (100 ) , MNINDI iOi  )  , 

CT  I M  E  N  (  1*  0) 

COMMON/ TERNN/NGONI 30 ) , NCOI , T IMEZ ( 1000) , CDS TZ ( 1000) ,NCOEZ( 1000) t 
CNCOUN T ( ^0  ) 

COMMON/ ARC  1/ ARC ( ^00 ) , INODE (500) »ONODE (500) , I TI ML  T< 500 ) , T ARG1 ( 500 
C T  ARG2 ( 5  GO ) » TARG3(500) ,COSTC( 500) ,COSTV<500) 

COMMON  /ARC/ /  T I  ME ( 500 ) , I STA T ( 500 ) , PROB ( 5 00 ) 

COMMON/ TERN  I /  I N OD I ( 1 0 ) » I N OD T 

DO  1  I  =  1 »  NODI 
JJ  =  I 

I  F ( NOON ( I ) • EQ. INSM )  GO  TO  2 

CONTINUE 

CALL  TERM ( 69 ) 

NCOUNT ( JJ )=NCOUNT( JJ ) +1 
LM=NCOUNT ( J  J ) 

I RR= I R  EP  I T  +1 
TIMEZ(  IRR) =  SM T I H 
NODEZ (  IRR )  =  I N SM 
COS TZ (  IRR) =0. 

DO  3  1=1, NARC 
LM= INODEt I ) 

I  F  (  I  ST  AT  (  I  UEQ.O)  GO 
I F ( T  IM F (  I )  ,GT.SMT IM) 


TO 

GO 


3 

i.  x  .  in.  .  i,  j  vj  l  •  jn  i  &  ri  v  v,  w  TO  4 
COSTZ (  IRR) =COSTZ( IRR) +COSTC (  I ) ♦ COS  TV ( I ) * ( T I  ME (I  )-T I  MEN ( LM ) ) 
GO  TO  3 

IF(TIMEN(LM).GT.SMTIM)  GO  TO  3 

COSTZ (  IRR) =COSTZ( I RR ) +C  OS  TC ( I )  +  COSTV(  I ) *(SMTIM-TIMEN(LM)  ) 
CONTINUE 
RETURN 
END 


MAT  07990 
M  A  T  0  8  0  0  0 
MATOS 010 
MAT08U20 
MATOS  030 
MA  T  08u  3  l 
MATOS 032 
MATOS 060 
MATOS  06 1 
r  MAT 0 80 80 
M  A  T  0  8  0  8  1 
MAT081CO 
MAT081 10 
MAT98120 
MATO 8  L  30 
MAT08L40 
MAT08150 
MATOS  1 60 
MATOS  170 
MATO  8 1 80 
MAT  08 1 90 
MAT08200 
MAT08210 
MAT082  20 
MAT08230 
MAT  08  2  40 
MAT  082  50 
MAT  08  260 
MAT  08  2  70 
MAT08280 
MAT  0  82  90 
MATOS  300 
MATO  83 10 
MATO  8  3  20 
MATOS  330 
MATO  8340 


SUBROUT INE  TE  RM ( I  ) 

THIS  ROUTINE  W  I  I  LL  ACT  AS  AN  ERROR  TERMINATOR 
COMMON/ ARC  1 /  ARC (500) , I  NODE ( 500) , ONODE ( 500 ) , I T I  ME T ( 500 ) , 
CT  ARG 1(500), TARG2(500),TARG3( 500) ,COSTC ( 500 ) , COST V( 500 ) 
COMMON/ PARA/NNODE, NARC 

COMMON  /ARC?/  T I  ME ( 500 ) , I STA T ( 500 )  , PROB ( 500 ) 

WR ITE ( 6 , 1 )  I 

1  FORMAT (  IX  ,  3 1HEX  ECUT I  ON  TERMINATED  FOR  REASON, 15) 


MAT  08  3  50 
MATOS  360 
MA  TO  83  70 
MAT00970 
MAT  00971 
M AT00990 
MAT01000 
MAT08380 
MAT08390 
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o  o  o 


MAT08400 


IF( I  • ECU  3030)  GO  TO  2 
CALL  EXIT 

2  WR I T  E ( 6  » 3 ) 

3  FORMAT (  25H  ARCS  ACTIVATED  THUS  FAR  ) 

DO  4  J  =  1 y  NARC 

IF(ISTATIJ)  • NE«  0)WRITE(6,5)  ARC(J) 

5  FORMAT (  IX, A4) 

4  CONTINUE 
CALL  EXIT 

END  MAT  0  84  20 


SUBROUTINE  R ANDU ( I  X , I Y , YFL  ) 
I F ( I  X  .NE.  65539)  GO  TO  10 
CALL  RANSET  (Y) 

10  CONTINUE 

YFL  =  RANF(Y) 

YFL  =  ABS(YFL) 

IY  =  1 

RETURN 
END 


MATO  84  30 


MAT  0  8490 
MAT  08  500 


SUBROUTINE  GAU SS ( I  X , S , AM , V ) 
A  =  O.Q 

DO  50  1=1,12 

CALL  RANDU( IX, IY,Y) 

I  X=  I  Y 
50  A=A+Y 

V= ( A-6*  0 ) *S  +  AM 

RETURN 

END 


MATG851C 
MAT08520 
MAT08530 
MAT  08  540 
MATO055O 
MAT08560 
MAT08570 
MATO850O 
MATO  8590 


SUBROUTINE  TR I ANG ( I  XT , A , B,  C,  X ) 

THIS  ROUTINE  WILL  CALCULATE  RANDUM  TRIANGULARLY  DISTRIBUTED 
VARIABLES 
I  F(  C  •  ECU  A)  GO  TO  1 
I F ( B«  ECU  A)  AM  =  0. 

IF(B.EC.A)  GO  TO  2 
AM»( B-A ) / ( C-A ) 

2  CONTINUE 

CALL  RANDU( IXT , I XY,VAL) 

I XT= I XY 

I F ( V AL#  LE# AM)  X  I =SQRT ( AM*V  AL ) 

IF(VAL.GT.AM)  XI*U-SORT(U-AM-VAL  +  AM*VAL) 

X*A+XI*CC-A) 

RETURN 
1  X  =  A 

RETURN 

END 


MAT08600 
MAT  0  8610 
M ATO062O 
MATO063O 
MAT  00  640 
MAT08650 
MATO066O 
M AT08670 
MATO060O 
MAT08690 
MAT  00700 
MAT087  10 
MAT08720 
MATO07  30 
MAT08740 
MATO075O 
M AT08760 
MAT08770 


SUBROUTINE  SGRAPH 
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MATO078O 


o  o  o  o 


DIMENSION  DOT ( I 2 r  ) 

MAT08S00 

0  I MENS  I  ON  DARK ( l  ) 

MATO 8 3 R9 

DIMENSION  ARRR  (  500r* ) 

MAT08900 

DIMENSION  ANN (10) 

MAT  0  89  10 

DIMENSION  NBLIP(30) 

MAT08940 

DIMENSION  AVAL! 10) 

MAT  089  50 

COMMON /TERNN/N  DON  (  20)  ,  NTDI  ,T  IMEZ  (1000) 

,  COSTZ ( lOoO ) , N  jDEZ ( 1000) , 

MAT  088  20 

CNCOUNT ( ^0  ) 

MAT  083  21 

PEAL  NODE 

MAT08340 

COMMON /NODI/  NODE! 100) , I ARC( 1 00, 10 ) » OARC ( 100 . 10 ) , PPU A( 100 , 10 ) , 

MAT  088  50 

C  I  ARC  I  (  i;  0) ,  ARC  I  (  IOC.)  ,  ITYPF(  100)  .OTYPE(IOO)  ,MNIND!  100) 

t 

MAT  0885 ) 

CTIMfcNl 100) 

MAT08R52 

COMMON/ ITERA/  ITER 

Mtt  T  0  8  R  9  ) 

COMMON/IDD/  R UN  i n ( 2l ) 

MAT0897l: 

DATA  AVAL/. If .2* .3.. 4,. 5,. 6 t.7,.8,.9. 

MAT08969 

DATA  DDT/4H..../ 

MAT08933 

DATA  ANN/1H1, IH2, lH3t 1H4,1H5.1H6,IH7» 

1H8 

,  1H9, 1 H  C / 

MAT08920 

DATA  DC  T/l  20*  1H  =  / 

MAT088 10 

DATA  TT/1HI/ 

MATO  8  790 

MAT08980 

TIRST  GENERATE  INDIVIDUAL  GRAPHS  BY  NODE 

MAT08990 

MAT09CCU 

MAT  090  IU 

DO  I  J  =  1 » N  D  I 

MAT0902G 

LM=NODN( J) 

MATO  9C  30 

LL=NCDUNT ( J ) 

MAT  09G40 

DO  3  I Z~  1 1  2 

MAT09050 

LLK  =0 

MAT09G60 

DO  2  K*l » I  TER 

MAT  09070 

IF(NODEZ (K I.NE.LM)  GO  TO  2 

MAT09J80 

LLK=LLK+ 1 

MAT09090 

IF(IZ.bU.l)  ARRR(LLK)=TIMEZ (K) 

MAT  0  9100 

IF(  I Z . t  0. 2 )  ARRK ( LLK ) =COSTZ ( K) 

MAT  091  lu 

2 

continue 

MAT09120 

IFILLK.NE.LL )  WRITE(6,109)  LLK, LL 

MAT  091 3C 

109 

FORMAT (IX,  25HVALUES  OF  LLK  AND  LL  ARE  , 

215) 

MAT  091 40 

IFILLK.NE.LL )  CALL  TFRM(444) 

MAT09150 

IFtLLK.EQ. j)  GO  TO  3 

MAT09I60 

CALL  GRAPH! ARRR tLL ) 

MAT  09170 

WRITE! 6» 10) 

MAT09I80 

WRITE! 10) 

MAT091 90 

IF(IZ.EO.I)  WR I T  E ! 6» 1 1 )  NODE(LM) 

MAT  09200 

IFtIZ.EO.2)  WR I TE ( 6 » 1 2 )  NODE ( LM ) 

MAT092 10 

WRITE! 6, 10 ) 

MAT09220 

WR I T E ( 6 »  25 )  RUNID 

MAT09230 

25 

FORMAT! 20X  t  2  5 A4 ) 

MAT09240 

10 

FORMAT! 1H  ) 

MAT09250 

11 

F0RMATI20X,  44HGRAPH  OF  COMPLETION  TIMES 

FOR  TERMINAL 

NODE  , A4 ) 

MAT  09260 

12 

FORMAT! 20X,  44HGRAPH  OF  COMPLETION  COSTS 

FOR  TERMINAL 

NJDE  v  A4 ) 

MAT09270 

3 

CONTINUE 

MAT  092  80 

1 

CONTINUE 

MAT09290 

CALL  GRAPH! T IMEZ, ITER ) 

MAT  09300 

WRITE(fc»lO) 

MAT093  10 

WRITE! o, 10) 

MAT  093  20 

WRITE! 6,65) 

MAT093  30 

65 

FORMAT ( 20X ,  39HGRAPH  OF  COMPLETION  TIMES 

FOR  ALL  NODES) 

MAT09340 

WRITE! 6, 10 ) 

MAT 093  50 

WRITE! 6,25)  RUNID 

MAT09360 

CALL  GRAPHICOSTZ, I  TER) 

MAT09370 

WRITE! 6,10) 

MAT09380 
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WRITE! 6, 10 ) 

MAT09390 

WRITE! 6,66) 

MAT094  00 

66 

FORMAT ( 2QX ,  39HGRAPH  OF  COMPLETION  COSTS  FOR  ALL  NODES) 

MAT094 10 

WRITE! 6 ,10) 

MAT0942U 

WR I TE ! 6,25 )  RUNID 

MAT09433 

WR I T  E ( b  » 111) 

MAT09440 

1  1  l 

FORMAT ( 1H1 ) 

MAT09450 

DO  200  JJ=1,N0DI 

MATO 94 6  > 

NCC  =  NC  UNT! JJ  ) 

ROGE  =  FLOAT! NCC )*1000.0/FL? AT! ITER) 

TOG  =  RCGF  -  FLOAT!  I F I X ( ROG E ) / 1 ) 

I F l DOG  .LT.  .5)  GO  TO  80 

ROGE  =  ROGE  +  1.0 

MAT  094  70 

BO 

CONTINUE 

200 

NBLIP(JJ)  =  ROGE 

DO  201  J  J=  1  *  NOD  I 

JZ  =  NBLIP! JJ ) /IDO 

JY  =  (NBLIP(JJ)  -  JZ*100)/10 

JX  =  (NfiLIP(JJ)  -  JZ*100  -  JY*lu) 

MAT  0949  0 

I F ( J  Z • F0«  0 )  JZ=10 

MAT09520 

IF(JY.EO.O)  J Y=10 

IF! JX  .1 0.  0)  JX  =  10 

MAT  095  30 

KK=NODN ( J J ) 

MAT09540 

LL=NBL  IP! JJ  ) 

MAT  095  50 

IF(LL.EO.O)  GO  TO  201 

LL  =  LL/10 

I F l LL  .60.  01  LL  *  1 

MAT  095  60 

WRITE! 6,  202 )  ( DOT! K ) , K  =  l, LL )  ,TT 

MAT09570 

WRITF (6,2  03  INODE ! KK ) , (DOT! K)  ,  K= 1 , LL )  ,TT,DDT, ANN! JZ) , ANN! JY ) 
1  ) 

WR I T  E ( 6, 204  ) 

f ANN( JX 

M AT09590 

WRITE! 6, 204) 

MAT  09600 

202 

FORMAT! 11X.1HI, 120A1) 

MAT09610 

203 

FORMAT (6X,A4, 1X,1HI ,120A1 ) 

MAT096  20 

204 

FORMAT!  11X, 1HI  ) 

MAT  096  30 

201 

CONTINUE 

MAT09640 

WRITE! 6, 13 ) 

MAT09650 

13 

FORMAT!  12X  ,  10  (  10H - 1  )  ) 

MAT09660 

WRITE! 6,14)  (AVAL (I), 1=1, 10) 

MAT09670 

14 

FORMAT (  13X, 10(7X,F3.  1)  ) 

MAT096  80 

WRITE! 6, 10) 

MAT09690 

WR I TE ( 6, 10  ) 

MAT09700 

WRITE! 6,15) 

MAT09710 

15 

FORMAT ( 20X ,  27HGRAPH  OF  NODE  PROBABILITIES) 

MAT09720 

WRITE (6,10) 

MAT09730 

WR I TE ( 6, 25 )  RUNID 

MAT  09740 

WRITE! 6, 111) 

MAT09750 

RETURN 

MAT09760 

END 

MAT09770 

SUBROUTINE  GR A  PH ( ARR, L I M ) 

MAT097  90 

DIMENSION  ANN! 10) 

MAT09810 

DIMENSION  ARR! LIM) 

MATO  9840 

DIMENSION  C AT ( 50 ) 

MAT09850 

DIMENSION  I  CAT (50) ,NUM(60)  , DOT (120) .AVAL! 10) 

MAT09870 

DATA  DOT/  120* 1H  =  / 

MAT09890 

DATA  AVAL/. 1,. 2,. 3,. 4,. 5, .6, .7, .8, .9, 1.0/ 

MAT09880 

DATA  TERM/1HI / 

MAT09860 
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DATA  BLANK/1H  / 

DATA  ANN/1H1,  1H2, 1H3, IH4, IH5, IH6, 1H7, 1H8, 1H9, 1H0/ 
DATA  DDT/4H..../ 

H  IG=0. 

small  =  i:ooo -OP.n 
SSUM=C. 

DO  1  l  *  1  *  L I  M 
SSUM=SSUM+ARR( I ) 

I  F  (  ARR  «  II.Gr.Bir>)  BIG~ARR(  I ) 

I  F ( ARR( I) .LT. SMALL!  SMALL  =  ARR ( I  ) 

1  CONTINUE 
RANGE=B IG-SMALL 
IF(RANGE.EQ.O)  GO  TO  900 
A INT  =RA*  GE/ 25» 

DO  4  K=  1*  50 
4  ICAT ( K ) =0 

DO  2  J  =  1 ,  L  I  M 
AH IGH  =  S  BALI 
DO  3  K= 1  *  25 
KK  =  K 

ALOtJ  =  AH  IGH 
AHIGH=ALOW+A INT 
I F ( K  «  EQ  »  25 )  AH IGH=B I G 

I F (  (  ARR ( J). 06. ALOW). AND. ( ARR ( JULE.AHIGH))  GO  Tj  77 
3  CONTINU 

77  ICAT(KK)=ICAT(KK)+1 

2  CONTINUE 
SUM=0. 

DO  90  1=1,25 

90  SUM=  SUM*  I C A  T  C  I  ) 

SCAT  =  C« 

DO  91  1=1,25 

CAT  (I)= FLOAT! I  CAT II) ) /SUM 
IFICATI  D.GT.SCAT)  SCAT  =  CAT(  I) 

IF! SCAT. EO. CAT ( I) )  ISCAT=I 

91  CONTINUE 
A INS  =  »  02 

DO  7  1=1,25 
LM=25-I +1 
NUM! LM ) =C 
AH  I GH=0 • 

DO  8  J= 1 , 50 
JJ  =  J 

AL  OW  =AH I GH 
AHIGH=ALCW+AINS 

IF( ( CAT ( I). GT. ALOW). AND. (CAT! I) .LE.AHIGH) )  GG  TO  20 
8  CONTINUF 

GO  TO  7 

20  NUM! LM) =JJ*2-1 
7  CONTINUE 

VAL  =  B  IG 
WR I TE ( 6 , 101 > 

WR I TE ( 6 , 822  )  VAL 
822  FORMAT! IX, F9.1, IX, 1HI ) 

DO  10  1=1,25 


LM 

=  25- 

1  +  1 

JZ 

=  C  AT 

(LM) 

♦1 

0. 

JY 

=  C  A  T 

(LM) 

*1 

oo.-jz*n 

JX 

=  C  AT 

(LM) 

♦1 

000.-JZ*100. 

F  ( 

J  X  •  E 

0,0) 

J 

X  =  10 

MAT09F30 
MATC9H2C 
MATC9800 
MAT0990' 
MAT099 10 
MATC992P 
MAT0993C 
MAT  G9940 
MAT0995'? 
MAT0996C 
MATC9970 
MAT  09980 
MAT09990 
MAT  1C jCO 
MAT  10  1 10 
MAT  10020 
MAT  10330 
MAT  10040 
MAT  10. 50 
MAT  100  60 
MAT  10070 
MAT  10  380 
MAT  10390 
MAT  10100 
MAT  10 1  10 
MAT  10120 
MAT  10 130 
MAT  1014  0 
MAT101 50 
MAT  10160 
MAT  10170 
MAT  10180 
MAT  10 190 
MAT  10200 
MAT  1 02  10 
MAT  10220 
MAT  102  30 
MAT 1024J 
MAT  10250 
MAT  10260 
MAT  10270 
MAT  102  80 
MAT  10290 
MAT  10300 
MAT  103 10 
MAT  10320 
MAT  10330 
MAT  10340 
MAT  103  50 
MAT  10360 
MAT  10370 
MAT  10380 
MAT  10390 
MAT  10400 
MAT  104  10 
MAT  10420 
MAT10430 
MAT  10440 
MAT  10450 
MAT  10460 
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IF(JY.EO.O)  JY  =  1C 
IF(JZ.EQ.O)  JZ  =  10 
VAL=VAL-A INT 
101  FORMAT!  IH1) 

NUM8=NUM( I ) 

I  F  (  NUMB  -»GT.  110)  NUMB=110 
IF(NUMB.EQ.O)  GO  TO  15 
NUMB  =  NUM8- 1 

WRITE! 6, 12)  CD0T(K)»K=1»  NU  M  8 ) .TERM 

12  FORMAT!  IX, 1DX, IHI , 120A1) 

W  R  I  TE(  6 ,  11  IVAL,  !O0T!K),  K=  1 »  NUM8  )  ,  TERM  ,  BLANK  ,  DDT  ,  ANN! JZ )  ,  ANN!  JY) 
C , ANN ! JX  ) 

ll  FORMAT! IX, F9.1, IX, IHI , 120A1) 

GO  TO  10 
15  WRITE(6»80) 

WRITE!6,81)  VAL 
80  FORMAT!  11X,  1HI  ) 

8  l  FORMAT!  IX,  F9.  1  ,  1  X  ,  1>!  I  ) 

10  CONTINUE 

WR I T  E ( 6 , 1 3 ) 

13  TORMAT!  I  2  X,  10!  10H - I)  ) 

WR  ITE! 6, 14)  ! AVAL!  I  ) , 1  =  1 , 10) 

14  FORMAT! 13X, 10! 7X.F3, 1 )  ) 

AMEAN  =  S  SUM/L I M 

SS0=0. 

00  800  1=1, LIM 

800  SSQ=SSQ  +  ( AMEAN-ARR! I ) )**2 
SDUM=0. 

DO  801  1=1,25 
SDUM  =  SDUM  +  C AT (  I ) 

I F ! SDUM . L  T»  0«  5 )  GO  TO  801 

AMED=SMALL+ !I-1)*AINT+AINT*( SDUM-CAT ( I ) )/SDUM 
GO  TO  8  02 

801  CONTINUt 

802  CONTINUE 

AMODE=SMALL+AINT*! I  SCAT- 1 ) +A I  NT/2. 0 

AVAR=SSO/LI M 

WRITEI6.810) 

810  FORMAT! 1H  ) 

WR  ITE 16 ,811)  AMC AN, AVAR, AMEO , AMC  DE 

811  FORMAT!  13H  THE  MEAN  IS  ,F10.2,  18H  THE  VARIANCE  IS  ,F1G.2, 

C  16H  THE  MEDIAN  IS  ,F10.2,  14H  THE  MOOE  IS  ,F10.2) 

RETJRN 

900  WRITEI6.101) 

WRITE! 6,902 )  RIG 

902  FORMAT!  IX,  49H  ALL  VALUES  IN  THE  ARRAY  ARE  IDENTICAL  .AND  ARE  =» 
C  F12.4) 

RETURN 

END 


MAT10470 
MAT10480 
MA  T 1 0490 
MAT  10500 
MAT  105 10 
MATIO520 
MAT10530 
MAT  10540 
MAT10550 
MAT10560 
MAT10570 
MAT  105  71 
MAT  10590 
MAT lOoOO 
MAT  106 10 
MAT l 0620 
MA  T 1 06  30 
MAT  10640 
MAT10650 
MAT  10660 
MAT  106  70 
MAT  1 06  80 
MAT  10690 
MAT  10700 
MAT  1 07 10 
MAT10720 
MAT  107  30 
MAT  107  40 
MAT  107  50 
MAT10760 
MAT  10770 
MAT10780 
MAT10790 
MAT  10800 
MAT108 10 
MAT10820 
MAT10830 
MAT10840 
MAT10850 
MAT  10860 
MAT  108  70 
MAT  10871 
MAT  10890 
MAT10900 
MAT10910 
MAT  10920 
MAT  1092 1 
MAT10940 
MAT  10950 


BLOCK  DATA  MAT10960 

COMMON/IDD/  RUN  ID ( 20 )  MAT10970 

COMMON/PARA/  NNODE , NARC  MAT10990 

C0MM0N/ARC2/  DUM ( 500 ) , I DUM I  500 ) , DUM A < 500 )  MAT11010 

COMMON/ ARC  1 /  BDUMI 500 ) , JDUMI 1 500 ) , CDUM ! 2 500 )  MAT11030 

COMMON/NODI/  KDUM ! 2 100 ) , Z DUM < 1 000 ) , LDUM < 500 ) , YDUM 1 1 00 )  MAT11050 

COMMON/ TERNN/LLDUM! 31 ) , A8RB! 2000 ) , KKKL ! 1030) 

COMMON/TERNI/  KJDUMI11)  MAT11090 
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o  o  o  o  r»  o 


K AT  1 1 1 00 


DATA  KJDUM/11*0/ 

DATA  NN ‘DE/j/fNARC/v / 

DATA  DUM/50  >*0*C/ 

DATA  I DUM/500*  0/ 

DATA  DUMA/5  30*0*0/ 

DATA  BD  UM /5 00* 0*  0/ 

DATA  JDUM/150^*C  / 

DATA  CD UM /2 500*0*0/ 

C ATA  KDUM/2 100*0/ 

DATA  L  DUM/n00*0.0/ 

DATA  LDUM/5  J0*0/ 

DATA  YDUM/1 00*0*0/ 

DATA  LL  DUM/ 31*0/ 

DATA  ABBB/2 000*0*0/ 

DATA  KKKL/1 330*0/ 

DATA  R UN  I D/20*4H  /  MAT1098C 

END  M  A  T 1 1 L 1 3 


SUBROUTINE  DUMP1 
COMMON/ P A RA/NNODE,N ARC 

COMMON  /ARC  2 /  TI ME ( 503 ) « I S TA T( 500) «PROB ( 500 ) 
WRITE! 6, 1 ) 

1  FO RM  AT  (  IX  ,  2 1 HC  3  MING  THROUGH  NODCHK  ) 

WR I T  E ( 6 , 2  )  (  ISTAT(K) ,K*I ,  N  AR  C ) 

2  FORMAT (  1015 ) 

RETURN 

END 


M  A  T 1 1 1  20 
MAT  1 1 1 3  C 
MAT  III  40 
MAT  11150 
MArill60 
MAT  1 1 1 70 
MAT  11180 
MAT  1 1 1  90 
MAT1 1200 


SUBROUTINE  DUMP? 

COMMON/PARA/  NNODE  f N ARC 

COMMON  /ARC  2/  TI ME ( 500 ) t I STA T ( 500 ) t PROB ( 500 ) 
WR  IT  E  (  6 , 1 ) 

1  FORMAT! IX,  21HC0MING  THROUGH  ARCCHK ) 

WRITE! 6, 2)  ! ISTAT(K) ,K  =  I,NARC) 

2  FORMAT (  1015) 

RETURN 

END 


MAT11210 
MAT  112  20 
MAT  1 1230 
MAT  1 1 2 AO 
MAT  1 1 2  50 
MAT  1 1 2  60 
M AT 1 1270 
MAT11280 
MAT  1 1290 


SUBROUTINE  UN  I FRM ! I  XT , A , B , C , X ) 

THIS  SUBROUTINE  WILL  CALCULATE  UNIFORM  DISTRIBUTED 
VARIABLES 

PA  ADDITION  TO  PROGRAM  AT  8  DEC  70 
CODE  NUMBER  IS  3 

IF!  C  .EG*  A  )  GO  TO  1 
IF!  B  *NE.  3.0  )  CALL  TERM!  1001) 

CALL  RANDU! IXT,IXY,VAL) 

I  X T  =  IXY 

X  =  A  +  VAL  *  <C  -  A  ) 

RETURN 
1  X  =  A 
RETURN 
END 
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APPENDIX  I  I 

GENERAL  FLOW  CHART  AND  LISTING  OF  RISCA 
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This  appendix  Includes  a  general  processing  flow  chart  for 
RISCA,  but  It  does  not  Include  a  description  or  flow  chart  of  the 
Individual  subroutines.  A  detailed  description  and  flow  chart  of 
the  Individual  subroutines  are  provided  In  Appendix  III.  In  addition 
to  the  flow  chart,  a  complete  program  listing  Is  provided.  This 
version  of  RISCA  Is  designed  to  run  In  a  batch  mode  only,  although 
the  Interested  user  could  certainly  adapt  this  program  for  a  time 
sharing  mode. 

RISCA  Is  written  in  FORTRAN  IV.  Even  though  FORTRAN  IV 
is  considered  to  be  a  standard  language,  adapting  the  program  for 
any  computer  wt I  I  probably  require  minor  program  modifications. 

These  modifications  generally  result  f  rom  pecu I iar i t i es  of  the  given 
system. 


Unlike  MATHNET,  the  RANDU  subroutine  In  RISCA  contains  its’ 
own  uniform  random  number  generating  routine.  Therefore,  no  modifica¬ 
tion  in  this  area  Is  required. 
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d> 


0- 


READ  IN 

READ  IN  ARCS 

ANAME,  AINODE, 

AND  ASSOCIATED 

ADNODE,  I  DIST,  L 

DATA 

PI  .P2.P3.C1 

C2,D 

K 


RETURN 


READ  IN  NODE 

V  * 

CALL 

NAMES  &  ASSOCIAT. 

/  * 

NODIN 

DATA 

. 

RETURN 


CALL 

CONTROLS 

RUNNING 

OF  NET 

RUNSYS 

_ 

^ 

-C 


RETURN 
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c 

c 


2 


C 

c 


c 

50 

1 


COMHON/RANC/  IXX 

program  risca 

THIS  IS  THE  MAIN  ROUTINE 


I XXX-0 

CALL  RAN DIM  I  XXX, IXYYfRVAL) 
I XX= I  X YY 

CALL  repip 
call  arcin 

CALL  NODIN 
CALL  SCAI 

call  runsys 
call  exit 

STOP 


END 

SUBROUTINE  RANDUI I , J,RAN) 

IF (I.EO.O) l-llllllll 

J*l*?5  *  J»J-( J/67108864)*67)08864 


J«J*25  *  J»J-( J/67106864 
J“  J*5  S  Jmj~( J/67108864 
RAN=FLOAT( JI/67108864,  $ 
SUBROUTINE  REPIP 
READS  IN  NAME  OF  NET  AND 


*67108864 
*67106864 
RETURN  $  END 

SETS  ITERATIONS 


MAT  00010 
MAI  00020 

MAT001 10 
MATOO 1 ?0 


MAT00f>nn 


COMMf.N/ITERA/ITER 
COMMON/  I  L»D/  RUN  I  D  (  20  ) 

REAP ( 5,2)  RUNID 
FORMAT ( 25A4 ) 

ITER-500 

RETURN 

END 

SUBROUTINE  ARCIN 

THIS  ROUTINE  WILL  SERVE  TO  READ  IN  ARCS  AND  ASSOCIATED  DATA 

C0MM0N/ARC1/ARC( 100) , I  NODE ( 100) , ONODEI 1 00 ) , I T I  MET < 1 00 ) , T AR81  ( 100 ) 
1 T ARG2 (  100)  ,TARG3(  100)f  COSTCUOO)  ,COSTV(  100) 

COMMON/PAPA/NNODErNARC 

COMMON  / ARC2/  TIME! 100) , ISTAT( 100) ,PROB( 100) 

COMMON/NODI/  NODE (1 00) , I ARC( 100, 10) ,OARC< 100, 10) ,PpOA< 100, 10) , 

1 1  ARC  I ( 100),0ARCI ( 100), ITYPE( 100) ,0TYPE( 100)  ,MNINd<  100) , 

2T I  MEN ( 100) 

REAL  NODE 

INTEGER  OARCI ,OARC 
INTEGER  ONODF. 

DATA  RETC/4HRETU/ 

READ (5, 1)  ANAME, AINODE, AONODE, 1DIST,P1,P2,P3,CI,C2,D 
FORMAT (3A4, U,3F10,0,3F10.0) 

IF(ANAME.Eu.RETC)  RETURN 

NARC-NARC+1 

ARC (N ARC)* ANAME 

ITIMET(NARC)»IDIST 

TARGKNARO-Pl 

TARG2(NARC)«P2 

T  ARG3 ( NARC ) *P3 


MAT  05180 


MAT  07920 
MAT  07940 
MAT  07950 


MAT07990 

MAT  080 1 0 
MAT  08020 

MAT  08040 
MAT  08050 
MAT  08080 
MAT  08060 
MAT08070 
MAT  08260 
MAT  08260 
MAT  0827  0 
MAT  08280 
MAT08290 
MAT08300 
MAT08310 
MAT08320 
MAT  08330 


SI 


o  o 


c 


9 

10 


1 1 


12 


20 

C 


30 


40 


COSTC(NARf.  )=C1 
COST V  I NARC ) *C2 
PROP  ( N ARC ) *D 

iF(Nnor.r:.co,n)  co  to  io 

DO  9  I «1, NNODE 
I  SAVE*  1 

I  SAVE  SIMPLY  KEEPS  VALUE  OF  I  TO  USE  0()T  Oh  JjU  l  oOp 

I K I  A  I  NODE .EQ, NODE  < I ) )  00  TO  ll 

CONTINUE 

nMonh>MNonr+i 
NODE  <  NNODE  )  *A  I  NODE. 

INODE"  (  N  A  P  C  )  *  N  M  0  D  E 
OaRC  (  NNODIL',  1 ) "NAPC 

oarci  <  NiinnE » « l 

G<>  TO  12 

INODE ( NARC )*I SAVE 

OARCI  (  I  SAVE  1*0  ARC  I  USAVD  +  I 

LM=OARO I ( ISAVE) 

OARCI  ISAVE,  I.MI-UARC 
CUNT  1 NUF 
DO  20  1=1, NNODE 
I  SAVE*  1 

ie  i  aonode ,eo,node  iim  go  to  30 

CONTINUE 

WE  HAVE  A  NEW  NODE 
NN0DE-NKOBEM 
MODE  1 NNODE ) *AONODE 
ONODE I NARC ) »NNODE 
I  ARC! NNODE, t I-NARO 
I  ARC  1 1 NNODE 1*1 
GO  TO  40 

ONODE ( NAPC )■ I  SAVE 
1  ARC  I  (  ISAVE)*I  ARC  I  (  ISAVF. ) +  1 
Lf1«l  ARCI  (  ISAVE) 

1ARCI ISAVF, LM)»MARC 

CONTINUE 

GO  TC  50 

END 

SUBROUTINE  NODIN 

THIS  ROUT  1 NF  WILL  READ  IN  NODE  NAMES  AND  DATA 

COMMON/NODI/  NODE! 100),IARCI100,10),OARC(100,10),PPOA(100,10), 

1  1  ARCI (  100) , OARCI I  100 ) , lTypEI 100) ,OTYPE( 100 > , MN I ND < 1 00 > , 

2T1MENI 100) 

REAL  NODE 
DIMENSION  DDUM(IO) 

COMMON/PARA/  nnodf,narc 

COMMON/TFRNN/  NODN ( 30 ) , NOD  I , T 1 MEZ ( 50  0 ) , COSTZ ( 500 > , NOdEZ ( 500 ) , 

1 NCOUNT ( 30 ) 

COMMON  /ARC2/  T I MC ( 1 00 ) , 1 STA T ( 1 00 ) , PROR I  1 00 ) 

COMMON/ ARCI /ARCI 100) , INOpFI 1 00) , ONODE ( 1 00 > , 1 T I  MET < 1 00 ) , T AR81 < 1 00 > 
1TAR02I 100) , TAR03 ( 100 ) , COSTC ( 1 00 ) ,  COST V  (  100) 

DIMENSION  ZNAMI 10) ,ZPROB( 10) , IZNAMI 1 0 ) , Y 1  I  1 0 ) , YO ( 1 0 ) 

INTEGER  OARC, OARCI 

COMMON/ TERN  I/I  NODI ( 10)  ,  1NODT 

integer  otype 

DATA  C7Z/4HZZZZ/ 

DATA  RETC/4HPETU/ 


MATCH 3 40 
MAI G83b0 
MATO  13 
MAI  083/0 
MAT  08380 
MAT083O0 
M  A  1 08  4  U  0 
MAT  084  1  0 
MAI  08420 
MAT08430 
MAI  08440 
MAT084S0 
MAI  08400 
MAT08470 
MAI  08400 
MA (08490 
MAT08500 
MAT08510 
MAI  08520 
MAI  08530 
MAI  08540 
MAT  08550 
MAT  08560 
MAT  085/0 
MAT  08580 
MAI  08590 
MAT  C8600 
MAT  086  1  0 
MAT08620 
MAI  08630 
MAT  08640 
MAT  08650 
MAT  08660 
MAT  08670 
MAT  08680 
MAT  08690 

MAT08710 
MAT  00600 
MAT  00610 
MAT  00620 
MAT00640 
MAT  00650 

MAT  00670 
MAT  00680 
MAT  00690 


MAT  00750 
MAT00770 
MAT  00790 
MAT  00800 
MAT  00760 
MAT  00780 
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MA 100810 


C 

20  r: l a r  ( b ,  a o )  adom,  i  mim,  tgdum 

40  FORMAT  (  A4,  1  1  ,  11  ) 

IF  (  ADUn.M'.RETC)  RETURN 

no  mo  i*i,nnode 

MDUK*  1 

IKNOIlt  (  1J  .EO.ADUM)  GO  TO  101 

100  CONTINUE 
NNODE-NNODt+l 
NDUH»NNOm 

101  I T  YRL (  NHItli )  *  1  DUll 
OTYF’K(NDOM)*IODUM 

I  F  ( 1TYFF (MDUM I . FQ. 4  I  GO  TO  90 

C  I  7  YPF  =  4  nrANS  INITIAL  NODE 

II  (  I T  VF-F  ( f icon )  .EQ.6)  GO  TO^  400 

C  1  TYPE-6  UFA  NS  WE  HAVE  A  t/1  NODE  WITH  A  NEGATIVE!  INPUT 

IT  (  ITYPF (ODUM) ,EQ.5I  GO  TO  60 

C  1  TYPE*  5  MEANS  A  1/1  NODE  WHICH  MEANS  wF.  MUST  READ  IN  ORDERINGS 

30  IF  (OT  YFF  (HDUH)  ti:0.2)  GO  TO  70 

C  GTYPL-2  MEANS  PROBABILISTIC  FIRINGS. .WE  MUST  REAp  IN  PKOBAB I L I  T  Y5 

IF  (OTYprU'bMMl  .E0.4)  GO  TO  80 
IF  (OTYPFdtDUHI  .EQ.7)  GO  TO  400 
GO  TO  20 
60  CONTINUE 

c  insert  reads  here  for  j / 1  nodes 

REAMS,  3031  MM,  1  VI  (  I  I  ,  YOU)  ,  1-1/Mm) 

303  FORMA  1  ( 12,  1  0 ( A 4 , A 4  > ) 

203  FORMAT ( 12, 10(A4,F6.3) ) 

DO  3 Ob  L* )  , MARC 
DO  306  k«i,mh 

] F  t  ARC ( L)tEG.YI(K)>  I  ARC < MDUM, K I «L 
IFtARCU  I.EO.YO(K))  OARC ( NPUM, K ) -L 
306  CONTINI'F 
305  CONTINUF 
GO  TO  20 
70  CONTINUF. 

C  INSERT  RTAd  FOR  PROBABILITIES 

220  RLAC(b,203)  NN, < ZU AN  II I , ZPROB ( I  I , I » 1 > NN I 

IF  (NN.KE.OARCI (MDUM) I  CALI  TERM ( 1  I 
DO  20b  1*1, NN 
DO  206  J-  1 ,  MFi 
JJ«J 

LM-OARCINDUM,  1  I 

IFIZNAMI J) .En.APC(LM) I  GO  TO  207 

206  CONTINUE. 

CALL  TERM ( 2  I 

207  PPOAINDUh, I I-ZPRODI Jj I 
205  CONTINUE 

GO  TO  20 

00  MODI-NODI+1 

MODEM  NODI  I  -  N  D  (J  M 
GO  TO  20 

90  INODT-INCIDT  +  l 

I  NOD  I ( INODT ) *NDUM 
GO  TO  30 
400  CONTINUE 

420  READ(b,303)  MM,  ( Y I ( I  I ,  YO ( I ) , I » 1 , MM  I 

DO  40b  L«1,NARC 


MA  I  0  1  l  4  (' 
MAI n  1  ISO 
MAT01  160 
M  A  1  0  1  1  /  0 
MAT  0  1  160 
M A  TCI  190 
M  A  I  0  1  2  0  0 
MA I  0  I  2  I  0 

MA I  0  1220 
M  A  I  C  1  2  3  0 
MA  )  01240 
MAI  01250 
M  A  I  0  1  2  M  i 
MAT  C12/0 
M  A  I  0  1  ?  o  0 
MA I C I  290 
MA  I  013UI1 

MA I  0  1  320 
MA I  01330 
MA I  0 1 340 
M  A  I  0  I  4  2  0 
MA I  0  1  430 
M  A  I  0  1  4  4  0 
MAT  0 1 460 
MA I  0 1 460 
M  A  I  0  1  4  /  0 
MA I  0  I 4»0 
MA  T  0 1 490 
MAI  01500 
MAT  0  1510 
MA  F  01520 
M.A1  01530 
MAT  01620 

MAT  01640 
M.AI  01660 

m  A 1  0 1 6  b  n 

MAT  01670 
MAT  0  1680 
MA  T  0 1 69  0 

MAI01710 
MAI  01720 
MAT01730 
MAI C1740 
MAT  01750 
MAT  01760 
MAT  017  7  0 
MAT  01780 


MAT01900 

MA101910 
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4  0  6 
405 


PO  4  c  6  F  r  1  ,  ’T Ti 

MA I  0 1 920 

If  (ARCd.l  .  E  G  ,  Y  I  (K)  )  I  ARC  (  MLUM,  K  )  »L 

M  A  I  0  1  9  3  0 

IFIAKCII l.EQ.YG(Kl)  OARC ( NPUM, K ) =L 

MAI  0 1940 

U  (  Y  I  (f  ).|  O.C7.7.)  !  ARC  l  NOUN, K)  =  1U0 

(  OFT  I  Ml  F 

MA  T  0 1 9hO 

CONTINUE 

PA1 019/0 

I  ARC  I  (  Miim  )  si  ARC  I  (  NPl)M  )  ♦  1 

MA | 0 1 9bO 

GO  TO  20 

PA  1 0 1 9R  0 

END 

MA I 02  1  «n 

SUBROIJT I ) F  SCAN 

MA 107230 

PRIMS  A  SUMMARY  OF  NET 

Ki  A  1  H703M 

CorU'OM/ITf  KA/ITER 

COMMON/ I  H)/  RUN  I  IK  20) 

C OMMOU/R ANC /  I  XX 

PA  1 0/250 

COMMON  /  A  f  C  2  /  TIME!  100)  ,  1ST  ATUGO  )  ,PROh<  100) 

COMMON/AFCI /AFC  (  100)  »  INOnFM  00)  #  OHODE  (  1  on  )  ,  ITIMfT  *  »  ooi  »TAR61  (  100  ) 
I  TAH62( 1  no ) , TARG3( 100 ) , COSTC< 100) ,COSTVr 100) 

COflMON/PAPA/  NNOPE,  It  ARC 

Cup.uon/i  npi/  NorE(ion),fARC(ioo,iO),OAhr<ion,io),PpOA<ino,in), 
llAKei<10C),OARCI<inO),ITYPFUOO),OTYpE<l<)0),MNlNn<lOO), 

? T  J  MFf4  (  100) 

COKMOf-  /  TfkNn/  MO|)U(  >J0  )  ,  NOP  I  ,  T I  MEZ  ( 500  )  ,  COST  Z  (  50  0  )  ,  NOntZ  <  500  )  , 
ltlCOI'f.T  (  40  ) 

Rt. a i  Nr.ru 

INTEGER  f  t  OOf , CT  YPE, OARC I 

DIMENSinil  TYpK(4), INRUt  <  7  )  ,  OUTRUL  (  7  )  ,  Rl>l  <  Si> 
tiATA  TYPT/ '  NORM  1  ,  '  TR  I  '  ,  '  l!K!  I  F  '  ,  '  CON  •/ 


PA  I 0729n 
P  A  1  07300 
MAI  0731  0 


M  A  107330 
PA  I  07340 


OAT  A 
OAT  A 
OAT  A 


I  URN!./  '  ANP 
OUTl'I  L/  *  ALl 

rim  /  •  r«n 


’  TR  I 
1 ,  ’OR  ’ 
• , • PROP 


•  I  N I  T  *  ,  • 
, ’TERM* , 


1-1 

’1-1 


,  •  1-18’ , ’PREF’/ 

•  ,  •  1  -  1  B  •  ,  ’PREF1/ 


•  YES’/ 


MAI  07350 


WRITE (6, 4)  RUN  I  I) 

I  rOKMATt  1M,25X,?5A4) 

WRITE (6, 7)  I TFR 

F  UPI’AT  (  50X*  I  rt,  2X,  •  ITERATIONS'  ) 

WRITE  (6,50 

WRITE (6,  1  ) 

format  (  mo,  ix,  •  arc  imp  note  Out  mode 
i  1  ARl>2  arcs 

2 1  COMP  '  ) 

WRITE (6, NO) 
no  2  1=1, MARC 
l M» 1 NCOF  t  1  ) 
t  N-OriOOE  (  1  ) 
t.T ■  I  T  I  MRT  (  I  ) 

WRITE! 6,3)  AFC ( 1 ) , NODE (t  M ) , NODE ( UN ) t T YpF ( I  T  ) 
1 , TAPR2 ( I ) , T  AKG3 ( I ) , COSTC ( I ) , COST  V ( I ) ,PROB< 1) 
FORRATtPX, A4,5X,A4,8X, A4,PX, A4  ' ‘  ^  c‘n  ° 1 
12H  T,7X,T  4.2) 

continue 


time 

COST 


DI  ST 


ARG  1  ' 
P  OF  ' 


MAI  07360 


MAT07390 
MAT  07  400 
MAT  07410 


,TARGl  ( I  ) 
,3(3X)F10.2) ,4X,F 10.2, 3H  ♦  ,FJ0.2, 


50 


51 


NOW  WRITF  OUT 
WRITE(6,bC) 

HR  I TC  t  6, 50 ) 
FORMAT ( I H  ) 
WRITE(6,51  ) 
FORMAT ( IX, ' 


ARCS  AND  ASS C Cl ATEP  DATA 


NOPE 


mo,  of  input  arcs  no.  of  output 


MAT  07460 
MAT  07470 
MAT  0746  0 

MAT  07500 
MAT  0751 0 
MAT07520 
ARCS' , M  A  I  07530 
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OUTPUT  PULE  ' ) 


61 

60 


c 

c 


20 


10 

1 


61 


1*  INPUT  RULE 
WRITE <6, 50) 

DU  60  I*I,NN0DE  MA>07550 

IK*  IT  YPF.  (  I  ) 

IUR-OTYPF  (  I  ) 

UK  I  TL  (6,t,l  )  NOnEU)  ,  I  ARCI  (  I  )  ,0ARCI  <  I  )  ,  iNRULI  1R>  ,0UTRUL(  IOR) 


PURE AT ( 3X, A4, 15X, I  2, 15X, I?, 19X, A4, 13X, #4) 

CONTINUE  MAT  075*0 

R1  TUX-'!  NAT  07500 

END  MA  I  076UO 

SUBROUTINE  Rl'NSYS  MAI  05400 

THIS  ROUTINE  CONTROLS  THE  PUNNING  OF  A  NET  MA1C5&00 

m a  i  055  i  n 


COl'MCIN/ JPR7/IRCPIT 
CUMI'ON/KK  )  ND/K  I  MD 

CON  ON/ APCl/ARCI 100 ) , I  NODE ( ICO) r  ONont ( I  00 ) , IT IMET I  100  )  , TARGI  ( 100  ) 
IT  A  KG 2 ( I  0  0 ) , TARG3 ( 100>,COsTC(100>,COSTV(100) 

COMNOM/ITf RA/ITER 

COMMON  / AFC2/  T I  ME C 1 0 0 ) , I  ST AT ( 1 00 ) , PROB (  1  00  ) 

COMMON/RUhER/  ITERMS 
CuUNON/UlliT/SMTIM,  INSM 

C0)'M0N/TrKNN/  NOON  (  30  )  ,  NOT  I  ,  T I  HEZ  (  500  )  ,  COSTZ  (  500  > ,  NODEZ  (  50  0  )  , 
INCOUNT (30) 

COMUON/Nuri/  NODE! 100) , 1  ARC!  160, 1 0) ,OArC< 100/ 10 ) ,PpOA (100,10), 

1  I  ARC  I  (  10(1)  ,  PARC  I  (  100  )  ,  I  TYPE  (  100)  ,OTVPE(  100)  ,MN1Nd*  100  )  , 


MA I C 5520 


mat  o  5  p  5  o 

MA | 055/0 
MA) 05!  80 


MA  I  06610 
MA ) 05620 


2T I MFN (  100) 

COMMON/PAE  a/nijodf,  MARC 

COMMON/ TI RNI /  INOPI(IO), I  M  G  D  T 
INTEGER  OAKC,OAPCI  ,OTYPE 

SIT  NUMBER  OF  ITERATIONS  SO  FAR  TO  0 

1 REPIT-0 

SMT I M-999999.0 

ABOVE  INITIALIZES  SMALLEST  TtPM I M  Al  FINISH  TIME  TO  A  LARGE  NUMBER 
THE  FCI.IOUING  CODE  SETS  THE  I  Nl)  I  CATER  OT  WHETHER  Or  NOT 
A  TERMINAL  NODE  HAS  BEEN  FILLED  TO  INDICATE  NO, 

I TLPMS-O 

MOW  FIRT  INITIAL  NODES 
DC  10  I - I , I  NOD T 
L  M  *  I  N  0  D  I  (  I  ) 

T  I  MEN ( LM ) -0 , 

IF  (OTYPF(l.H)  .EO.  1  )  CALL  AI.LF  I  R  (  LM  ) 

IF  ( OTYPF (LM) ,E0.2)  CALL  PROF  I R ( LM ) 

CONT I NUE 
CALL  ARCCKK 
CALL  NODCHK 

IF (KIND. EG. 0)  GO  TO  60 
IF ( ITERMS, EQ.O)  GO  TO  1 
CALL  END  I  T ( KEY ) 

IF ( KEY  .EG. 0  >  CO  TO  1 

call  ptf.rm 
I  REP  I T- I  REP  I T+I 

IF ( IREPIT.EO.ITER)  CALL  SPRAPH 
IF( IREPIT.fcO.ITFR)  RETURN 
PO  40  I -  l  , N ARC 


MAI  056^0 
MA  105650 
MA I  05600 
MA 1056/0 
MAI  05660 
MAT  05690 
MAT05710 
MA 105730 
MAI  05720 
MAT  05740 
MAT  05750 
MA  f  057  f>0 
MAT  05770 
MAT  05790 
MAT  057R0 
MAT  05600 
MAT  058 1 0 
MAT  05820 
MAT05840 
MAT  05850 
MAT05860 
MAT05870 
MAI  05880 
MAT  05890 
MAT  05900 
MA  T05930 
MAT  05940 
MA 105950 
MAT  05960 
MAT  05970 
MAI  05980 
MAT  05990 
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T  !  me:  ( I )  =  o . 

I  MAI  (  I  >=(i 
40  CONTINUE 

OP  50  1  =  1,  HNUJif 

T I M F r J  (  I  )  =  0. 

50  CONTINUE 
00  To  20 

C  ConF  HERE  IF  HO  NODES  HA vF  FIREl) 

60  If  C ITERI1S.EQ.0i  CALL  TERm<4> 

C  DO  AHOVF  IF  HO  TERMINAL  Ni'I-TS  HAVE  BEEN  FILLED 

GO  TO  0] 

C  BO  A 0 0 V F  IF  TfiETE  15  A  TERfIML  ‘'OliE-EnP  ITF.KATIoN 

Et  P 

50DF0UTII  r*  ALLF  I  R  (  I  ) 

C  FIKF  A(  I  OOTPIIT  ARCS 

COMMON  /  ARC'?/  TlMC(  100),lSTAT(10f)),PkOp<  100) 

COMfON/MOpi/  NCDEdOO),  lAf?C(10(!,10)rOAkC<lOO,10),PpOA(lnO,ini, 

)  I  ARFI  (  I  00  )  ,OARCI  (  )  00)  ,  I  TYPE!  1  00  )  ,0TYPE(  1  On  )  ,K.NI  Np(  ion  )  , 

2TIME.NC 1 00) 
real  no  nr 

INTEGER  f*AKCl,OARC 
I  i  *  (.  A  R  C  I  (  I  ) 

Ot  1  K  =  1  ,  I  I 

kF  =0 ARC ( I , K) 

I  ST  AT  t  KK )  =  1 
1  CONTINUE 
RETURN 
END 

SUBROUTINE  P  R  0  F  I  R  <  I  ) 

FIRE  NODE  I  I  ■  $  I  f  i  6  STUCHASTIC  CONSIDERATIONS 
common  /  ARC2/  TIMEC1U0),  I $T AT ( 1 00 > , PROB < 1 00 1 

COMMON /HOT)  J  /  MOpE ( 100  > , I APC ( I  00, 1 0  > , OApCC 100, 10) , PpOA ( 100, 10  I , 
)  I  ARC  I (  100) ,OaRC  I  ( 100 ) , ITYTf  ( ICO) , 0 T  Y P E ( 1  00 ) , MNINp ( 100), 

2TIMCNC 100) 
peal  nook 
common/rauc/  ixx 

INTEGER  OARC/OAFCI ,OTYPF 
IXXX=IXX 

CALL  RAfini'I  IXXX,  IXY,RVAL1 
I  X X* I  X  Y 
AHIGH»0. 

I  I =0 A RC I ( I  ) 

DC)  1  K  *  1  ,  I  I 
KK-OARCC I ,K 1 
ALOW*  At!  I  Gfi 

AH  I  ft H ■  A l  iH'-FppOAI  I  ,  K  ) 

IF ( (RvAL.GE.ALOW) ,AHD. (RVAL.LE.AHIGH) )  CO  TO  2 

1  CONTINOF 
CALL  TERM  (  3  I 

2  I  ST  AT (Fk 1*1 
RF  TORN 
END 

SUBROUTINE  A RCCHK 

IF  INITI ATEPBCHECK  PROBABILITY  UF  COMPLETION 
IF  COMPC  F TEPBC ALCOLATT  TIME  AND  COsT 


MA106010 
M  A  I  0  6  U  2  n 
MA I  06030 
MAT  06040 
MAI  C  60  5  0 
MAT  060/0 
MA  T  06 OP •) 
MA I 0609P 

FI  A  I  0  6  I  1  0 
MAI  06120 
MA | 06 130 
M  A  T  0  6  1  4  0 
MA  I  02BH0 
M  A  |  0  2  H  1  P 


MA | 02830 
M  A  (  0  2  P,  4  0 

MA  1  (.2ftpn 
M  A  I  0  ?  o  /  n 
MA I  02860 
MAT  02000 
M  A  I  C  2  0 1>  0 
MAT  0201 0 
MA I  02020 
MA I C2950 
MAT  02960 
MA 102550 
MAT  02560 
MAI  07750 

MA I 025HC 
MA  T  02500 

MAT  02610 
MAI  02620 
MAI  02630 
MAT  02.640 
MAT  02650 
MAT  02660 
MAT  02670 
MAT  026b0 
MAT  02600 
MAT  02700 
MA  T  027  1  0 
MAT  02720 
MAT  02730 
MAT  02740 
MA  T  02750 

MAT  02770 
MATG27H0 
MAT  02790 
MAT  02100 
MAI  02200 
FIAT  022  1  0 
MAT  02220 
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c 


2 

1 


C 

c 


50 


CDMUON/RAI'C/  I  XX  M.A  1  G223I' 

common  /af.c2/  t i hi: <  looj,  istati  100)  ,prob<  i ou ) 

COMMON/ ARC  1/ APC ( 100 ) , I NOdF <100 ) , ONODE < lOO) , IT  1  MET < 100) , TARG1 (  100  ) , 

1TAKG2I 100) , T  ARG3 ( 100), COSTCI  1 00)  , COST V ( 10  0) 

common/para/  n  n  o  d  r ,  n  a  r  c  MAio?2/n 

COMMON/UOpl/  NODE  (  100)  ,  I  ARC<  100, 10)  ,OARC<  100, 10)  ,PpOA<  100,  10)  ,  M  A 'I  02200 

1 I  ARC  I  ( 100) ,OAPCI ( 100), I T  YPE ( 100) ,OTYpE( 100),MN|Np( 1 0  0 ) ,  MAI  02290 

2  T 1 MEN!  100) 

RIAL  NODE  MAI0231  O 


nc,  1  1*1,  MARC 
IXXX*1XX 
A  *  T  A  R  G  1  (  1  ) 

B*TAPG2( 1 ) 

C  *  T  A  R  G  3  (  I  ) 

IF  (  ITIMfcTU  )  .EO.l)  CALL  G  AUSS  (  I  XXX  ,  B  ,  A  ,  T  V  AL  ) 

I  F  (  I  T  I  MF.T  (  I  )  ,F0 . 2  )  CALL  TR  I  A  NG  (  1  XXX  ,  A  ,  B  ,  C ,  Tv  AL  ) 
IMITIHFT(  D.EG.3)  CALL  UNI F < I  XXX, A, b, TVAL ) 

If  (  ITIMFT ( I  )  ,FQ.4)  TVAL-A 
I XX=I XXX 
Lfi*  I  N'ODF  (  I  ) 

TIME! I )*TVAL+TIMEN(LM) 

If  (ISTAT(I).NE.l)  GO  TO  I 
if  iprofi  i ) . ro, i . ) go  to  2 

IXXX-IXX 

CALL  RANDIU  ixxx, IXYV,RVAL) 

I XX*I XYY 

IT (PVAl,LE.PROP< I  ) )  GO  TO  2 

I  S T  A T  t I ) *3 

GO  TO  1 

I  ST AT  t I ) *2 

CONTINUE 

RETURN 

EM) 

SUBROUTINE  G AIJSS (  I  X ,  S,  AM,  V  ) 

CALCULATES  NORMAL  DISTRIBUTED  VARIABLES 


MA 102320 
MA  I  02330 
MA I  0  233  0 
MAI  o?35n 
MAT  0  2  36  0 
MAT  023/r 
MAI C23M0 

SFP  7  0 

MAI  02390 
MAI  02400 
MAI  024  1 0 
MAI  02430 
MA | 02440 
MA I  02450 
MAI  02460 
MAT  02470 
MA I  02480 
MAT  02490 
MAT  02500 
MAT  0251 0 
MAT  02520 
MAI  02530 
MAT  02540 
MAT054UP 


A  *  0 , 0 

DO  50  1-1,12 
CALL  RANDIMIX,  I  Y,Y> 

I  X  *  I  Y 
A-A  +  Y 

V*< A-6,0)*S*AM 

RETURN 

END 

SUBROUTINE  TRIANG(IXT,A,B,C,X) 

THIS  ROUTINE  MIL  CALCULATE  RANDUM  TRIANGULARLY  DISTRIBUTED 
VAKI  ABLF.S 

IF  (  (C-A).EQ.O,)  CALL  TF.RM  (  9 ) 

Am*(R-A)/(C-A) 

CALL  RANPIX  IXT,IXY,VAL) 

I  XT* I XY 

IF ( VALiLE. AM )  X I ■SORT ( AM* VAL ) 

IF ( VAL.GT.AM)  XI »1. -SORT  I  1 ,-AM-VAL*AM*vAL> 

X«A*XI* IC-A) 

RETURN 

END 


MAI  054 1 0 
MAT  05420 
MA I  05430 
MAT  05440 
MA  T  05450 
MAT  05460 
MAT  05470 
MAT  05480 
MAT  07  1 1  0 
MAT  07190 
MAT  072U0 
MAI071B0 

MA  T  07 120 
MAT  07 1 30 
MA  I  07140 
MAT  07 150 
MA  T07 160 
MAT07170 
MA  T  072  1  0 
MAT  07220 
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on  o  o  o 


subrnuTii  r  mu f  t  ix,on  ,pes,vau 

C  CALCULATES  UNIFORMLY  DISTRIBUTED  VaRIAbIES 

c 

CALL  RANDM( TXr IY,Y) 

I  -  =  I  Y 

VAL=HRT+Y* (prS-OPT) 

RETURN 

End 

SUhRliUTIliK  NflJJCHK 

C  5LE  WH/T  Ui-DKS  ARE  READY  TO  FIRE, TIRE  THOSE  THAT 

APE  READY 

CGMl'ON/PAFWUliODE,  MARC 

CUhVHN/NODI/  MOpEI 100> , J APC< 100, 10),OArCC 100, 10) ,PpOA( 100, 10) , 
1 1  akci  (  mm  ,oarci  ( loo) ,  ityrei  ino)  ,otype<  loni  ,mniijd<  ioo)  , 

2T  1  MDl  t 1 nc ) 

CIJM’ION  /APC2/  TlHEt  )00),  lsTAT(  100),PROB(100) 

INTEOEP  OTTPK,OARC!,OARC 
COMMON  /I KIMJI/KIN!) 

REAL  MO nr 


K  i  N0*U 

DO  1  1*1,1 MOW 

1 1  =  I 

J='j 

IF  (  I  T  YPF  (  I  )  ,1.0.1  ) 
IP  (  ITYFEt  I  )  .Eft, 2) 
IK ( IT  YPF  <I),K0.5) 
IK  I  T  YF  1  (  I  )  ,K0.6) 
111  ITYPfl  I  )  ,EQ. 7 ) 


CALI.  AnDTST  (  IL,  J) 
CALI  ORTr.TUL,J> 
OR  TO  30 
00  TO  3) 

GO  Til  32 


IMJ.EC.O)  Cfl 

K!NB*KIND«-1 


TO  1 


22 


30 

31 

32 


GO  TO  22 
ISTAT(IM): 
ISTATU.Mh 
ISTAT(LM)' 


MOW  ZERO  OUT  INPUT  ARCS 
I RK  =  I  ARC  I  t  I  ) 

DO  22  I  J* 1 , I RK 

Lf!  =  I  ARC- 1  I,  I  J) 

I  F ( I STAT ( LM ) ,EO. 0 ) 

I F  t I  ST  AT ( L  Ml ,E0.2) 

IF (  ISTAT ( LH ) ,EQ. 3 ) 

I  F  t I  ST  AT  1 LM ) . EG • 1 ) 

CONTINUE 

inoTYpfd )  ,i:o,i ) 

II  (OTYPF  t  1  )  .110,2) 

IF  (OTYFF.t  I  )  ,K0.4) 

IF  (OTYPF!  I  )  ,r.n,5) 

IE  (flTYPFt I ) ,E0.6 ) 

IF  (OTYPF  (  1  )  ,F.G . 7  ) 

IT  (OTYPF.t  I  )  .EG. 5) 

I  F ( OTYPF ( I ) ,RQ,6) 

If  ( OTYPF  t I > , F Q . 7 ) 

CONTINUE 

RETURN 

END 

SUBROUT  I  NT  ANDTST ( I , J ) 

COHMCN/MODl/  I'ODEI  100),IARC(100,10),OArC(10D,10), PpO A  (  100,10), 
1  I  ARC I ( 100) ,nAPCI ( 100) , I  TYPE ( l 00) ,OTYPE  < 100) , MNINpl 1 0  0 )  , 

2TIMEM  1  00) 


CALL 

CALL 

CALI 

CALI. 


ALl FIR(  1  ) 
PROF  I R  C I ) 

I T  ALL  t I ) 
ONFOiin  il,  J) 
CALL  ONFBAPl  lL,.f) 
CALI.  PREFER  ( I  L ,  J ) 
KIND-KIND+J 
K I ND»K IMD+ J 

kind-kind+j 


Sll  7  0 


SF.P  7  0 
SEP  70 
SIP  70 
SEP  70 
SF.P  70 
MAI  02970 
MA I  0290  0 
MA I C299D 
M  A  I  0  3  D  U  0 
MA I  030 1 0 
MAI  03020 
MA 1 03030 


M  A  I  0  3  0 t>  0 
MA  f  030/0 
MA I  0  30  6  0 
MAT  05700 
MAT  03090 
MA I  0310  0 
M  A  I  0  3 1  1  0 
MAT  03120 
M A  T  03 1 30 
MA 1 031  40 
MA  T  03  160 
MAT  03170 

MAT  03180 
MAT  03190 
MAT  03200 
MAT03210 
MAI  03230 
Pi  A  I  03240 
MA I  03260 


MAI  03270 
MAT  03280 
MAT  03290 
MAT  03300 


MAT  03330 
MAT  03340 
MAT  03350 
MAT  04950 
MAT  0497  0 
MAT  04960 
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M  A  I  0  5  0  0  n 


c 


1 


2 


c 

c 


c 


1 


c 

c 


c 


Rl.AL  NC<r»t 

COMMON  /  A  R  C  2  /  TIME<lon>,IsTAT<100),PROB(iOO) 

I  1 » I  ARC  I ( I > 

T  I  fl  *  0  . 

no  i  K“ i / 1 1 

KK*IAKC( I,K> 

If  (  ISTAT (KK > ,NE.2>  GO  TO  2 
If  (Tl'lf  (KK).GT.TIM)  TIM-TIME(KK) 

CUNT  I  NUK 
T  I  -FN  (  I  )  =  T  I  M. 

J=  1 

RE TURN 
J»(> 

RETURN 

cun 

SUBROUTINE  ORTST<I,J> 

THIS  ROUTINE  TEST  'OR'  MODES 


MAT  07490 
M  A  I  0  5  D  3  0 
M  A  I  0  5  0  4  ft 
MATn5050 
MAT  G506n 
MAT  05070 
MAT  0  50(i0 
MAI  05100 
MAI051 10 
MAI  05120 
MA I  05130 
MA  I  05140 
MA I  05150 
M  A  I  0  5  1  f>  n 

M A  I  04550 
MA I  0  4  5/ 0 
MA  I  066»  0 


COMMON/NOni/  MODE ( 100 ) , I ARC( 100/ 10  > ,OArC( 1 00 , 1 0 ) , PpO A ( 1 00 ,  1 0 )  , 

II  ARC  I < 1001,0 ARC  I ( 100) , 1  TYPE! 1 00) ,OTYpE{ 1 00 ) ,  MN I  Np < 1 00 )  , 

2TIMEM  100) 

REAL  NODE 

COMMON  / ATC2/  T I  ME ( 1 0 0 ) , I  ST AT ( 1 00 ) , PROB < 1 00 ) 

COMMON/ ARCl/APCI 100) , INODE! 100) , OHODEt 100),ITIMET(1U0),TARG  1(100) 
1TARG21 1  DO  > , TAPG3( 1  0 0  > , COSTC ( 1 00 > ,COSTV( 100) 

I  I»I  ARC1  (  I  ) 

T  im*i  onr.no  o. 

J»0 

DO  1  KbJ,II 
KK  *  I  ARC (  1,K> 

If  ( ISTAT (KK  > ,ME,2)  GO  TO  1 
THIS  RGUTTNF  TESTS'AND'  NODES 

J»  1 

IF (TIME<KK),LT,TIM)  TIM-TIME(KK) 

CGNTINIIF 

IF(J.EO,0)  RETURN 
TIMEN1 I >-TlM 
RETURN 
END 

SUBROUTINE  OMEONE( ILK, J) 

TESTS  ONF-ONE  NODES 


MA 1 045M0 
MAT  0461‘0 

MAT  047  1  0 


MAI  05670 
MA 1 04750 
MAI  04760 
MAT  04770 
MAT04780 
MA  IC4790 
MAI  04800 
MAT  04960 
MA  T  03220 
MA I  04810 
MAT  04820 
MAI  04830 
MAT04840 
MA  T  04850 
MAT  04930 
MAT  04940 


COMMON/PARA/NNODE, NARC  MA103370 

CUMMON/NOD 1 /  MODE (ICO), IARC(100,10),0ArC(100, lO),RpOA< 100,10),  MAT  03380 

1  I  ARC  I ( 100 > ,0ARCI ( 100) , I  TYPE ( 100),OTYPE( 100) ,MNINfi< 100) ,  MAT  03390 

2T I  MEM (100) 

COMMON/ ARR1 /ARC  (  100)  ,  I  NODE  <  100  )  ,  ONODE(  1  0  0  >  ,  IT1MF.T<  1  00  ) ,  TARG1  (  1  00  )  , 

1TARG2( 100) ,TARG3( 100>,COSTC( 100) ,COSTV( 100) 

COMMON  /ARC2/  T I  ME ( 1 00 ) , I STAT ( 1 00 ) , PROb ( 1 00 ) 


INTEGEP  OTYPE,OARCI ,OARC  MAT03440 

REAL  NODE  MAT 03450 

J«0  MAI  03460 

I  I  *  I  ARC  I (  I  LK  >  MAT  0347 0 

T! M® 1000000,  MAT  034e0 
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u  u 


Fu'  1 

*;  arc  (  n  k  ,  k > 

IM  JST*T  (KK  )  .HE  .2)  GO  TO  | 

If  (  T  1  •'!  ( I  K  )  .(?P  .T  J II )  GO  TO  1 
T  I  M-T  1  ,"f  (KK  ) 

Tr-S 

)  Cl) NT  i  Ml  'f 

itiT.Li, .C)  Pf; Tiir-N 

G  7t.KO  Ill. T  All  INPUVS,  nrp  r  T||  i-otput  arc. 

no  7  k  -  i ,  1 1 

l  f’=l  ,R<  (  u  K  /  K  ) 

!  r  (1ST  AT  !  L(! )  .E'0.0  )  r.O  Tf!  2 
! ST  AT (  U. . =  4 

2  CCNTIfiur. 

I  il-OARC  (  II  K,  J1 
IM/Tll  M  )  »  I 
t :  i  pi. ini  )  =i  i  ii 
“IT  HKfJ 
rr  i' 

siip-k jut  nit  nr  r.nARi  n.K,.n 

Tl  STS  ONE  -UfJK  GAR  SOUPS 


f* a  i  034yu 
H  A  I  0  AN i  n 
haioanjo 
MA  1  CAN?'! 
M  A  I  OANAO 
ma  i  OAS.so 
MA I f ANNO 
MA I f  ANt  r 
MA  I  CAN/’  ft 
MAI  ft AN NO 
MA I  CANO 0 


MA I  ft  Aft  1  ft 
MA  I  ft  A  f  7< 
o  A  |  ft  A  6  A  ft 
M  4  I  f  A.ift.n 
M  A  I  it  3  7  7  ft 
M  a  I  ft  3  7  A  ft 


C 


1 


c 


2 


5 


Cl'M'UN/  PAR  A/Ntinrt ,  MART 

ruMi'fii./Muitt/  t^nrt  (  ic'd)  ,  i  Apni  i  on,  t  ft  i ,uApci  i on,  io )  ,ppu/. uoft ,  in ) , 
iuwri(iitoi,nAf;ci(ino),iTYpEnoof»OTYpE(ion>»MNiNnMDO>» 

?  T i MPn< ] oo) 

c (.'much/  a i  r  i  /  arc  1 1  on  i ,  troop  i )  no  i ,  miont  1 1  on ) ,  it  diet  ( vuco .  t  arui  t  lorn 

IT  ARCm  i  on  I  ,  T  AOGHirO  )  ,cost  1 1 1  On  )  ,  CORTVi  UlO  ) 

COMMON  /ARC?./  TIME(  1001#  1CTAT  (  10lM,Pn0b(  IOI)> 

INTEGER  OTYPf  /CARCI / 0 A RC 

pf  * l  hc  nr 


M  A  I  tlAOf.O 
MA I  0 Aft/ ft 
MA I  ( l  A  ft  8  ft 


M  A  I  ft  4  ft  A  ft 
M  A  1  ft  4  ft  A  ft 


J-0 

I  1  =  I  A RC  I  (  UK  ) 
t i n=i ofton no. 
no  l  k* i  f  1 1 

KE»I  ARC!  1  l.K  ,  K  ) 

ir (KK.Eo.iom  ir.Ak-K 

If  ( KK  .  F.0, 100)  on  Tn  l 

If  ( lSTAT(IK) .Up. 2)  GO  TO  1 

If  (TIME (KK)  .GP.TIM)  go  To  1 

T  I  fl  =  T  INF  (KK  ) 

.1“  K 

Gi'NT  1  MUF 

I F ( J.Eft.O  )  GO  TO  5 

ZERO  OUT  ALL  INPUTS,  FIRF  J  Til  OOTplIT 
DO  2  K*  1  ,  !  1 
Lh“ I  ARC ( II  K,K ) 

IF(  ISTAT(LM)  .EO.ft)  GO  To  7. 

I  ST  AT ( LM ) -4 
CONTINUE 
LM«OARt( II  K,. I) 

1STAT1I M)«*l 
T  I  MEN  (  ILK) -Till 
RETURN 
T I M*0 . 
no  6  k  =  i ,  1 1 
K),  ■  I  ARE  (  II  K,K) 


M  A  1  ft  4  0  N  0 
MA I  04000 
r A  I  04  07  0 
M  A  |  o  4  0  8  0 
MA I 040ftft 


MAI  04120 
MA ! 04 1  AO 
M A  I  04  I  40 
MAI 041D0 
MAI  04160 
MAT  04170 

ARC  MA104180 

M A  T  04 1 90 
MA  I  0421)0 


MAI  04220 
MA  T042A0 
MAI  04240 
MAT  042t>0 
MA  I  0434(1 
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MAT  04370 
MAT  04380 


MAI  044110 


MAI  04420 
MAT  044J0 
M  A  |  0  A  4  4  fl 
MAI 044b0 

M  A  I  0  4  4  0  0 
MATC44/0 

MA I 044N0 
MAI  0  4 ‘>(<0 
MAI  0  4  6  3  0 
MA  I  C46« 0 
MA  I  046N0 


commun/nodi /  nope<  ioo), i arc< ioo, 10) ,oarc< ioo, io>  ,proa(  ioo,  l o > , 

1 1  ARC  I < 100) ,0ARCI (100), ITYPEl 100) ,OTYPE( 100),MNlNn<100), 

2T  I  MEM (100) 

COMMON/ ARC1 /ARC ( 100) , I NOpE ( 1 00 ) , ONODE< lOO),ITI ME T<100),TARG 1(100), 
1TARG21 100) , T ARG3 ( 100 ) , COsTC( 100)  ,  COSTV< 1 0 0 ) 

COMMON  / ARC2/  TIME( 100) , ISTAT( 100) ,PROB< 100) 

INTEGER  OType,0ARCI,0ARC 

real  nope 

J-0 

I  I  “  I  ARC  I (  ILK) 

00  1  K-  I  ,  1 1 
KK  — I  ARC ( I LK, K  ) 

I F  ( KK . F.O,  1 00  )  GO  TO  1 

IF  (  ISTAT(KK)  ,EO.(l)  GO  TO  30 

IF(ISTAT<KK),EQ.l)  GO  TO  30 

1F( ISTAT(KK) .EG. 4)  GO  To  30 

CONTINUE 

J-l 

IF  WE  OFT  HERE  THE  NODE  WILL  BE  FIRED 
FIRE  FIRST  ARC  PAIR  WITH  2  STATUS,  IF  THERE  IS  ONE 
TIM-0, 

I  I  1  =  11-1 
DO  2  K-l ,  I II 
KK- 1  ARC ( ILK/K) 

KKK-K 

IF  (  T  I  HF.  ( KK  )  ,GT . T I  M  )  TIM-TIME(KK) 

I F (  ISTAT(KK) ,E0,2)  GO  To  4 
CONTINUE 

IF  WE  GET  HERE  FIRE!  BAR  ARC 
L  M-OARC ( ILK, I  I ) 

CO  TO  5 

FIRE  THE  KKK  ARC 
LM-0 ARC ( ILK, KKK) 

I  ST  A  T ( LM ) - 1 
00  40  K-l, I  I 
LM-IAFCI ILK,K) 

I STAT ( LM ) -4 

continuf: 


II  (KK.LQ.10tj)  RO  TO  6 
IF ( ISTATIKK ) , NE . 3 )  GO  TO  7 
l  F  (  T I  ME ( KK ) .LE.T I M 1  GO  TO  f> 
T  I  M-T  I  MF.  ( KK  ) 

CONTINUE 

LM-OARC ( ILK, IDAF ) 

I  STAT  ( l.M  1 « 1 
J=  1 

no  12  K-l,  1  I 
LM= I  ARC ( II K,K ) 

ISTAT (LM)-4 
CONTINUE 
T I  MEN ( ILK  I -TIM 
CUNT  I NUE 
RETURN 
F  I1!  I1 

SURR0U1  lilt  PREFER!  ILK,  J) 
TESTS  PREFERENCE  MOPES 

Common/ papa/ n nope, narc 
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oo  oooo  o  ooo 


30 


c 


1 

2 


TIMET. I  ILK  >  =T  I  h 
CONT  1 1  Ilf 
KETURK 
L  N  D 

SUKROUllt.L  I  T  ALL  (II) 

THIS  ROUTINE:  WILL.  HANDLE  A  TERMINAL  NOdE  BEING  Fll.ltD 
IT  WILL  Sft  IF  THL'  T  IMF  IS  SMALLER  THAN  ANY  OTHER  TERM  I  NAl 
NODE  TIME  AND  IF  SO  SWAP  TIME  AND  COST  INDICATORS 

COMMON/MINT/SMTIM,  INSlI 
COMMON/RUNfc  P/  ITERM5 

COMfON/NOUl/  mode ( 100) , I ARC( 100, 10>,UArC< 100, 10  > ,PpOA( 100, 10  )  , 

1  I  ARC  I  (  1  00)  ,OARCI  (  (1101,  I  Type  <  too  >,  OTYPE<  100),MNIN)i<  IDO)  , 

?T  I  MEN  (  100) 

I  I ERHS*1 

AlJOVf  INDICATES  a  TERMINAL  node  has  bErN  filled 
IK  (THEM  1  n.OE.SMTIril  rfturn 
SMI  I M*T I  HEN < l  I  ) 

INSM«I I 

RtltJRN 

END 

SUBROUTINE  ENDIT(KEY) 

this  routine  checks  to  see:  if  there  are  completed  arcs  with  times 
smaller  than  the  smallest  terminal  node 

IE  SO  SET  KET-0,  IF  NOT  SFT  KEY»1  AND  TERMINATE  THIS  RUN 

COMMON  /  ARC2/  T I  ME  ( 1 00  )  ,  I  ST  AT  (  l  uO  >  ,  PROB  <100 1 

COMHON/MIEIT/  SMT I  M,  I  NSM 

COMMON/PAR  A /NNfipF, NARC 

DO  1  HI, MARC 

IM  ISTAT(  I  ),e;o,0>  GO  TO  1 

I F (  iSTATt I >.EQ,3>  GO  TO  1 

IF  (  I  ST  AT  (  I  )  ,E(5.4)  GO  TG  1 

IF(TINE(  n.LT.SHTpt)  GO  Tfl  2 

CONTINUE 

KE  Y=1 

RETURN 

KCY  =  0 

RETURN 

END 

SUBROUTINE'  PTFRM 
DETERMINES  COST  Of  ITERATION 


MA  1  07740 
MAI  0  7  7  6  0 
MA  I  077/0 
M  A  I  0  7  7  8  0 
M A  I  C77N0 
MAI078O0 
M  A 1  0  7  8  1  0 
MA 1 07820 
MA 1 07830 


MA  I  0  781,0 
MA I  0 7  8 f •  0 

M A  I  078/0 
MAI  0/880 
MAI  0/8^0 
MA I  07 ROO 
MA I  0701 0 
MA  1063/0 
MAI  06380 
M  «  T  0 1  >  3  c<  0 
MAI  06480 
MA  T  064 1 0 

MAT  06430 
MA I C6440 
MA I  0 6 4 b 0 


MA 1064/0 
MA 106480 
MAI  06400 
MAT06500 
MAI  06510 
MA  T 06520 
MAT  06530 
MAT  06150 


commom/para/  nnode,narc 
COMMON/ ir<R7./I  REP  IT 

COMMON/MINT/  SHTIMflNSM  MAI  06160 

COMMUN/UOJ'l/  node ( 100) , I APC( J  00, I  0 1,0 ARC ( 100, 10) ,PpOA< 100, 10) ,  MAT06170 

1  I  ARC  I ( 100 ) ,OARCI (100),ITyPE(100),OTYPE(10  0),MNINp<100),  MAT06180 

2T I  MEN ( IOC) 

COMMON/TFRNN/  NODN  (  30  )  ,  NOD  I  ,  T  I  MEZ  (  50  0  )  ,  C.OSTZ  (  50  0  )  ,  NODEZ  (  50  0  )  , 

INCOUNT ( 30  ) 

COMMON/ a RC1/*PC<  100), INODE (100), ONODE ( 1 00 ) , 1 T IMETI 1 00) , TARG 1(100), 

1TARG2I 100 ) , TARG3 ( 100),COSTC(100),C0STV(100) 

COMMON  / ARC2/  TIME! 100), 1STAT( 100) ,PH0b( 100) 

COMMON  / I DD/  RUN  I D ( 20 ) 

COMMON/TFTN1/  IN0D1 ( 10) , INODT  MAT  C6220 

MAT  06230 
MAT  06240 


C 


DO  1  HI,  MODI 
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o  o 


J  J=  1  MA  I  06250 

IMMilDM  I)  .EO.INSH)  GO  TO  2  MATG62M) 

1  CONTINUE  MAI  0627  0 

Call  TFRMI69)  MAT062F.0 

2  NCOUNTF JJ) -NCOUNT<  JJM-1  MAT0629O 

LM-NCOUNT ( J J >  MAI  06300 


IKR-IREPIT+1 
T I MTZ ( IRHI-SMTIM 
6  NOntZI IPR)«INSM 
COSTZ ( IRR)»0. 

DO  3  I » 1 ,  N ARC 
LM«  I  NPDL  (  I  ) 

I  F  (  I  ST  AT  (  I ),EQ,0)  GO  TO  3 
IF  <TIMF.<  I  )  .RT.SMTIM)  GO  TO  A 

COSTZ ( I RR ) “COSTZ ( IRR)*COSTC( I ) +COSTV( I ) * ( T I  ME ( I ) -T I  MEN ( LM ) ) 

00  TO  3 

4  IF  (TIItFMI  LM)  »GT . SMT I M )  GO  TO  3 

COSTZ  (  I  RR  )  «L'OSTZ  (  IRRH-CORTCI  1  )  ♦COST  V  C  )  )  *  (  SM T  I  M-T  I  MEN  (  LM  )  ) 

3  CONTINUE 


RtTURN  MAI  06360 

END  MA| 06360 

SUBROUT  INK  SORAPH  MAT  06640 

DIMENSION  DOT (120) 

COMMON/TLRNN/  NODN  (  30  )  ,  NOD  I  ,  TI  MEZ(  500  >  t  COSTZ  (  50n )  ,  NODtZ  ( 500 )  , 

INCOUNT ( 30  ) 

REAL  NUDF.  MAT  0657  0 

CUMMON/NOJU/  NODE ( 1  CO ) , I  ARC ( 100, 10 ) , OArC( 1 00, 10 ) ,PpOA ( 100, 10 ) ,  MAT06560 

J  I  ARC  I  (  100) ,OARCI  (  100  )  ,  ITrpF.I  »  00)  ,  OtYPE(  100)  ,  MNIN|1(  100)  ,  MAI  06590 

2T  I  MEN ( 100) 

COMMON/ I  TER  A/  ITER  MAT06620 


DIMENSION  ARRR ( 500 ) 

DIMENSION  ANN (10) 

DIMENSION  NBLIP(IO) 

DIMENSION  AVAL(IO) 

COMMOM/IDD/  RUN  I D ( 20 ) 

DATA  TT/1HI/ 

DATA  DOT/120* 1H»/ 

DATA  ANN/4H1111,4H2222,  4H3333, 4H4444, 4H5555, 4H6666, 4H7777, 4H8888, 
14H9999,4H0000/ 

DATA  DDT/4H. .  .  ./ 

DATA  AVAL/. If .2, ,3, .4, ,5, .6, .7, .8, ,9, J .0/ 

MAI  06640 

FIRST  GENERATE  INDIVIDUAL  GRAPHS  BY  NOdE  MAT06650 

DO  1  J-lfNOni 
LM«NODN( J) 

LL-NCOUNTC J» 

DO  3  17*1,2 
LLK*=0 

DO  2  K-l, ITER 

I F ( NOPFZ ( K ) .NE.LM)  GO  TO  2 
LLK-LLK+1 

IF(IZ.EQ.l)  ARRR(LLK)«TIME2(K) 

IF ( 1Z.EQ.2)  ARPR(LLK)«COSTZ(K) 

2  CONTINUE 

I F ( LLK.NE.LL )  WRITE(6f109>  LLKfLL 
109  FORMAT! lXf 'VALUES  OF  LLK  AND  LL  ARE  *,215) 

IF (LLK.NE.LL)  CALL  TERM ( 6 ) 

CALL  GRAPH ( ARRR, LL) 
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wRITE<6,60) 

IFUZ.FO.l)  W  R I  T  E  (  6 , 6 1  )  MODE  <  LM  ) 
lF(IZ.En.2)  WRITE<6#62)  MODF(I.M) 

WRITE (6/ 60) 

WRITER,  26)  RUNID 

60  FORMAT ( 1H  ) 

61  FURf*AT(20X, 'GRAPH  OT  COMPLETION  TIMES  FOR  TERMINAL  NODE  ',A4> 

62  FORMAT! 2 OX,  'GRAPH  OF  COMPLETION  COSTS  FOR  TERMINAL  NODE  ',A4) 

3  CONTINUE 

1  CONT  I  NI'F 

CALI.  CRAFH(  TIMEZ,  ITEP) 

WKITE<6,10> 

10  FORMAT  I 1H  ) 

WRITEI6, I  0) 

WRITE! 6,65) 

65  FORMAT! 20X, 'fiRAPH  OF  COMPLETION  TIMES  FOR  ALL  NODES') 

WR  I  TE  <  6, 10) 

WPJTEI6, 25)  RUNID 
CALL  GRaPHUOSTZ,  ITER) 

WRlTE(6f 10) 

WRITE!6,  1.0) 

WRITE<6,66) 

66  FORMAT120X, 'GRAPH  OF  COMPLETION  COSTS  FOR  ALL  NODES') 

WRITE<6, 10) 

WRITE! 6, 25)  RUNID 
25  FORMAT (20X,25A4) 

WKITE(6, 1 1 1 ) 

111  FORMAT  <  1111 ) 

DO  20F  JJ»  1 , NOD  I 
NCC**N  COUNT ( J J ) 

200  NBLIP!JJ)bFIOAT(NCC)*100./FLOAT(1tER) 

DO  ?01  J>l,NOI)I 
JZ-NFLIPl JJI/10 
JY-NBLIP! JJ)-JZ*10 
I  F  ( JZ.FO.O)  J7-1D 
IF(JY.EO.O)  JY-10 
KK«N0DN< JJ) 

LL*NBLIP!JJ) 

IF (1.1. EG. 0)  GO  TO  201 
WRITE(6,202>  ( DOT  I K ) , K* 1 , LL ) , TT 

WRITE <6, 203)  NODE < KK >,< DOT <K  > ,K» 1 , LL >, TT , DDT , ANN ( JZ ), ANN < JY ) 
WRITE<6,204) 

WRITE(6,2.04) 

202  FORMAT! 1 IX, 1HI, 120A1) 

203  FORMAT  <6X»  A4, IX,  1HI,  1 20 A 1 ) 

204  FORMAT  (  UX,  1HI ) 

201  CONTINUE 

WRITE!6, 13) 

13  FORMAT  !  12X,  10!  10H - I)) 

HR  I TE  <  6, 14)  ( AVAL! I ) , I»l,  10) 

14  FORMAT! 13X, 10!7X,F3.1) ) 

WRITE<6, 10) 

HR1TE(6, 10) 

WRITE (6, 15) 

15  FORMAT(20X, 'GRAPH  OF  NOcF  PROBABILITIES') 

WRITE !6, 10) 

WRITE! 6,25)  RUNID 
WRITE! 6,  i  1 1  ) 


MAT  06B;r) 

MAT  C67AM1 
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MAT  07090 
MAT  07 1 00 


RF  TURN 


EM  0 

SUBRO'JTIliE  GRAPH  (  ARR ,  I  I M  » 

OlMEMSIOh  ANNtlO)  MAT09070 

DIMENSION  ARR(I.IM)  MAT09100 

COMMON /NODI/  IJODE( 100) , 1 ARC( 100, 10) ,0ARC( 100, 10) , PpOA( 100,10), 
1IAKCI(100),0aRCI(100),ITYPE(100),0TYPE(100),MNINdUU0), 

2T )  mt.n  (loot 


COMMON/ IDl)/  RUN  1  D  (  20) 

D 1  MENS  1  OH  C AT ( 50 )  MAT09110 

D 1  MENS  I  OH  I  CAT  (50)  ,NUK(6o  ) ,  DOT<  120  )  ,  A  V  At.  1 10)  MA  T  091.50 

LATA  DDT/4H. . . . / 

DATA  ANH/4H1  1 1 1, 4M2222,  4H3333,  4H4444 , 4H5555,  4H6ftft6,  4H7777, 4H8888, 


t  4H9999, 4(10000/ 

DMA  BLANK/ IH  / 

DATA  TERM/1H1/ 

OATA  AVAL/. 1, .2, .3, .4, .5, .ft, .7, .8, .9, 1.0/ 

DATA  DCT/120MH-/ 

R  1  G  =  0 . 

SMALL* 10000000. 

a  v  e. »  o .  o 

S  D  *  0 , 0 

IF  ILIM.Ftj.O)  GO  TO  900 
DO  1  1*1,1  I M 

A VL  =  A VF+ ARR ( 1  ) 

IK  ARR<  I  )  .GT.RIG)  B I  G-  ARR  (  I  ) 

IF ( ARR ( I ) .LT, SMALL)  SHALL  * ARR ( I ) 

1  CONTINUE 
PANOE-MG-SMALL 

II (RANOF.FU. 0. )  GO  TO  900 
AV£  -  AVF/FLOAT (LIHI 
A1NT-RAHGE/25, 

ItO  4  K- 1 , 50 
4  I  C  A  T  <  K  )  *  0 

DO  2  J- 1 , L I  M 

SD-SD*< (AVE-ARR(J)  )**2) 

AUlGH-SMALL 
DU  3  K  =  1 ,26 
ALOW-AMJGH 
AH  I  GH- Al.OW  +  A  I  NT 

II  (  (ARRU)  .GF. ALOW)  , AND. <ARR(J)  .LE.AHIgH)  )  I  CAT  ( K  )  ■  I  C  AT  (  ft  )  +  1 
3  CUNT  I NUF 

2  CONTINUE 

IF ( LI M.GT. 1 1  VAP-5D/( FLOAT ( LI M)-l  .0) 

I  F  (  LI  M.l.E,  1  )  VAP-0 
SD-SQRT ( VAR ) 

SUM-C  , 

DO  90  1-1,26 

90  SUM-SUM* (TLOAT ( I  CAT  ( I  )  ) ) 

INDEX-1 
DO  92  1-1,26 

92  CAT( I ) -FLOAT  < I  CAT ( I ) )/SUM 
GO  TO  03 

94  INDEX-2 

DO  91  1-2,26 

CATL  I  ) -CAT ( I J*CATtI-l ) 

91  CONTINUE 

93  A I NS=  ,  02 


MAI  09090 
MAT  09 l 20 
MA I  09140 
MA  T  09 1 50 
MA  T  091 00 


MAI  09180 

MAT09190 
MAT  09200 
MAT09210 
MA  T  09220 


MAT  09240 
MAT  09250 
MA  T  09200 
MAT  092/0 

MAT  09200 
MAT  09290 
MAT  09300 
MAI  09310 
MAT  09320 
MA  T  09330 
MAT  09340 


MAT  09350 
MAT  09360 


95 


r.r  7  1*1,  26 

m  a  i  n  9  4  i  n 

|.  M  *  2  -  1  ♦  l 

M  A  l  0  9  4  ?  n 

NHM  ( I.M  )  *0 

MAI  00430 

A  H  l  G  H  *  C  « 
r.tj  n 

M  A  1  0  0  4  4  0 

JJr  J 

M A  1  0  94  Mi 

A i  0W“AK  1  Gil 

MAI  094/0 

AH  I  GH*  A LOU  +  A  I  liS 

M  A  1  0  0  4  r.  0 

If  (  (CAT  (  I  )  .  GT  , ALOW )  .  Af'D  .  (CAT  (  1  )  .LC.AHIgM)  )  GO  TO  2u 

M  A  1  0  o  4  o  n 

a 

COMT  1  NUF 

MAI  0  9*S  u  n 

GO  TO  7 

MA  T  095  1  r, 

?o 

Ut*l<LM«J,r*2-l 

m a  i  0C|52  0 

7 

( .U*JT  I  MUF 

M  A  1  (;  9  b  3  0 

V  A|.*h  I  ft  +  A  J  NT 

MAI ft 

wkiti inn 

M  A  I  0951*0 

ro  in  i=i,2f 

l  •■1  =  26-1  +  1 

J Z*C  AT  ( I.M  )  *  1  of 

JV=CAT  f|Ji) « 100. -FLOAT (J7) • Id. 

jx=cat(  1 1!)  *  in oo. -flu at  ( j7» •  i oo. -Float (  jy>*  to. 

MA ) 09560 

IF ( JY.EGl.P)  JY-10 

iKJZ.r.o.oi  .iz»in 
in jy.ni.o i  jy=io 

MA  1  ('9500 

V  A  L  =  V  A  t  -A  J  NT 

MA  1  096  in 

lot 

FORMAT  (  1 M  1  ) 

M A  T  0 96 ? 0 

NUMP=NUM( 1  ) 

MAT  09630 

II  (NUME-.GT.  110)  NUMB*  1 1  0 

MAI  0964ft 

IF  (miHF«F:O.OI  GO  TO  15 

M A  1  0  96b  ft 

NUMf>=  NUMB  —  1 

M A  1 09660 

WK1TF  (6,  12)  <nOT(K>  ,  K»l, NUMB)  ,  TERM 

MAI  09670 

12 

FORMAT < IX, 12X, lHI, 120A1 ) 

WRlTElftf 11 ) VALf (D0T(K1 f K«J ,MUhb) , TERM, BLANK, DDT , ANN ( JZ >, ANN ( JY  ) 

1  ,  A  N  M  (  J  X  ) 

11 

FOkMATt lXf H 1.3, IX, lHlf 120A1 1 

GO  TO  1  n 

MA ) 097 ) 0 

15 

WRITE(6,80) 

MAT  09720 

WKITF(6,01)  VAL 

MAT  09730 

80 

FORMAT ( 13X,  IH!  ) 

81 

FORMAT ( lXf F 11.3, lXf 1HI ) 

10 

CONT 1 NUF 

MAI  09760 

WRITF(6, 13) 

MAT  097/0 

13 

Format ( 12X,  io(  10H - 1 > ) 

MAT  09780 

WKITE<6,14>  ( AVAL  (11,1*1.10) 

MAI  09790 

14 

FORMAT l 13X, 1 0 ( 7X,F3. 1 ) ) 

WkITE(6,70)  AVE, VAR/ SB 

MAT  09800 

70 

FORMAT  (  1  Fin,  16X,  '  MEAN  ■  '  ,  F  1  1  ,3,  15X  ,  '  VAR  I  ANCE  ■  '  ,  F  1  l  .  3,  1  bX  ,  '  ST  AND  AR  0 

1  LFEVI  ATIOM  ,  FI  1.3) 

IF<  INDEX. F.Q.l  1  GO  TO  94 

return 

mai 098in 

900 

WKITF-U.,  101  ) 

IF(LIM.NE.O)  WRITE (6,002)  BIG 

IF (LIM.FO.O)  WRITF(6,903) 

MM  09820 

902 

FORMAT  (  IX,  '  ALL  VALUES  IN  THE  ARRAY  ARf  IDENTICAL  ,ANI)  ARE 

MAT09840 

1  F  A  2  .  A  ) 

MAI  09860 

903 

FORMAT ( IX, ' THI S  NODE  WAS  NEVER  A  TERMINAL  NODE’) 

RETURN 

MAI  09860 

END 

SllbPOUl  ILF  TFPM(  I  ) 

MA I 0987  n 
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C  THIS  SUBROUTINE  prints  error  messages 

c 

CuMMnN/AR01/ARC<  100>  ,  INOnrt  100>  ,  OMODE  (  1  OO  >  ,  l TIMET! 100)  ,  TAR61  (  100  )  , 
it»rg2<  if»0) ,  targ3(  100 ) ,  r.nsTC(  ino>  ,costV(  ioo> 

COMMON/IRRZ/IREPIT 

COMMON  / AC C?/  TIME! 100) , ISTAT(IOP), PROP! 100) 

COMMON/ P APA/NNODE/ N ARC 

00  TO  (  J  ,2,3,4,5,6,24,25,27) ,  I 

1  WRITE (0,7) 

7  FORMAT  (  IX,  ,***ERROR***PRDRABJLISTlC  OUTPUT  NUDE  -  WRONG  NUMBER  OF 
) ARCS  S T  A T T  P *  *  *  N 0 0  I N  '  ) 

call  lx  it 

2  WRITE  (N,  HI 

8  Fl'Rl'AT(  IX,  ’***EPROR***PRORABILISTIC  OUTPUT  NODE  -  WRONG  ARC  LISTED 
1  »  * •NOD  I N • I 

CALL  EXIT 
.5  WRITE <6, 9) 

9  FORMAT! IX, ‘ *«*FnROP«**pRnRABILISTIC  OUlPUT  NODE  -  NO  ARC  INITIATED 
l***PROriK'  ) 

Ca|  L  EXIT 

4  WRITE! 6,10) 

10  FORMAT!  IX,  '  ***ERROR***NO  TNpliT  RULES  WERE  SAT  I  SF  1  Eii***Ri)NSYS  '  ) 
WRITE! ft, 13)  I  REP  IT 

13  FORMAT! 1HD, 17X, ' ITERATION  ’,13) 

WRITE! ft, 14) 

14  FORMAT! 1CX, 'STATUS  OF  ALL  ARCS  FOLLOWS',/) 

PC  IP  I* 1 , NARC 

lb  WK I Tt ( ft, 17)  ARC ( I ) , I  ST  AT ( I ) 

17  FORMAT ( I7X, A4,5X, 1 1 ) 

WHITE (ft, 18) 

18  FORMAT  (  10X,  '  WHEREO  0  MOT  I  N  IT!  ATED  '  ,  /  1  9X  ,  '  1  I  N  I  T  IATF.D '  ,  /  1 9X, 

1'2  COMPEETED  SUCCESSFULLY' ,/19X,  '3  FA  I  LED ' , / l 9x , 

2' 4  COST  ♦  TIME  VALUES  ALREADY  CONSIDERED') 

CALL  EXIT 

5  WRITE (ft, 11) 

11  FORMAT  (  IX,  '  •••ERROR***COul.D  NOT  DF.TERMINE  TERMINAL  N0dE**«PTERM  '  > 
CALL  EXIT 

6  WRITU6,  12) 

12  FORMAT  t  IX,  '  ***ERROR***LLK  AND  LL  MUST  bE  E(JU  AL*  *  *  SGR  APH  '  ) 

CALL  EXIT 

24  WR 1 TE ( 6 , 23 ) 

23  FORMAT ! IX, ' ***ERROR***NO  TERMINAL  NOdEs  ARE  PUNCHED  TO  PRINT') 

CALL  EXIT 

25  WR I TE ( ft, 26 ) 

26  FORMAT! IX, ' ***ERROR***CAN  NOT  FIND  TERMINAL  NODES  FOR  SCAN') 

CALL  EXIT 

27  WR  I  TE ( ft, 20 ) 

28  FORMAT! IX, '***ERROR*«*CHECK  ALL  ARC  CARDS, ' /I  IX,  '  AT  LEAST  ONE  SHOW 
IS  FIRST  AMD  THIRD  ARGUMENT  EQUAL  WHILE  TIME  DISTRIBUTION  TYPE  IS 
2TRI  ANGULAR'/l IX, 'CHANGE  TO  CONSTANT') 

CALL  EXIT 

STOP 

END 

BLOCK  DATA 

COMMON/IDD/  RUN  I D ( 20 ) 

COMMON/PARA/  NNODE,NARC 
COMMON/ ARC2/  DUM ( 1 00 ) , I DUM ( 1 00) ,DUMA( J  00) 

COMHON/TERNN/  LLDIJM  (  31  )  ,  ABDl'M!  1000  )  ,  KlDUM!  530  ) 


MAI  07720 
MAI  08720 
MAT08730 
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Cl'Ml'ON/  Al-X  1  /  BpUn  (  10O  )  ,  JnilK(3no  )  ,  CDUM(f)On) 
r.OMMOK/I.DDl  /  KPUMI  21  00  >  ,  7DUH(  1UO0  )  #  LDllM(  HOO)  r  YDuMHOO  > 
Cur.MON/Tf  PM  1  /  K J D U M (  1  i  ) 

DATA  f  III  1  [,/20 *4H  / 

n  A  T  A  “If  (i pi;/ 0/ /  MARC/O/ 

ii  at  a  imn/ionoo.o/,  I  pun/ 1  oo*o/»ouha/ioo*o.o/ 

data  I.Lnlin/31*0/r  ABDIIM/  1  uCO*O.U/,KLDUM/S30*0/ 

DATA  B  pill!/ ion*  0.0/z  JDUM/ 30(1*0/,  CPUM/bO 0*0.0/ 

data  K  mill  /  2 100*0/, ZUUM/10 00*0.0/, LDUM/50 0*0/, YD  ljM/100* 0.0/ 

data  K.mnM/ 1 1*0/ 

EM' 

*  UST(STIlP) 


MA  I  OfJft.TO 


t<  A  I  0 


r  A  I  0  8  A  4  0 

MAT  fi  h  l\  M II 
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APPENDIX  I  I t 

SUBROUTINES  THAT  MAKE  UP  MATHNET  AND  RISCA 
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Next  page  is  blank. 


This  appendix  includes  a  detailed  flow  chart  for  each  of 
the  25  subroutines  in  RiSCA**.  Table  I  contains  a  listing  of  the 
25  subroutines  in  the  sequence  of  their  appearance  in  the  appendix. 
In  addition  to  the  name  and  sequence  number,  a  brief  description  of 
the  subroutine  function  is  given  in  Table  I. 

TABLE  I.  SUBROUTINE  FUNCTION 


Name 


Functions 


I  .  Subrouti ne  iDIN 

2.  Subroutine  ARCIN 


3.  Subroutine  NODIN 

4.  Subroutine  SCAN 

5.  Subroutine  RUNSYS 

6.  Subroutine  REPSET 

7.  Subroutine  ALLFIR(f) 

8.  Subroutine  PROFIR(I) 


Reads  in  run  identifier  card  name 
of  network. 

Reads  in  the  Arc  cards.  Stores 
and  orders  Arc  data  into 
corresponding  arrays. 

Reads  in  the  Node  cards.  Stores 
and  orders  Node  data  into 
correspond i ng  arrays . 

Summarizes  the  input  data  and 
prints  out  an  input  listing. 

Runs  the  simulation  and  prints 
the  results. 

Sets  the  number  of  iteration. 

Fires*  nodes  using  "AND"  output 
arcs.  Fires  all  output  arcs 
from  fix  node. 

Randomly  selects  an  arc  that 
exits  from  a  probability  node. 


*  Initiates 

**Most  of  these  subroutines  are  the  same  as  those  used  in  MATHNET. 
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TABLE  I.  SUBROUTINE  FUNCTION  (Continued) 


Name 

Funct i ons 

9.  Subroutine  ARCCHK 

Checks  all  of  the  arcs  to  see  if 
they  have  been  initiated.  Checks 
the  probability  of  successful 
completion  of  initiated  arcs. 
Calculates  the  time  for  comple¬ 
tion. 

1  -  normal  distribution 

(sub  GAUSS) 

2  -  triangular  distribution 

(sub  TRIANG) 

3  -  uniform  distribution 

(sub  UNIF) 

10.  Subroutine  GAUSS 

Calculates  normal  random  variables 

II.  Subroutine  RANDU 

Generates  uniform  random  numbers. 

12.  Subroutine  TRIANG 

Calculates  the  triangularly 
distributed  random  variables. 

13.  Subroutine  UNIF 

Calculates  the  uniform  random 
var i ab 1 es. 

14.  Subroutine  NODCHK 

Executes  the  output  rules  of 
those  nodes  whose  input  rules 
have  just  been  satisfied  by  the 
arcs  completed  in  subroutine 
ARCCHK.  Determines  what  nodes 
are  ready  to  fire.  Fires  those 
nodes  that  are  ready. 

15.  Subroutine  ANDTST 

Tests  MAND”  nodes.  Tests  the 
time  (t)  taken  by  each  arc  and 
stores  this  information.  Saves 
the  time  (t)  of  the  longest  arc 
to  comp  1 etion. 

16.  Subroutine  ORTST 

Tests  "OR”  nodes.  Calculates 
cumulative  time  needed  to 
satisfy  input  rule.  Stores 
time  (t)  taken  by  each  arc. 
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TABLE  [.  SUBROUTINE  FUNCTION  (Continued) 


Name 

Functions 

i7.  Subroutine  ONEONE 

Tests  and  runs  one-one  nodes  (i/i). 

18.  Subroutine  ONEBAR 

Tests  and  runs  one-one  bar  nodes 

< I/I > . 

19.  Subroutine  PREFER 

Tests  and  runs  preferred  nodes. 

20.  Subroutine  ITALL 

This  subroutine  will  handle  all 
Terminal  Node. 

2i.  Subroutine  ENDIT 

Stores  information  on  completed 
arcs . 

22.  Subroutine  PTERM 

Determines  which  terminal  node 
was  reached  first  in  each 
iteration.  Calculates  how  many 
times  each  terminal  node  was 
selected.  Computes  the  total 
cost  and  completion  time  for  each 

I te ration . 

23.  Subroutine  SGRAPH 

Determines  what  is  to  be  graphed 
and  supplies  titles  for  each 
graph.  The  actual  graphing  is 
done  by  subroutine  GRAPH  except 
for  the  final  graph  which  lists 
the  probabilities  of  reaching 
the  various  terminal  nodes. 

24.  Subroutine  GRAPH 

Does  the  actual  printing  of  the 
time  and  cost  graphs. 

25.  Subroutine  TERM 

Prints  all  error  messages  in 

MATHNET  and  RISC A. 
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Subroutine  ARCIN 
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© — H 


DO  100  \ 

I  =  i ,  nnqde  \ 

NOUM-I  / 


(  RETURN  ) 


ITYPE  (NDUM)SIDUM 
OTYPE  (NDUM)sIODUM 


tNODTiINODT+1 

1NOD1UNOD1>NDUM 

Subroutine  NODIN 
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Subroutine  NODIN  (Continued) 
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L  =  N  ARC 


YES 


NO 


DO  405 
L  =  1,NARC 

DO  406 
K  :|,MM 


YES 


YES  / 

*  \A  RC  ( L) 


OARC  (NDUM,  K) 


I  ARC  (NDUM,  K ) 


NO 


YES 


I A  RC  NDUM,  K)  = 


100 


I  ARC  I  (NDUM)  = 
IARCI (NDUM) + 
1 


READ: 
MM  ,  ( Y I  (I), 
YO  (I),  1=1, 
MM) 


Subroutine  NODIN  (Continued) 
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FIRE  NODE  I 

USING  STOCHAS- 

TIC 

CONSIDERATION 

IXXX=IXX 

CALL  RANDU  (IXXX) 
IXY.RVAL) 

11  =  OARCI  (1) 

IX  X  =  1  X  Y 

AH  IGH  =  0 

RETURN 


YES 


CALL  TERM  (3) 

YES/ 

IF 

\NO> 

/if\ 

/  RVAL  >  \ 
ALDIO& 

xrvahahigh/ 

KK-OARC  (I , K  ) 
ALOW1  AHIGH 

ISTAT  (KK)  =  1 

K:ll 

AHIGH=  ALOW  + 
PP0A  (I,  K) 

NO 


■W  K=K+1 


Subroutine  PROFIR  (l) 


in 


Subroutine  ARCCHK 
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Subroutine  GAUSS  (IX,S,AN,V) 


GENERATES 
UNIFORM 
RANDOM 
DIST. 


c 


J  =  IX  25 
J  =  J-  { J/67108- 
8 64  )X  67108864 


I 


1=11111111 

J  =  JX  2  5 

J  =  J  -  (J/67108- 
864)*  67108864 

1 

f 

RAN  =  FLOAT  (J  ) 

/ 6 7108864* 

J  =  JX  5 

J  =  J  -  ( J/67108  - 
86 4  )X  67108864 

/  * - 

Subroutine  RANDU  (I,J,  RAN) 
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Subroutine  TRIANG  (1XT.  A,  B  ,C,  X) 


Subroutine  UN1F  (IX, OPT,  PES.  VAL) 


© 


SEE  WHAT  NODES 

V-  ■  , 

ARE  READY  TO 

kf  i  yri-n 

/  * 

FIRE  .FIRE  THOSE 
THAT  ARE  READY^ 

MNU-U 

DO  ? 

I  =  1  ,NIMODE  j 


Subroutine  NODCHK 


115 


YES 


^KOTYPE  (I )  =  5 


OTYPE  (I)  =6 


p TYPE  (I)=  7 


I  =  NNODE 


KIND=  KIND  +  J 


YES 


KIND  =  KIND+J 


YES 


KIND  =  KIND+J 


NO 


YES 


RETURN 


Subroutine  NODCHK  (Continued) 


I  16 


Subroutine  ANDTST  (I , J ) 
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Subroutine  ORTST  (I,J) 
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Subroutine  ONEONE  (ILK,  J ) 
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1 18  A 


TIM  =  0.0 


DO  6 
K  =1,11 


u 

OZ 

< 

ii 

(ILK.  K) 

1 

TIMEN  (ILK)  =  TIM 

r  ( 

^  \ 

YES 

Subroutine  ONEBAR  (ILK,J)  (Continued) 
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122 


23 


Subroutine  ITALL  (II) 


Subroutine  ENDIT  (KEY) 


124 


22 


DETERMINES 
COST  OF 
ITERATION 


DO  1 

1  =  1, NODI 


NO 


COSTZ  (IRR)  = 
COSTZ  (IRRJ  +  COS- 
TC  (I)  +  COSTV(I)  * 
(TIME  (I)-TIMEN(LM) 

i 


NCOUNT  (JJ)  = 
NCOUNT  (JJ)  +  1 
LM=  NCOUNT  (JJ) 
IRR  =  IREPIT  +  1 
TIMEZ  (IRR)  =SMTIM 


NODEZ  (IRR)  =INSM 
COSTZ  (IRR)  =  0.0 


COSTZ  (IRR)  =  COSTZ 
(1RR)  +  C0STC(I)  + 
COSTV(I)  X(SMTIM- 
TIMEN  (LM)  ) 


Subroutine  PTERM 
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Subroutine  SGRAPH 
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Subroutine  GRAPH  (ARR,UM) 
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AINS  =.02 

_ J  D°7  \__ 

26-1  +  1 

L  1  1  1  I  I  JlA  1  _  A 

[ 2 1,26  ji 

NUW  (LM|  *  0 

AHIGH  =  0.0 
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Subroutine  TERM  (I) 


In  MATHNET,  the  error  message  refers  to  an  integer  I  and  is 

stated : 

"EXECUTiON  TERMINATED  FOR  REASON  !" 


The  fol lowing 
i  nd i cate : 

is  a  list  of  the  possible  integers  and  the  errors  they 

a. 

(x  is  stated) 

I  =  99.  the  number  of  stochastic  outputs  from  node  x 
is  not  correct  -  (NODIN) 

b. 

1  =  3030.  no  input  rules  were  satisfied  (RUNSYS) 

c. 

1  =  1210.  I/I  output  is  not  correct  (NODCHK) 

d. 

1  =  1211.  I/I  output  is  not  correct  (NODCHK) 

e. 

1  =  1212.  preferred  output  is  not  correct  (NODCHK) 

f.  I  =  100.  probabilistic  output  node  -  wrong  arc 
listed  (NOD  IN ) 

g.  t  =  69.  probabilistic  output  node  -  wrong  number  of 
arcs  stated  (PTERN) 


h. 

(PR0FIR) 

1=1.  probabilistic  output  node  -  no  arc  initiated 

i . 

1000  (SGRAPH), 

1  -  444.(6)  number  of  iterations  set  is  greater  than 

in  RiSCA,  the  errors  which  result  in  termination  of  the 
execution  are  printed  as  statements. 

The  following  is  a  list  of  the  RiSCA  error  statements: 

a.  ***ERR0R***  PROBABILISTIC  OUTPUT  NODE  -  WRONG  NUMBER 
OF  ARCS  STATED***NOD IN 

b.  ***ERR0R***PR0BABI LI  ST  1C  OUTPUT  NODE  -  WRONG  ARC 
LISTED***  NODIN 

c.  ***ERR0R***PR0BAB I L I  ST  I C  OUTPUT  NODE  -  NO  ARC 
I N I T I ATED***PR0F I R 


d. 

***ERR0R***N0  INPUT  RULES  WERE  SAT  1 SF 1 ED***RUNSYS 

e. 

***ERR0R***C0U LD  NOT  DETERMINE  TERMINAL  NODE***PTERM 

f. 

***ERR0R***LLK  AND  LL  MUST  BE  EQUAL***SGRAPH 

132 


g.  ***ERROR***NO  TERMINAL  NODES  ARE  PUNCHED  TO  PRINT 

h.  ***ERROR***CAN  NOT  FIND  TERMINAL  NODES  FOR  SCAN 

i .  ***ERROR***CHECK  ALL  ARC  CARDS  AT  LEAST  ONE  SHOWS 
FIRST  AND  THIRD  ARGUMENT  EQUAL  WHILE  TIME  DISTRIBUTION  TYPE  IS 
TRIANGULAR  CHANGE  TO  CONSTANT. 
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Next  page  is  blank. 


APPENDIX  IV 

MATHNET  OUTPUT  FOR  THE  GOING  TO  WORK  PROBLEM 
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YOU  ARE  IN  ENTER  NOOE  MOOE 

ENTER  NODE  NAME* INPUT  RULEtOUTPUT  RULE 

FORMAT  IS  A4* 1 1  *  1 1 

INPUT  AND  OUTPUT  RULES  ARE  AS  FOLLOWS 

RULE  NUMBER  INPUT  RULE  OUTPUT  RULE 
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GRAPH  OF  COMPLETION  COSTS  FOR  TERMINAL  NODE  N016 
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THE  MEAN  IS  43.09  THE  VARIANCE  IS  1.72  THE  MEPIAN  IS  43.00  THE  MODE  IS  43.36 
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THE  NEW  ORIVING  TO  WORK  PROBLEM 
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THE  MEAN  IS  48,  20  THE  VARIANCE  IS  1.30  THE  ME  C I  AN  IS  47.71  THE  MODE  IS  50.40 
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THE  MEAN  IS  51.81  THE  VARIANCE  IS  5.62  THE  MEDIAN  IS  51.32  THE  MODE  IS  52.15 
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